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The Transient Reaction of an Airfoil Due to 
Change in Angle of Attack at Supersonic Speed’ 


CHIEH-CHIEN CHANGT 
The Johns Hopkins University 


SUMMARY 

Near the destination of a flying guided missile, it is sometimes 
necessary to change the angle of attack of the wing or any control 
surface from one value to another in a short interval of time so 
that the trajectory can be modified to intersect the target. In 
order to calculate such a trajectory accurately, the transient reac- 
tion of the lift, moment, and wave drag of the wing has to be 
evaluated as a function of time. 

For the preliminary investigation, it is assumed that the angle 
of attack of an airfoil moving at constant supersonic speed starts, 
at a certain instant, to increase linearly with time until it reaches 
a specific value, after which it stays constant. For the analysis, 
both the angle of attack and its changing rate are expressed in the 
form of Fourier integrals. 

The linearized differential equation for the unsteady flow of a 
compressible fluid is used along with Possio’s solution for a simple 
source. Unljke the solution of a harmonically oscillating airfoil 
where numerical integrations of finite integrals of the Bessel func- 
tions of the zero and first order are necessary, the lift and moment 
of this case can be explicitly expressed in a closed form as a func- 
tion of time. Some examples are shown. 

As expected, the lift, wave drag, and moment will be reduced to 
Ackeret’s results of steady two-dimensional supersonic flow as 
soon as the transient reaction damps out after the change in angle 
of attack ceases. 

The problem is also approached by finding the unit solution of 
an arbitrary change in angle of attack of constant rate with time. 
By superposition of a proper number of such unit solutions, the 
pressure, lift, and moment can be evaluated for the angle of 
attack with any arbitrary function of time. Some integrals are 
derived in the general case if the second derivative of the angle of 
attack with respect to time is continuous. 


Presented at the Aerodynamics Session, Sixteenth Annual ” 


Meeting, I.A.S., New York, January 26-29, 1948. 

* This investigation was initiated at The Glenn L. Martin 
Company last year. The author wishes to express his deep 
appreciation to Prof. Francis H. Clauser for his discussions and 
suggestions. His thanks are also due to Misses Vivian O’Brien 
and Patricia Clarken for completing the numerical calculations, 
drawings, and topographical work. Acknowledgment should be 
made to Prof. H. Ashley and Mr. G. Zartarian of M.I.T. for 
their advice on some of the misprints in the I.A.S. Preprint No. 
134 of this paper. 

+ Associate Professor of Aeronautics. 


(1) INTRODUCTION 


Shes THEORY OF THE NONSTATIONARY MOTION of an 
airfoil has many important applications in flight, 
such as flutter and wing loading due to gust. There 
have been many investigations in this field. In incom- 
pressible cases, notable work has been accomplished by 
Birnbaum, Glauert, Wagner, Kussner, Theodorsen, 
Jones, Reissener, von Karman, and Sears. In the two- 
dimensional subsonic cases, some of the pioneers are 
Possio, Borbely, and Kussner. As an extension of the 
subsonic cases, Possio,! von Borbely,* Garrick and Rubi- 
now,” Miles,* Temple and Jahn,'* Schwarz,‘ and Hénl"® 
have treated the supersonic cases. Both in the sub- 
sonic and supersonic cases, the theory is based on the 
linearized concept of the thin wing where the local 
downwash is small in comparison to the. flight speed 
and/or sound velocity. Except for the case of vertical 
gust which was treated by Kussner and von Karman 
and Sears® in incompressible cases and Schwarz‘ in 
supersonic cases, most of the theoretical investigations 
have treated the case of steady harmonic oscillations 
with respect to time in a number of degrees of freedom. 

There is another application of the theory to the 
calculation of the transient reaction on an airfoil due to 
change in angle of attack or to deflection of the control 
surface attached to it. Although not treated much in 
incompressible or subsonic cases, the application cer- 
tainly becomes increasingly important in the supersonic 
case. 

The present treatment investigates the transient re- 
action of a two-dimensional wing in supersonic flight 
when the angle of attack of the wing changes from one 
fixed value to another according to a given function of 
time. This is based upon certain current needs in the 
design of supersonic guided missiles and pilotless air- 
craft. 
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_ In the performance of a guided missile flying near its 
moving or stationary target, its trajectory must be 
precisely predetermined as a function of time in order to 
reach its destination accurately. Sometimes, however, 
after the missile is fired, certain unexpected factors may 
spoil the scheduled trajectory, and some correction of 
the path of the missile may be found necessary, par- 
ticularly during flight near the target. One feasible 
means of performing such a correction is to deflect the 
wing or any other control surface to another appropriate 
angle of attack during the course of time in such a 
manner that the missile will not miss the target. Ina 
certain sense, this may be considered as an analogy of 
von Karman and Sears’ case of a sharp-edged gust in 
incompressible flow, although the methods of approach 
in the two cases are entirely different. 


In the linearized two-dimensional theory of super- 
sonic nonsteady flow, the vortex theory can no longer be 
applied because the flow in the wake cannot influence 
the pressure on the wing. However, Possio’s! solution 
of acoustic sources can be used with interesting results. 
In the two-dimensional case the strength of the acoustic 
source at any point on the surface is linearly dependent 
on the downwash at that point and at that instant and is 
independent of the downwash of the neighboring points. 
Furthermore, as shown by von K§4rman," to the 
approximation of the linearized theory, the effect of the 
vertical displacement from the neutral plane on the 
downwash at the solid boundary may be neglected if the 
vertical displacement is small in comparison to the air- 
foil chord. Consequently, the downwash can be con- 
sidered to occur in the plane of the undeflected thin 
wing where the corresponding source strength is dis- 
tributed to replace the airfoil. Besides, the disturbance 
due to an acoustic source lies entirely in the downstream 
Mach cone and cannot be beyond it. These facts make 
the linearized theory simpler in the supersonic than in 
the subsonic case. 


The present treatment is based on the superposition 
of Possio’s solution of simple sources. The angle of 
attack is first constant, then changes linearly with time 
for a short time interval ¢, and then remains constant. 
A Fourier integral is used for the downwash to evaluate 
the velocity potential. In contrast to the case of har- 
monic oscillation, in which numerical integrations of 
finite integrals of Bessel functions are involved, the 
pressure, lift, and moment can be expressed explicitly as 
a function of time. 


If the angle of attack is an arbitrary function of time, 
it is shown that the pressure, lift, and moment can be 
obtained approximately by superimposing a number of 
the downwashes that linearly increase with time. If 
the second derivative of the angle of attack is a con- 
tinuous function of time, the lift and moment can be 
expressed as a group of integrals. 


(2) NOMENCLATURE 











¢ = wing chord 
Cp = transient wave drag coefficient 
Cyr = transient lift coefficient 
Cue = two-dimensional lift coefficient in steady case 
Cu = transient moment coefficient about xo 
Cua = two-dimensional moment coefficient in steady case 
Cp = 2(p — p:)/p:U?, pressure coefficient 
Co == 2a,/(M? — 1)1/2 
¢ = sound velocity in free stream 
F( ) = function concerning the lift coefficient 
f( ) = function concerning the pressure coefficient 
G( ) = function concerning the moment coefficient 
Jo( ) = Bessel function of the zero order 
M = = Mach Number 
m = (MJ? — 1)/M, parameter of Mach Number 
fi = Static pressure in free stream 
, =e Me 8? OF - Dig — 9? + & - 94) 
(M? — 1) 
Pe A(x — 8)? — (MP 1)(y — 9)?/(M? - 1) 
t = time 
: = t/t, nondimensional time 
i = C/U 
‘ = (x — xo)/U 
te = (§ — x0)/U 
w = z-component of velocity 
U = free-stream velocity 
u = x-component of velocity 
v = y-component of velocity 
2 = axis along chord direction 
x’ = x/C per cent of chord 
Xo = axis of rotation of the entire airfoil 
xe’ = = %0/C 
y = vertical axis 
z = spanwise axis 
a = angle of attack 
a = total increment of angle of attack from f to t; 
c = source location along z-axis 
n = source location along y-axis 
r = (x — £)/c(M? — 1) 
p’ = fj’ — 7’ 
% = x/U(t—) = x'/(t' — &’) 
Vio = x0/U(t — ti) = xo'/(t’ — ti’) 
vie = C/U(t — &) = 1/(t’ — k’) 
é = source location along x-axis 
pi = free-stream density 
T = time interval 
r’ = 1r/t; nondimensional time 
6 = $/[—a/(M* — 1)'”] 
o = velocity potential 
po = velocity potential of a simple source line 
w = frequency parameter as an integration variable 


(3) STATEMENT OF THE PROBLEM 


The present problem is to find, at any time, the force 
and moment acting on a thin airfoil in a supersonic 
stream of constant velocity, when the angle of attack 
changes as follows: (a) After being constant for a long 
period of time, it increases linearly for a short time in- 
terval and then remains constant; (b) it changes 
arbitrarily with time. 


The assumptions are: (1) The flow is isentropic and 
consequently irrotational; (2) the fluid is nonviscous; 
and (3) the angle of attack and its rate of change are 
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sufficiently small that the downwash is small in com- 
parison to the sound velocity. 


(4) LINEARIZED DIFFERENTIAL EQUATIONS AND 
ELEMENTARY SOLUTIONS 


Under the above assumptions 1 and 2, the differential 
equation for unsteady three-dimensional flow can be 
obtained and is nonlinear. Little is known about its 
solution for any particular case. However, if the 
disturbance velocities are sufficiently small compared to 
the sound velocity and to the difference between the 
free velocity and the sound velocity that their squares 
can be neglected, the differential equation can be 
linearized as 


; . 2M 1 
Ow + bz = (M?* — bez + oes gon + a du (4.1) 
c 2 


where 1f = U/c, c is the sound velocity of the free 
stream, and ¢ is the velocity potential such that V¢@ = 
(U+ uji+ vj + wk, where i, jh R are unit vectors in x, 
y, and z directions, respectively. Through this process, 
three advantages are gained. First, the particular solu- 
tion of a simple source or of a simple doublet is rather 
easy to derive by means of the classical mathematical 
analysis. Second, being linear, its solutions are super- 
imposable. Third, the general solution consists of one 
or both of two parts, odd and even. The even part of 
the solution can be uniquely determined by superposi- 
tion of simple sources and the odd part by superposition 
of simple doublets. Thus, we can superimpose the 
particular solutions to fit particular required boundary 
conditions. Then, the whole problem is reduced to 
determining the strength of the source and/or of the 
doublet with the given boundary conditions. 


To find such a simple solution, we consider the effect 
at a point P caused by an acoustic source at Q in the 
presence of a supersonic flow with uniform velocity U. 
The points P and Q have coordinates x, y, z and &, 9, £ 
which are fixed in space. Fig. (la) shows the view on a 
plane which passes through P and Q and which is 
parallel to the supersonic stream. As a sound pulse 
issues from the source, the spherical waves moving with 
the fluid will increase in diameter with acoustic velocity 
c while being transported downstream with velocity U. 
The envelope of such a system of spherical waves is the 
stationary Mach cone. If P is outside of the Mach 
cone—i.e., in the so-called forbidden region—the 
velocity potential at P must be zero. However, if P is 
inside, there are two spherical waves simultaneously 
passing P at time ¢. Both issued from Q at earlier 
times. These time intervals can be determined from 
the cosine law. 
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TRANSIENT REACTION 








OF AN AIRFOIL 











(c) For t<% 





Fic. 1. A point P and an acoustic source Q are both in a uniform 
flow of supersonic velocity U. The strength of Q is zero when t < 
0. The acoustic waves which are issued by Q and blown down- 
stream reach P simultaneously at t as shown 

(a) P and Q are on parallel streamlines. 

(b) P and Q are on the same line. 

(c) P and Q are on the same line. 


where 





1 —_ 
r= ——— V(x — 8)? — (M? — 1)[(y—2)? + (2-9) 


M*—1 
(4.3) 
It has been shown by Garrick and Rubinow’ that the 


velocity potential at P due to a source at Q whose 
strength, A, varies with time is 


A(é,n, f,f ie 72”) + A(é, n, ct ma 71") 


r 





do’ (x, y, 2, t) = 
(4.4) 


This can be extended to give the potential for a two- 
dimensional line source by summing up the potential of 
sources uniformly distributed along a line = constant, 
n = constant—i.e., by integrating Eq. (4.4) with respect 
to ¢ holding £ and 7 fixed. ‘ 
ott godt 

In this integration, 71”, 72”, and 7 are, of course, func- 
tions of ¢. Points on the line source whose Mach 
cones do not contain P cannot affect the potential at P. 


Consequently, the integration’ may be carried out 
along the source line only between the two points: 


do(x, y; Zz, t) 7" 
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Fic. 2. Distribution of a(t) and a&(t) vs. ¢ for any point on a 
straight airfoil. 
f= 2—- V(x — 82/(M? - 1) - 
and 





f= 2+ W(x — 8)2/(M? — 1)] — »? 


If we use these two limits and assume that at a given 
time the pulses emanating from the line source are 
everywhere of the same strength, then this integral be- 
comes Possio’s solution, which is 














do(x, y; t) _ 
"? A( ’ at a ) , 
J oo = dr (4.5) 
V8? — (x — § — Ur)? — (y — 9)? 
where 
M(x — &) fr’ M(x — &) 
See oP oe, eR —— Ae — 4. 6 
Ee fg > ¢(M? — 1) - (4.6) 
where 
hae M2 — VG 2— (M? — 1)(y — »)? 


7, and 72 are now the shorter and longer time intervals 
for the cylindrical waves to travel from the line source 
at Q(é, 7) to the point P(x, y). [See Fig. (1a).] 

We see that in the two-dimensional case, the spherical 
waves of an infinitely long acoustic line source scatter 
over the whole of the Mach wedge of the source. 
Therefore, P is affected not only by the pulses emitted 
at previous times 7; and 72 but also by all pulses be- 
tween these times. 


1948 


With the principle of superposition, the total poten. 
tial at P due to an infinite source sheet on the positive x. 
axis can be obtained by integrating the potential of the 
above source lines with respect to é 


= So" godt (4.7) 


where & = x — y VM? — 1 is the nearest element of 
the source sheet whose disturbances will reach the point 
P. 

To facilitate the integration, we can substitute c? x 
(M? — 1)(r — 71)(t2 — 7) for c?7? — (x — — — Ut)? - 
y*. Expressing Eq. (4.7) in another form we have 


G(x, y, t) = 


o(x, y, t) 


onda -[ A(t, +0,¢ — 7) 
Wir a a a Ti V(r = 71) (T2 —_ T) 


This expression has been given by Possio! and 
Garrick and Rubinow.* Now the problem is how to 
determine the source strength from the boundary condi- 
tions. For the two-dimensional case, it is known? that 


A(x, +0, t) = *(c/r)v(x, +0, 2) (4.9) 








dr dé (4.8) 


where v is the vertical velocity immediately above or 
below the source sheet. Thus, the negative sign of » 
is to be used as y approaches zero from the positive side, 
and the positive sign of v is to be used as y approaches 
zero from the negative side. It is interesting to note 
that, if the lower wing surface is equivalent to a positive 
source sheet, the upper surface corresponds to a nega- 
tive source sheet of the same strength. Each side can 
be considered independently as a source sheet for the 
present problem, because the influences of the two will 
not meet except in the region behind the trailing edge 
with which we are not concerned in this problem. If 
the flow behind the wing were considered, the problem 
would have to be treated as a doublet sheet of finite 
chord length. 

Substituting Eq. (4.9) into Eq. (4.8), we have the 
velocity potential on the upper half plane, which is 


o(x, y; t) - 











vat — oe ls 
v Tee pea gia 
(4.10) 


If the point considered is on the airfoil but below’the 
surface, Eq. (4.10) becomes 











(x, —0, t) = 
1 1 * £9 v(g, t — 
er eee < z dr df 
we (M*—1)" Jo Jn V(r — 11)(12 — 1) 
(4.11) 
where 
x —H 1 ) x —é 
(2S idliee oo ee _ ps 
Cc M+ 1 U+e 
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(5) THE DowNnwasH DISTRIBUTION 


Assume that the angle of attack and its rate of change 
are both zero at t < 0 and that the rate of change of the 
angle of attack remains constant from time +0 to 4. 
Beyond this time interval ¢,, the angle of attack is kept 
at a constant value a;. Figs. (2a) and (2b) show the 
distribution of the angle of attack a(t) and its derivative 
a(t) at a fixed point (x, +0). With Dirichlet’s integral 
itcan be shown that 


a [s wt — sin w(t — t) i 
T 0 Ww 
ay © dw 


2irJ/Jo w 





a(t) = 


a gH) 


Il 


[e(1 ad cr _ oH 


(5.1) 


where &(t) = da(t)/dt and & = a:/t;. It can be shown 
easily that this equation satisfies the proposed physical 
condition at any instant of time, positive or negative. 
Similarly, we can express a(t) as 


yt i * dw 
ot a = = (cos w|t| — cos w|t — 4]) 
2 0 


T w? 


~ fe“ (1—e-*") te (1-eT ™")] 


(5.2) 


where df; is introduced for a. In order to represent 
a(t) as specified, the absolute values of both ¢ and t — 4 
are introduced. 

No loss in generality is occasioned by taking the angle 
of attack before deflection to be zero. Because of the 
linear character of lift and moment, any steady state 
contribution can be added to the transient values com- 
puted here. The effect on the wave drag will be con- 
sidered later. In case of constant dp (i.e., da att < 0) 
and & ~ 0, the above equations are still true if a is 
replaced by &: — a, & by Ad, and aby Aa. In other 
words, only the increment in & at ¢ = 0 is considered. 

From Fig. 3, we see that, within the limits of our 
approximations, the vertical velocity of any point of a 
flat plate airfoil can be expressed in terms of a and & as 
follows: 
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Fic. 3. Geometry of the airfoil at time t (0 < tf < 4). 








Fic. 4.- Downwash distribution (—v/U) in the t-x plane. Note 
the discontinuities of downwash at ¢ and 4. 
v(x, t) = —{U[a(t)] + (x — x)[a(é)]} (5.3a) 


where xo is the location of the hinge axis about which the 
wing rotates. 
If we put (x — x) 


v(x, 0, 2)/U = 


Introducing Eqs. (5.1) and (5.2) into Eq. (5.3b), we 
can obtain 


SS O, t) Er cat 
ay U 
| © dw ° ° 
cos w|t|) + t, — [sin¢ — sin (¢ — 4)] 
0 w 
_ Th 4 Lf Shel (*- i) ” 
2 2/70 w? 1 

My] ae il (4 i)a _ omy | 
1 ) 


(5.4) 


o/U = t,, then, ‘approximately, 


—a(x, Ye, t) — te[a(x, ye t)] (5.3b) 


i} 
= (cos w|t — 4 -- 


Fig. 4 shows the proposed downwash in the ¢-x plane. 





—0,t) = 





o(x, 








The velocity potential will be used to determine the pressure distribution on the airfoil. 
the and x derivatives at y = Oare necessary. Since y can be placed equal to zero before or after differentiating ¢ 
and x, we simplify the work by using ¢ evaluated at the airfoil. 


a(M? — toa ££ V(r — 71) = — 7) 


(6) CALCULATION OF THE VELOCITY POTENTIAL 


For this purpose, only 


From Eq. ats we have 
v(&, t — 





dé dr 
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Substituting ¢ — 7 for t and é for x in Eq. (5.4), we have 


mw v(é, t — 7) 


ay U 2 Z 0 w? 1 


where ¢, = ( — xo)/U. Then 
. “= OR i Ud J 
o(x = 0, t) = «(M2 — on. 1) 1\'/2 ‘Se = {a se = 
ne *drdé 


roe ae ee -[ i i a 


Introducing a new variable 6, which is defined by 








27 = (t2 — 71) 


we can evaluate® 


So’ SI, dédr V(r _ 11)(t: — 1) = 


and 


et dr 








ape —— e Tiel +71)/2) ft —1:)/2]} cos @ id= 
V(r — 71) T) 0 


where \ = (r2 — 7)/2 = (x — £)/c(.W? — 1) and then 
(T2 + T)) 2? 
Substituting Eqs. (6.4) and (6.5) into Eq. (6.2) 


gi _1 
p(x lan py” 


= (x — 


U. = 


1)a = 





V (r — 71) (T2 = T) 


lf a fo “I ates = i)ci — e~ iat 
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ey os 


‘f “ae (8 +1)( ~ eye | (6.1 

2J0, @* | 
‘- - a — % 

=+ f° f ee patch = (1 ~¢ i) x 


et"drdt 


—twif — 4. th —_ xy ie | 6.2 
a Th V(r —_ 71) (T2 = 7) " 


dédr 











cos 6 + (72 + 71) (6.3) 


(6.4) 


So’ So” dedé = xx 


iwr!) (a 
” ba a Jo(\w) (6.0) 


£)M/c(M? — 1) = Mr (6.6 


—iw(i— Mr 
“a i Thr) x 


(+ 1)a et) fo hw) (6.7 
1 


For convenience of future integration, we can express the above equation in terms of trignometrical functions and 


change tod. Then 


: —] bg ~~ x/e(J oe —1) 
o(x, —0, (asa = nas Ucl. wal a’ {> (cos wit— tt — Mn — coswit — MN) 


Se 


Jo(Aw) + 


where ¢; = (€ — xo)/U = t, — mx, and m = (AM? — 1)/M. 


in Appendix A where the velocity potential is discontin 


In order to understand those lines of discontinuities 
some explanation of the physical relation of the distribu- 
tion of the source and the velocity potential ¢ should be 
treated before the analytic expressions of ¢ are given. 

(a) Zone of Influence.—The acoustic disturbance due 
to the source Q (which is proportional to downwash) at 
any point x on the airfoil at time ¢ influences the pres- 
sure or the velocity potential at other points on the air- 
foil behind x at a time later than ¢. This can be ex- 
pressed in the ¢-x plane as shown in the shaded down- 
stream wedge in Fig. 5a. This wedge vertex is at the 


x/c(M*~1) 
(t, — mdr)dr [? 
Tv 0 0 w 


— > [sin w(t — Wr) — sin w(t — hh — MA)]Jo(Aw) = (6.8 


The integration of the above equation is carried out 
uous along certain lines in the domain of the ¢-x plane. 





source Q(x, /), and the slopes of the wedge sides are 
(U —c)and(U+ c). Outside of this wedge, the dis- 
turbance at Q no longer influences the pressure. This 
wedge area may be called the zone of influence of Q in 
the t-x plane, and it is similar to the downstream Mach 
cone in the steady supersonic flow. 

(b) Zone of Action.—Conversely, the potential at any 
point x on the airfoil at time ¢ is obtained by a surface 
integral [Eq. (4.11)], which sums up the sources that 
have dccurred at a point forward of x and at a time 


previous to ¢. These sources are contained in the 
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graded upstream wedge in Fig. 5b. This wedge is at 
the point P(t, x) and the slopes of the two sides are also 


U—c)and(U+ c). Outside of this area, the acous- 
tie source produces no action on the potential at P. 
Therefore, this area may be called the zone of action 
upon P, and it is similar to the upstream Mach cone in 
the steady supersonic flow. 

Now for the present problem, the downwash is zero 
iorx < Oand¢t < 0, and, consequently, the downwash 
distribution in Fig. 4 is discontinuous along x = 0 and 
‘= 0. Att = 0, an abrupt increment in 4 (i.e., a1) is 
made, and & remains constant hereafter. We want to 
know the effect of this change on the velocity potential 
of a point x on the airfoil at different times. First, if 
)<t<x/(U + oc), the zone of action upon P; is shown 
by the shaded area in Fig. 5c. Later, ifx/(U+¢) Z 
1Zx/(U — c), the zone of action upon P, is shown by 
the shaded area in the figure. Furthermore, if ¢t > 
/(U — c), the zone of action upon P; is shown by the 
shaded triangular area. The potential at x, being a 
surface integral of the shaded area, must be different in 
nature for differently shaped areas. This may be con- 
sidered as one reason for the discontinuities of the 
velocity potential in the ¢-x plane. 

Furthermore, the strength of the source at ¢ = 0 and 
t = Ois discontinuous both int andx. The zone of in- 
fuence of this point source is a wedge with side lines 
t=(U —c)xandt = (U+c)x. The two lines divide 
the region bound by x > 0 and ¢ > 0 into three zones, 
called I-1, I-2, and I-3 as shown in Fig. 6a. Here | 
refers to the first region produced by an increment in @ 
ie., a1) att = 0. In zone I-3, the leading edge exerts 
no influence, and, consequently, the potential at P3 is the 
same as though an airfoil extending to infinity ahead of x 
had started abruptly to rotate. In zone I-2, the in- 
fuence of both the abrupt start in the rotation and the 
eflect of the leading edge is present. In zone I-1, the 
potential is the same as though a finite airfoil were 
mdergoing a steady rotation. The leading edge effect 
is present but not the abrupt start in rotation. Thus 
each zone, the source strength is different in nature 
and the potential is different in expression accordingly. 

The above physical picture occurs because there is an 
abrupt change in & at ¢ = O (i.e., Ad = a), and & re- 
mains constant for ¢ > 0, as shown by the dotted line of 
Fig. 2b. Now at? = ¢, another change in & (i.e., Ad = 
 — & = —d, & being zero) is made as shown by the 
slid line from the dotted line. The second change in & 
tan be treated like the first one. Thus region II, 
where ¢ — 4; > 0, can be introduced. Three zones II-1, 
l]-2, and II-3 can be similarly defined. 

Introducing the new parameters vp = x/U(t — to) (to 
bing 0) and », = x/U(t — t), we can show that, in 
gion I, zone I-1 corresponds to0 Z » Z 1 — (1/M); 
wne I-2 corresponds to 1 — (1/M) ZmZ1+(1/M); 
and zone I-3 corresponds to 1 + (1/M) Z w. 

Similarly, in region II, zone II-1 corresponds to 0 Z 
nZ1— (1/M); zone II-2 corresponds to 1 — (1/M) 
£m Z 1+ (1/M); and zone II-3 corresponds to 1 + 
(/M) Zn. 
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Fic. 5. (a) Zone of influence of source Q. (b) Zone of action 
upon P. (c) Zones of action upon P;, Pe, and P; due to a —Jay 
att=0. Each zone of action is different in shape from the other. 
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TABLE 1 
¢ VS. vo and »; 


























~ Range of » a 
Range of ~~ POST Bin Rene” 
4" Bina | 
Region I - (I-1) (I-2) (I-3) 
ce t—h <0 ¢ = on & = dor ¢ = dos 
(I-1) + (11-1) a 
1 > = oa + ou t 
OfnZ21- M b Useen Impossible Impossible 
— (M2? — 1) a 
Regions 1 an 
I+II t—-—2Zy,Z1+— (I-1) + (II-2) (I-2) + (II-2) Impossible 
M M ? = gn + on > = go + diz 
1 (I-1) + (11-3) (1-2) + (11-3) * GC) + OES) 
1+ a ZyZo a des bes @ = do + on > = go3 + $13 
= —a,(M? — 1)/ tex — %) 








Zones in region I are shown by solid lines and zones 
in region II are shown by dotted lines in Fig. 6b. 

Using the above concepts, the integrations for the 
velocity potential have been carried out in Appendix A. 
It is clear that, once the velocity potential of the first 
region is obtained, the contribution due to the second 
region can be written down immediately. In region I, 
do1, Goo, and gos are introduced to represent the velocity 
potentials in zones I-1, I-2, and I-3, respectively. 
From Appendix A, if Eqs. (A10), (A11), and (A13) are 
substituted into Eq. (A3), the velocity potential in each 
zone can be written as follows: 

(a) Zone I-1,0 Z » Z1 — (1/M). 








—4 Uxt; — 2d 
Ai cco oot + ee ) ne 
"OP — a a 
a alt 2. l a (6.9) 
(M2 —1)"L » «+ 2(M?—1) 


where Eq. (A14) is introduced to reduce jy to $1, and 
where vo = xo/U(t — to) (xo being arbitrary). 
(b) ZoneI-2,1 — (1/M) Zu 21+ (1/M). 


— dy U. xt, 
22 OS + Dy + ® ) 
oo: ava D7 = 2 ll 


— ax? 1 — Voo 
accall ® = 17 ee 
(M aT, Sun u + Su) (M? — =a(| % 


zB 1 -(" a M)+ ‘|- (M? — 1) 
2(M? — =I! = ne 2 M x 
vo — vof 1. _, oa -3]- a =) 
{ = FE sin M(1 — ») — 5 red ata)x 
Vi — M1 — 7) (6.10) 


(c) Zone I-3, v = 1 + (1/M). 











b= ation ( E+ Om — bu + Ou) 


— dy 





— &x*/ vo — voo 
Qin = 
1) 03 M ( ya? ) (6.11) 


0 


Oe = 





It is easy to show that, for t < 0 (i.e., » < 0), the 
velocity potential is zero from Eq. (A3). This checks 
with the given condition. 

Similarly in region II, the contribution to the velocity 
potential due to the second change in @ (i.e., Ad = 
dj — & = —d, d& being zero) can be shown to be 
diy = — oj, Where v and vo in oy are replaced by 1 and 
Vio in oi; (j = 1, 2rd). Here, vio = xo/ U(t — th). 

With the above relations and Appendix A, a table can 
be formulated to represent ¢ for all cases in Eq. (A-3) 
(see Table 1). 

An important feature of Table 1 is that in the over- 
lapping regions ¢ = ¢ox + ¢im corresponds exactly to 
the overlapping zones (I-k) + (II-m), where k, m = 1, 
2, 3. Since go, is the velocity potential in zone I-k 
due to change in & at ¢ = 0 and ¢1,», in zone II-m due to 
change in @ at ¢ = t, we see that ¢ can be obtained by 
superimposing the solutions $0, and dim for t > fy. 

In the overlapping zone (I-1) + (II-1), a becomes 
constant at aj, and all of the transient effect damps out. 
The velocity potential becomes identical to Ackeret’s 
solution of the steady case 

It is apparent that the effect of the first change in 4 
and a at ¢ = 0 cannot influence the pressure fort <.0; 
therefore ¢) = O fort <0. Similarly ¢, = 0 fort <4 
in the second change at ¢ = 4. 

When the zone II-3 overlaps the zone I-3, the poten- 
tial is again independent of time. This means that the 
influence of & discontinuity (Ad = a) at ¢ = 0 just 
counterbalances that at ¢ = t(Ad@ = —a) and this 
potential must be equivalent to that of an infinite chord 
airfoil at constant a. 


(7) PRESSURE DISTRIBUTION 


To the approximation of the linearized theory, the 
pressure coefficient C, is connected to the derivative of 
the velocity potential ¢ as follows: 

Cie Rel at AE wi ee ee 

> pU?/2  Ude Uta Ude 

where ¢ is given in Table 1 for various regions and 
For region I (¢ > 0) the pressure coefficient in 
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TABLE 2 
Cp Vs. vo and 





Range of v9 
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I+ I1| 1 i 1 (I-1) + (11-2 2 Y ea 
\— M Enil+ M Cp = Con + Core | Cp = Cpor + Coiz | aaccahitinall are 
. 3 3 
eS a (I-1) + (11-3) (1-2) + (11-3) (1-3) + (11-3) ~ 
| * oad | Cp = Con + Cos Cp = Cooz + Chis Cp = Coos + Cos = ai 
SS ee ee ne Ricthicannieeniaa a EEE i — 
each zone can be defined as Criy = —(2/U)(doj;/dt) 
Cyoy = —(2/U)(doo;/dt) (j = 1,2,3) (7.2) = (a2 — &)(x/U)f,(r1, 10) (j = 1, 2,3) (7.6) 
Cyoy = (1 — a) (x/U) f;(v0, voo) where f;(1, vio) are the same as f;(v0, vo) given in Eqs. 
a , 7.3), (7.4), and (7.5) except vo, voo are replaced by 
where & — dp is introduced to replace a in order to ( —- 7.5) re Oy. E oy" 


extend the result to a more general case, a& being zero. 
Here f; can be obtained by differentiating Eqs. (6.9), 
(6.10), and (6.11), 


= 2 ddo, 











ie ae 
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filo, Yo) = 





= 2(? an. © 1), tn’ P4= 
M M 


Vo 
(7.5) 
Similarly, in the region II, (¢ > 4), we can define the 
pressure coefficient as 


and V10. 
In Eqs. (7.3) and (7.4), it is easy to verify that, at 
v; = 1 — (1/M) and at the same time, 
Airy vo) = folvy vo) (¢ = 0, 1) 
In Eqs. (7.4) and (7.5) we can get at », = 1 + (1/M) 
and at the same time, 
Sro(vi, v0) ws fs v0) 
Therefore along the two boundary lines of the zones 
of the ¢-x plane [ie., », = 1 — (1/M) andy, = 1 + 
(1/M)] the C,’s are continuous in value, but are dis- 


continuous in the derivatives. This is true for any one 
of a number of regions (i.e., ith region) as shown later. 


(7.7a) 


(¢ = 0, 1) (7.7b) 


Therefore the pressure coefficient C,(t, x) can be 
evaluated if the point is located in one of the three 
zones in region I where 0 < ¢ as shown in the first row 
in Table 2. Of course, only the increment in & at / = 0 
is considered. Nowif the point (¢, x) lies in both regions 
I and II, where ¢ > 4, the increments in & both at ¢ = 0 
and ¢ = & will influence the pressure coefficient, and 
both should be considered accordingly. This is clearly 
shown in Table 2 when ¢ 2 4. 


For ease of treatment, a new time unit ? = C/U is 
introduced. 7? is the time required for the free-stream 
flow to travel a chord length. A nondimensional time 
t’ = t/t can be used for the time scale. With the new 
time unit, Eqs. (7.2) and (7.6) can be written ds 


Coig = (Qet1 — Gey) x'f (v4, v0) («= 0,1) (7.8) 


where x’ = x/C, & is & immediately ahead of ¢,’, a+: is 
& immediately after ¢,’, and a;+1 — a, is the increment 
of &att,’. In Table 2, this new time unit is used. 


In the table » < 0, (¢ < 0), C, = 0. In the over- 
lapped zone (I-1) + (II-1) where the transient reaction 
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(b) 
Fic. 7. Three zones formed in region i + 1(¢> #;) due to &+1 — 
a; at t;’. 


damps out, C, = 2a;/(M? — 1)“ checks with Ackeret’s 
result. It is interesting that C, = 2a,/M in the over- 
lapped zone (I-3) + (II-3). Note the superposition of 
the parts of C, contributed by region I and region IT. 


(8) Lirr, MOMENT, AND WAVE DRAG 
It is generally known that the lift coefficient of an air- 
foil with a chord C is 


* 
C.() = 1 f dx[C,(—0) — G)(+0)] 
C Jo 


Cc 
= 2, [aee, 
CJ/o 


where C,(—0) = C,(x, —0, ¢) and C,(+0) = C,(x, +0, 
t). They are equal in magnitude but opposite in sign. 
C, = C,(x, —0, #) is introduced in the latter equation. 
The moment coefficient is 


Oe NF 
Cu(t) = 2f (x pa x)dxC, 


(8.1a) 


9 £¢ 
° a sf x dx C,—%9'Cy (8.2a) 
The wave drag coefficient is 
Colt) = (2/C) o° de Coa (8.3) 
where a = a(x,?). Fora flat plate 
Cr(t) = aC, (8.3a) 


Therefore, as soon as Cz; is known, it becomes a simple 
matter to calculate Cp, and, consequently, Cp will not 
be treated in this paper. However, in the case of 
a = a ~ 0, Cy + 2ay/(M? — 1)? and a + a, should 
replace C, and a, respectively, in Eq. (8.3). 


(A) Lift and Moment for the Present Case 


Now we may study C, and Cy in region I due to & - 
a att’ = 0. As C, is conveniéntly expressed in terms 
of v= x’ /(t! —_ to’) and vo = xo’ /(t’ —_ to’), we may 
transform Eqs. (8.1) and (8.2) into the forms 


9 OC 1/(t’ — to’) 
Cult’) =‘ f dC, = 2(t’ on a) fi Cyd 
Voc 0 0 


(S. 1b) 
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(c) 
(a) Arbitrary a(t’) vs. ¢’ approximated by a series of 
tangent line elements. (The curve is omitted.) 

(b) All the overlapped regions ahead of t’ are divided into three 
groups: first group, 0 to te’; second group, te+1’ to tm’; and 
third group, tm+1' to’ts’. 

(c) Corresponding & vs. t’ curve (dotted) approximated by 
steps. 
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4 Fre. 10a. Fj’ vs. t — t;’at M = 1.5,2,3,4. Left solid lines— 
| F;’ vs. i’ — t;'(t — t:’ < M/(M + 1)]. Center dotted lines— 
| F,' vs. t’ — t;'{M/(M +1) st’ —t;' S M/(M —1)]. Right 

a F solid lines—F,’ vs. t’ — t;’{M7/(M — 1) St’ — ti’). 
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Fic.9, (a) avs.t’ (abeingcontinuous). (b) Three character- 4 2 — } 
istic integration intervals shown: first interval fram 0 to &’; r 
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ae Se | 
ee a cs 
where vg = 1/(t’ — bh’). ; 4{ 
ap et . 
Se aaa 0 rY bepaa 3 
Cult’) = —— a volo, — Xo Cr 4 
(Yoo) “70 t-t,—- 


1/(t’—to’) Fic..10b. AFj’/xe’ vs. t’ — 4’. Left solid line—AF;’/x9’ vs. 
= 2(t' — tp’)? . 1 ned 2b) i’ — t,;’. Center dotted line— AF,’/xe’ vs. t’ — ti’. Right solid 
( 0) a vodlvoCy *o' Cz, (8.2b) line— AF,’ /x,’ vs. t’ — t;’. 
where Cz, and Cy; are lift and moment coefficients in the 

region I due to & — dat fy’. are three possibilities: (a) t’ — to’ 2 M/(M — 1); 
As shown before, C, is discontinuous along two lines, (b) M@/(M + 1) Zt’ — &’ Z M/(M — 1), where one 
yy = 1 — (1/M) and » = 1+ (1/M) in the #’-x’ plane. discontinuity » = 1 — (1/M) will be encountered; and 
The above integrals have to be split up once ifv >1— (c)0OZt'— ty’ Z M/(M + 1) where two discontinuities 

(1/M) and twice if » > 1+ (1/M). Therefore, there will be encountered. 





(a) t’ — to’ = MW/(M — 1). From Eggs. (7.4) and (8.1) and Table 2, ¢ 
Cru(t’) = 2(t’ — be’) fo"? (da — Go) x’fi(vo)dvo 
Substituting vo for x’, we have 
Cun = 2(t’ — to’)2(dr — Go) Jo” vodvofi(vo) = (a1 — do) Flt’? — bo’) (8.4) 


3 of where Cz, means the contribution of C, due to é; — & and Cz, refers to Cz; when t’ — fo’ 2 M/(M — 1). With 
Eq. (7.4) we can show 


and —? > et ee 4 (7 - es ‘la F’— AF’ (85 
F(t — to’) = 2(t — bo ) 3 vodvo f1(v0) = ata | — bh a (M2 — 1) — Ze =o = A 1 (8.5) 


where F,’ = F, when xo’ = 0 and where AF,’ is the contribution of x’ to F. 
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(6) M/(M — 1) 2 — th’ = M/(M +1). From Egs. (7.4), (7.5), and (8.1) and Table 2, the lift coefficient is 
Cre = 2(t' — bh’) [So I/D (dy — Go) fi(vo)dvo + al rt x’ (de — Go) fo(v0)dvo] 
= 2(t' — t0)?(ér — do) [J~9/ vodvofilvo) + So” v0 fo(vo)dvo] 
= (d — do) Fo(t’ — te’) (8.6) 


where F, can be evaluated by integrating Eqs. (7.4) and (7.5). 


1—(1/M) 1/(t’ —to’) 4 t’ sald te’ 
F,(t’ — to’) = 2(t’ — to’)? Fi vodvofi(vo) + rofdan | = r 7 ([ + 
0 1-(1/M) (M? — 1)” 2 
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; Tv 




















4(M? — 1) 4M?n M? - 1 
(M2 — mF eo f td ( 1 )) 
yi SING a BSS. ke? HD conan ff — fe’ -- — cones —1 ae nce Pe 
V(t! — ’)? — Mt’ — hh’ — 1)? + roa i+ 5+ peo Mm - se 
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(uw? — 1)'*\2 x lm ‘ Mr peer) Be. 
(M? — 1)” ‘ , 
a V(t! — ty’)? — M(t’ — th’ — i) = F,’ — AF,’ (8.7) 


where fF,’ = F, when x’ = 0; AF)’ is the part of F, contributed by x0’ ¥ 0. 


(c) OZ t' — th’ Z M/(M +1), From Eggs. (7.4), (7.5), (7.6), and (8.1) and Table 2, the lift coefficient is 
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1+(1/M) 
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0 1—(1/M) 1+(1/M) 
1 4 2 Axy’ 
ee ee ye Re ey an ee Se YX 8.9 
mm oe + a ee we . om 


where F;’ = F; when xo’ = 0 and AF;’ is the part of F; contributed by xo’ ¥ 0. If we introduce j = 1, 2, 3, F;’ 
and AF,’/xo’ are given in the above equations. F,’ and AF;’/x9’ vs. t’ — #,’ are plotted in Fig. 10 for M = 1.5, 2, 
3, and 4. 


(d) Similarly, we can have the moment coefficient for the three cases. 
Gany = (61 — do) [G,(t’ — to’) — x0'Fj(t’ — to’)] (j = 1, 2, 3) (8.10) 
When t’ — th’ > M/(M — 1),j = 1, 
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(t! — &’ +1) V(r — &)? — Mt’ — &’ — i) = G,’ — AG,’ (8.12) 
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In the above, it should be noted that G,’ = 
(j = 1,2, 3). 


It is interesting and useful to note that AF;/xo’ and 
AG,/x0" can be used to calculate Cz and Cy of an airfoil 
or a wing that starts suddenly to descend at constant 
vertical velocity h. The procedure is: Let xo’ = 1 in 
the above expressions of AF; and AG;. Then AF; = 
C,;/(h/U) and AG, = Cy,/(h/U) of the wing because of 
the flapping degree of freedom. Of course, some integ- 
ral relations may be easily formulated for the case of h 
being an arbitrary function of time. 
discussed in a later paper. 

(e) Decomposition of Cz: into three parts.—Physically 
Cy(t’) can be decomposed into three parts—namely, 
ACzs, ACzrs, and ACze. ACz,4 is the contribution of 
instant angle of attack Aa = (& — do)(t’ — to’) and 
can be calculated from Ackeret’s formula 


Cra = 4Aa/(M? — 1)” (8.14) 


The second part ACzz is the contribution due to rotating 
velocity @, — da and is proportional to & — do but 
independent of t’. Its value can be shown as 
4(a, — do) (t’ — to’) 
(M? — 1)” 
4(t’ — to’) 
(Mm? — 1)'7 


= a! [a — do) 


(8.15) 
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The third part ACz,¢ is the contribution due to transient 
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G; when xo’ = 
G,’ and AG;,’/xo’ vs. (t’ — t,’) are plotted in Fig. 11 for M = 1.5, 2, 3, and 4. 


The details will be ° 


yA 
= to’)? oe =| 
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0 and AG,’ is the part of G; contributed by xo’ ~ 0 


(8.13) 





effect that occurs when 0 Z ?’ 
If0 Zt’ — th)’ Z M/(M + 1), 


— th’ Z M/(M — 1) 


ACzre = Cris me Cru = (ay ai do) (F — F,) (8.16a) 
and if M/(M +1) Zt’ —’ Z M/(M — 1), 
ACic = a Cru = (a _ do) (Fe — F,) (8.16b) 


where Fi, F:, and F3 are functions of t’ — fo’ and are 
given in Eqs. (8.5), (8.7), and (8.9). 
Similarly, Cy4 can also be decomposed into three 


parts, 


ACwwa = [4Aa/(M? i = xo'] (8.17) 


1)'/*] (1/2) 


which can be evaluated from Ackeret’s formula if the 
instant increment of a is used. 


4(é — do) | 
(M? — 1)” 


M?—2 _ x9’ 
J es =| —_ xo’ ACrp 
3(M? — 1) 2 





ACus = 
(8.18) 


which is proportional to (a — do) but is independent of 
time, and the transient part is 


ACuc = Cyr — ACys — ACup (8.19) 
Mis © Vy ~~ GQ ~ 2 ~ FN 

if M/(M+ 1) Zt’ — tb’ Z M/(M — 1) or 
ACuc = (G1 — &)[G3 — Gi — x0’(Fs — Fi)] 


if0 Zt’ — ZM/(M +1). The discussion of the 


results will 2 treated later. 
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TABLE 3 
Cr/a;, and (Cu + xo’Cr)/a as functions of t’ — ft)’ and t’ — t,’ 
M | MM M a M 
Range of t/ — f’ i =k > 2% , ean <i Oo S St -4' ss = 
vane unge oO ly’ > Wni | wet <t bh’ < Wo | O<t bh’ < Worl 
’>o Zones crossed I-3 I-1 and [-2 I-1, I-2, and I-3 
Cr/a = F,(t’ — to’) F(t? — to’) F(t? — t’) 
(Cu + x0’Cr)/a = Gilt’ — by’) Go(t’ — to’) G3(t! — to’) : 
: ae + cua Pier 
‘= hh! ‘-' > ZF fhe oe ik 2 wo 
oo otf | M ra 1 " VW a M 3 
v— t’ | t—t'> : ——— ¢ t’ — th’ F< ———_ CSF —4' $5 — 
- 7s i M+1 M—1 U+i 
Sones erensed I-1 and II-1 I-1, II-1, and II-2 I-1, II-1, II-2, and Ig 
Cy fan = | At = af F(t’ — hh’) Flt! — &') — Flt’ — 4’) F(t’ — &’) — F(t’ — ty’) 
| ee 
| ane 
Regions [ (Cy + 'Ci)/m = | Gilt’ — ’) — Gilt’ — 4") | Gilt’ = 6) — Gilt’ — t') | Gilt’ = be") — Gat’ = 4/) 
I+ II 2t, 
> 0,’ = Of — ih 
th a: Ly SOT ae a ey Load eek ae 
fe —& | wes ~ 8 Sa | Weis SRT | OS - HST 
Range of M WU | Ww va 
ited tae 20 +4 S— <7 — 4 << <0 —4' <— 
M+1 =i | M+1 M+1 
Zones crossed I-1, 1-2, II-1, and I1-2 | I-1,1-2,1I-1,11-2,andII-3 | I-1, 1-2, 1-3, 11-1, II-2, and 
Cja= | Ai’—-#)- Ft’ —4') | A’ — 6’) — A’ —4’) | Al’ — 6) — Alt’ —4) 
(Cu + x0’Cr)/a = G.(t’ — te’) — Gilt’ — b’) G2(t’ — bh’) — Gilt’ — ty’) Gilt’ — bh’) — Gilt’ — &’) 
. 











(f) The second region of the given downwash.—In the 
case of another change in a slope at ty’ (i.e., d& — &), 
another new region, II, is formulated. If t’ > 4’, both 
regions will contribute to lift and moment. As the con- 
tribution due to & — da on lift and moment has been 
given above, the additional contribution due to @& — a 
in region IT is 
Craj(t’) = (@2 — &) F(t’ — hh’) (fj = 1,2,3) 


é,)(G,(t! — th’) — 
xo’ F(t’ — t’)] (8.21) 


Cyyaj(t") = (a, — 


where F;(t’ — t;’) and G,(t’ — t,’) are the same as given 
in Eqs. (8.5), (8.7), (8.9), (8.11), (8.12), and (8.13) ex- 
cept to’ is replaced by ty’. 

For the case d& = a = 0, the result is given in Table 
3. There are a number of interesting points. 

(1) Inregion I (0 Z ¢ Z t,’) where disturbance due to 
the second change in & at t = ¢, has not started, C, and 
Cy are due only to the first change in & at t’ = 0, and it 
is obvious that the second change can contribute noth- 
ing. 

(2) In regions (I + II), (t’ S 4’), changes of a& at 
both ¢’ = 0 and at ¢’ = 4,’ contribute to lift and mo- 
ment. The contribution. of the second change in 4, 
however, may be considered as independent of that due 
to the first change. In regions (I + II), one contribu- 
tion superimposes on the other. This shows that the 
principle of superposition can apply not only to the 
velocity potential in a linear differential equation but 


(8.20) - 


also to its derivatives or to the finite integrals of its 
derivatives such as C; and Cy. 

(3) If the chord crosses the zones I-1 and II-1, the 
transient effect dies out and the angle of attack is a. 
Cy = 4a/(M? — 1) and Cy = 2a:/(M? — 1)" 
check with Ackeret’s formulas for the steady case. 


(B) Extension to an Arbitrary Change in a at t,' 


Suppose that at ¢,’ there is a sudden change in a 
slope—i.e., &;+1 — &, as shown in Fig. 7— and then a 
slope remains constant at a+, after ¢,’. From the 
above analysis, it is not difficult to see that the contribu- 
tion to C, and Cy due to &;+; — &, at t,’ must be of the 
same kind as that due to & — dp at fo’ in region I. 
Similarly, we can write 


[Crit y(t’) ], = (Qi+1 — a) F(t’ _ 2) (8.22) 
and 
[Carig+ y(t’) ] 5 = (Qi+1 — a;)(G,(t’ —t,/)- 
xo’ F,(t’ > t,’)] (8.22a) 


where j = 1, when ¢’ — t,’ > M/(M — 1); 7 = 2, when 
M/(M+1) Zt’ —t! ZM/(M — 1); andj = 3 when 
OL —t' ZM/(M+1). F(t’ — ty’) andG,(t’ - 
t,’) are the same as F(t’ — fo’) given in Eqs. (8.5), (8.7), 
and (8.9) and G,(t’ — ¢’) given in Eqs. (8.11), (8.12), 
and (8.13), respectively, except ¢,’ replaces fp’. Obvi- 
ously, the contribution of &;4; — &, to Cz and Cy is zero 
et <8’: 
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(€) Total Lift and Moment at t’ for a Given a(t’) 


Before extending known results to this general case, 

some behavior of a(t’) and &(t’) curves must be clearly 
understood. 
» (1) a(t’) vs. t’ is decomposible into elements linear 
with time. For a given a(t’) where a is a continuous 
function of ¢’, the curve can be approximated with a 
series of infinitely small tangent line elements. If.we 
extend those line elements in the downstream direction, 
a series of wedge areas is formed. As an example, 
Aa(t’ — t,’) starts at the conjunction ?,’ as shown by the 
dotted lines. The contribution of this Aa(t’ — ¢,’) to 
the force or moment must be proportional to that a(t’) 
in region I of Fig. 2. In that case, its contribution is 
independent of later changes in &. Similarly, its con- 
tribution is also independent of early changes in &. 
a(t’) at any time can be obtained by summing up all of 
the wedge elements. 

(2) &(t’) vs. t’ is decomposible into step functions of 
time. The corresponding a(t’) vs. ¢’ and its approx- 
imation step curves are plotted in Fig. (Sb) It is also 
decomposible into a series of step curves which extend 
to infinity at constant value as shown in Ad; = a +1 — 
ay at i 

(3) There are three types of regions ahead of ?’. 
Suppose we are interested in C, and Cy at?’. Two lines 
can be drawn from the trailing edge of the chord at ¢’, 
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one with slope 1//(M — 1) and the other M/(M + 1) 
as shown in Fig. (8c). It has been learned from earlier 
results that any region ahead of time t’ — M/(M — 1) 
is covered by the chord in only one zone of type I-1. 
Similarly, any region between t’ — M/(M — 1) and 
t’ — M/(M + 1) is covered by the chord in two zones 
of types I-1 and I-2, and any region between /’ — 
M/(M — 1) and ?’ is covered in three zones of types I-1, 
I-2, and I-3. 

From above, we know that the contribution of each 
change in &- (say Ad; = a;+1 — a) can be treated inde- 
pendently of changes in & behind or ahead of ¢;. The 
total contribution must be the sum of the individual 
contributions. If t’ > M/(M — 1), the regions ahead 
of t’ are divided into three groups as follows: 

First group: k& + 1 regions with only one zone cov- 
ered by chord at ?’. 


ty’ < t’ — [M/(M — 1)] < tety’ 


Second group: m — & regions with two zones 


covered. 
tm’ < t’ — [M/(M + 1)] < tnt’ 


Third group: » — m regions with three zones 


covered. 


b,’ pA t’ 4 bata’ 


Based upon the principle of superposition, we have the following three categories: 


(a) t’ > M/(M — 1). 


k+l1 m+1 n+1 


Subscript III refers to existence of all three groups. 


Cr(t") on = > Cra Ma 2», Cree 5 p Cras 


i=1 ti=m-+2 


—- > (ai+1 — G4) Fy a 4  aets — a) PF, + > (ait+i — a4) Fo (8.23a) 
i=0 i=m+l1 
where all F’s are functions of ¢’ — ¢,’ as given in Eqs. (8.5), (8.7), and (8.9). Similarly, 
k+1 m+1 n+1 
Cult) mn = ) Cun + 2 c m2 >» Cua 
i=] i=k+2 i=m+2 
k m 
ag p (ai+1 — a)Gi + > (aj+1 — &)G2 + es Gi41 — &)Gs3 — X0'Cr(t) ny (8.24a 
i=0 i=k+1 pts 
(6) M/(M — 1) 2 t’ > M/(M + 1). In this case, the regions in the first group do not exist, and it can be 
shown that 
CLht)n = 7 (ai+1 — a) Fe + DO (aeeti — Gy) Fs (8.23b 
i= i=m+1 


Cul(t)n = ba (Giga — G)G2 + > (Geta — e1)G3 — Xo’ Cr(t’) 11 


#=m+l1 


(8.24b) 


where the subscript II refers to the existence of only two groups. 
(c) M/(M+1) 2 >0. 


Cit’); = > (as+1 — Gy) Fs 


Cult’): = du (ait1 — &1)Gs 


In this case, the regions in both the first and the second groups do not exist, and 


(8.23c) 


(8.24c) 
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Now for the case where a(t’) is continuous to the second derivative, the above can be expressed in terms of in- 


tegrals. 


If the number of the tangent line elements increases indefinitely and each decreases to infinitely small 


length, we find in the limits, t,’ — t’ — [M/(M — 1)], tm’ >t’ — [M/(M + 1)], and except a: — a at t = 0, 


which is discontinuous, all a;+1 — &can be approximated by &, At,’ (where At,’ = t,+1’ — t,’). 


ing to the above three cases, we have: 


(d) t’ > M/(M — 1). 


¢’— (M/(M —1)] 
Cr(t’)ma = (ds — do) Fi(t’) + 4 &(r’)Fi(t’ — r’)dr’ + 


where 7’ is introduced to replace ¢;’. 


t’—[M/(M—1)] 
Cult) = (da — do)Gi(t’) + 4 &(r’)Gy(t’ — r’)dr’ + 
0 


(ec) M/(M +1) <t’ S$ M/(M - 1). 


1’ —[M/(M +1)] 1 
Cit’) 1 = (ay _ ao) F(t’) + f a(r’) Fe(t’ = t’)dr’ +f 
0 t’—(M/(M +1)] 


1’ — [M/(M +1)] 1’ 
Cult’)1 = (da — do)Ge(t’) + i &(7’)Go(t! - r’)dr’ +f 
0 t’—[M/(M +1)] 


Thus, correspond- 


t’—[M/(M +1)] ; 
&(r’) Fo(t’ — r’)dr’ + 


t’—[M/(M—1)] 


t’ 
f a@(r’)F3(t’ —-r’)dr’ (8.25a) 
t’—[M/(M+1)] . 
t’—[M/(M+1)] 
&(r")Ga(t!-— r’)dr’ + 
t’—[M/(M —1)] 
1’ 
f &(r’)G3(t’ = t’)dr’ — xo’ Cr(t’) 1 (8.26a) 
t’—[(M/(M+1)] 
a@(r’) F(t’ — 1’)dr’ (8.25b) 


&e(r’)G3(t’ — t’)dr’ - xo’ Cr(t") 1 


(8.26b) 


where in the first term F;(¢’) and G2(t’) replace F,(t’) and G,(t’), respectively. 


(f) MM +1) >t >0. 
Cr(t’); = (a1 — do) F3(t’) = Sar’) Fs(t’ — 7’)dr’ 


(8.25c) 


Cult’); = (é1 — do)G;(t’) = 


+ fol a(r\Gt’ — 1’)dr’ — x0'Cz(t’), (8.26c) 


It is easy to see that C,(t’ = 0) and Cy,(t’ = 0) are not 
zero in general but 


C,(0) = (a: — do) F3(0) (8.25d) 


Cuy(0) = (a1 — d&o)G3(0) — x0’Cz(0)  (8.26d) 
The above integrals are similar to Duhamel’s integral'® 
or to Carson’s integral. 

’ All the above formulas are true if a = 0 for a length 
of time > C/(U — c) ahead of t’ = 0, because the past 
transient effect will die out before ¢ = 0. ap is not 
necessarily zero; its contribution to Cy is 4ao/(M? — 
1)’ and to Cy is 2a/(M? — 1)’*. Both should be 
added to the above if ao ¥ 0. 

In case of & being continuous except at a finite num- 
ber of discontinuities in a@, the above can also be 
applied. 

The above analysis shows one important physical 
principle. If the solution of an elementary boundary 
condition is known and if the given boundary condition 


can be obtained by superimposing the elementary bound- 
ary condition, the general solution can then be ob- 
tained by superimposing the elementary solution. 


(D) Some Applications 


As shown before, Eqs. (8.5), (8.7), and (8.9) can be 
applied to determine the contribution to C;, due to 
ai+1 — a; at t,’. For this purpose, Fig. (10a) shows 
F,’ vs. t’ — t,’ for M = 1.5, 2, 3, and 4. Herej = 1, 
when ¢’ — t,;’ > M/(M — 1); j = 2, when M/(M + 1) 
Zt’ —t,) Z M/(M — 1); andj = 3, when 0 Z t’ - 
t’ Z M/(M + 1). From the definition Fj’ = C,;;/ 
(ait+1 — &;) if xo’ = 0. If xo’ is not zero, AF;’/xo’ vs. 
t’ — t,’is givenin Fig. (10b). Similarly, G; and AG;,/x9’ 
in Eqs. (8.11), (8.12), and (8.13) are plotted in Fig. 11. 
Both figures are useful for calculating C; and Cy for an 
arbitrary a(t’) which is continuous and can be approx- 
imated by straight line elements. 

As an example, the following case of the downwash 
distribution in Fig: 2 is illustrated: 


C = 6 ft. c = 1000 ft. per sec. 
t, = 0.02 sec. (t;’=5) xo’ = Oandx’ = 0.5 


Figs. 12, 13, and 14 show C;/Czyo, Cy/Cr, and center 
of pressure from the leading edge vs. (t’ — fo’). For 
illustration, the steady lift coefficient Cy = 4a;/(M? — 
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Fic. lla. G;’ vs. t’ — t;’ (xe’ = 0). 24 
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Fic. 1lb. AG’/x9’ vs. t’ — ti’. 




















Fic. 14. Center of pressure from leading edge vs. t’ — ¢;’._ In 





itial angle of attack is assumed to be zero. 





oe 






Fic. 12. (a) Cr/Crevs. t’ — t&’ for the case a(t’) as shown in Fig 







2. M = 1.5, 4 = 0.02 sec., C = 6 ft., c = 1,000 ft. per sec. 
Solid line—xs’ = 0. Dotted line—x’ = 0.5. Note: ACza/- 





Cre, ACitB/Czre, and ACric/Crw in both xe’ = Oand xo’ = 0.5, 
(b) ACic/Cre vs. t’ — to’. Note: ACrc/Cr at t’ — ty’ is 
equal and opposite to ACrc/Cz at t’ — t)', if t' — bh’ = 0’ — be’. 
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Fic. 15a. 
2, 3, 4. . Solid curves—O st’ — 2;' S M/(M + 1). 
curves—M/(M+1) st’ —#;' S M/(M — 1). 


ACic/(a@i+; — a) vs. t’ — tj’ at xo’ = Ofor M = 1.5 
Dotted FIG. 


5b. ACre'’/x9’(ai+1 — a) vs. t’ — ti’ (ACzc’ is additional 


contribution of ACzc due to xo’). 












































[~ —a 
| TABLE 4 
Some geueral features in Figs. 12, 13, and 14 
Region 1, 0 Zt’ < ty’ | Region I + II, t’ = &’ 

Xo! | a: | 0.5 0 0.5 
Cia/ Cro | ze | fe } 1 1 
h’ 7 
Cra/Cio | 0.02 —().08 0 0 
Crc/Cuo Fig. (12b) | Fig. (12b) Cricl(t’ — th’) = —Cre(t’ — to’) | Cre(t’ — th’) = —Cic(t’ — b’) 
Cua/Cro | Ya | 0 1 0 
Bite eo Es. See Sf id Nelo OE, A SE Re ey Ie 
Cus/Ciuo | 0.0133 0.003333 0 0 
Cuc/Cuo | Fig. (13c) | Fig. (13c) Cuc(t’ — th’) = —Cuclt’ — to’) Cuc(t’ — th’) = —Cuclt’ — to’) 
Center of | «—— ———_____-___-_—— -~-—-Fig. 14—- ——-—__-______- —____--____ ________-_- > 
pressure 





1) is substituted for a so that this result can be used for 
any a. Of course, this result is equivalent to C,/d 
and Cy/é with &; = 15.625 rad. per sec. A more com- 
plete explanation will be found in Table 4, and some 
discussion will be given in the conclusion. © Fig. (15a) 
shows ACzc/(ai+1 — &;) as a function of t’ — #t,’ for 
M = 1.5, 2, 3, and 4 for x)’ = 0. Fig. (15b) gives the 
* additional contribution to ACze due to xo’ ¥ 0, and this 
additional contribution may be called ACz¢’. 


(9) CoNcLUSIONS 


(1) For an airfoil in a flow of constant supersonic 
speed, the transient reaction will damp out in a time 
period C/(U — c) immediately after the change in angle 


of attack ceases. After that period, the lift, wave drag, 
and moment become the same as given by Ackeret in 
the two-dimensional case. The transient reaction be- 
comes more pronounced and lasts longer if the free- 
stream velocity approaches the sound velocity or if the 
Mach Number approaches 1. 


(2) The effect on C, (or Cy) due to an increment of & 
at a certain time can be decomposed into three parts: 


ACra> ACzz, and ACie (or ACwa; ACys, and ACuc). 
(a) ACzs, or ACya is directly proportional to 
the instant increment in a and can be evaluated 
from Ackeret’s formula. The value of ACz, is 
independent of the hinge axis, but the value of 
ACya is not. 





1 

(5) 
great 
imme 
(Refe 
sure 1 
forwa 
more 
analy 
at the 
of pre 
(4) 
tionst 
ent at 
plicat 
throu, 
condit 


Wit 


and 


With 1 


Uxt, 
3 4 





ional 





TRANSIENT REACTION OF AN AIRFOIL’ 653 


(b) ACzs or ACysz is directly proportional to’ Aé 
and independent of time. The value of ACz, de- 
creases linearly with increase of xo’, the hinge axis. 
(c) ACze or ACye, the transient part, is propor- 
tional to A& but is a complicated function of time 
and Mach Number. ACy¢ decreases linearly with 
increase of xo’. ACyc decreases with increase of 
xo’ but not linearly. The effect of Mach Number 
on ACzye or ACye is more pronounced as M ap- 
proaches 1. (See Fig. 15.) For extremely large 
Mach Numbers, the contribution of ACze or ACye 
is small and may be neglected in certain cases. 

(3) Immediately after increasing a(t), C, and Cy are 
greater than the values that are linear with time, and 
immediately after decreasing a(t), C, and Cy are less. 
(Refer to Figs. 12 and 13.) Besides, the center of pres- 
sure moves backward after increasing a(t) and moves 
forward after decreasing a(t). The smaller the lift, the 
more pronounced the shift. (See Fig. 14.) From this 
analysis, &(t) should be kept small at the beginning and 
at the end of changing a, so that the travel in the center 
of pressure may not impair the stability of the missile. 

(4) In the linearized theory, there is a unique rela- 
tionship between C,(f) [or Cy,(¢#)] and a(t). The pres- 
ent analysis shows that in the case of a(t) being a com- 
plicated function of ¢t, this relation can be established 
through.superposition of solutions of a simple boundary 
condition—1.e., Aa(¢) being constant. 


(5) It is not difficult to calculate the lift and moment 
for the case of an airfoil in harmonic oscillation by 
means of the generalized integrals of Cz, and Cy. It 
will be treated later. 


(6) This result can be applied to study of the transi- 
ent reaction of the deflecting control surface if the main 
wing remains stationary. With the principle of super- 
position, the present analysis can be readily applied to 
more degtees of freedom. The details will be given in a 
later report. 


(7) With this analysis, the lift and moment of an air- 
foil can be calculated for any arbitrary a(t) without 
difficulty. More treatment seems necessary for the 
case with more degrees of freedom, the case of a wing 
with finite span, and the case with accelerated flow. 


(8) This analysis is based upon the linearized theory. 
It cannot be applied when M is near 1, and when is 
larger than |U/o). 

(9) It is obvious that the result can be applied to an 
infinite swept wing if the leading edge is ahead of the 
Mach line. In that case M should be replaced by M 
cos y and a by a/cos y in all of the above equations. 
(y = swept angle.) 

(10) The expressions of AF, and AG, can be easily 
used to calculate C, and Cy of a wing due to flapping 
degree of freedom. 





Appendix A 


DERIVATION OF THE VELOCITY POTENTIAL 


With the aid of reference 6, it is not difficult to show that the following integrals are true: 


( sin-' (b/lal) + Va? — 6? — |a| _—for a} = Id 


> dw . 
ri — (1 — cos bw)Jo(aw) = (Al) 
aes ( (/2)|0| — la for |a| Z |b 
and 
eile sin—! b/|a| for |a| = |d! 
f — sin bwJ (aw) = (A2) 
ite (x/2)(b/|d!) for |a| Z |d| 


With Eqs. (Al) and( A2), we can integrate Eq. (6.8) to the form 


— ay =-3 
g(x, —0, t) | | = 
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(: ie ire cng, Py pt, ir) 
—i+tt - M dX) sin " ae” + 


ci x/c(M?—1) SEs | 
ae aad anx< Vx? — (@—h — Mr? when! — MZ 1 (A3) 








lt -—t4 — MX 142") So. ET, I 
Mi horacut j—~% — Wr a oe 
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where two substitutions are used—namely, \ = (x — £)/c(M* — 1) andt, — mX = t;. If tis replaced byt — 4, 
in the first integral, the first integral becomes the second with opposite sign. Eq. (A3) is a discontinuous integral, 
The integration can be carried out only between points of discontinuity. 

For convenience of integration, introduce Iy for the upper line of the first integral with indefinite limit, I, for the 
lower; II yfor the upper line of the second integral with indefinite limit, and II, for the lower in the above equation. 
It is easy to show that the values of the indefinite integrals are: 


oe ee 
Ivy = (: +t, — =a ta . r) \ sin= (; - at) nal mae (ts) 1/ sin=! [a = 2 Ds) + 
M N (Mw? — 1)'7 t 


é? 1 r (M2 — 1) }? t t 
i La Slay ig *,/ — [yO for 4 tq, 
2(M? — 1) “A ge M . [ar i | M+1 u-i 


[> ZrZz t—kt 























IIly = —livif tis replaced byt — ht; wai 71% 4 ipo (A5) 

I, = +2(t+ te - ya for d <— ;and t > @ (A6) 

II, = wae if t is replaced by ¢ — forrA < oo andi—t >0 (A7) 
Introducing 

I,’ = -I, for \ > : “a <0 (A8) 

iif oie wiry 4<8 (A9) 





In order to express ¢ more explicitly, we introduce ®,; (¢ = 0, 1, and 7 = 1, 2, 3) to denote the definite integrals 
in different cases of Eq. (A3). Thus, 


&o, = cU(M? — 1) fy 1,}™" x/e(M3—=1) _ =|? tt. |) 1 = | weer (A10) 


T Vo 2(M? — 1) 





where v9 = x/U(t — to) and voo = x0/U(t — to) (to being zero). Note that ast 0,» — ~ and ast— ©,» 0. 


cU(M? — 1) 





—Gy = ———_— il, > tactile for vo < 0 (All) 
T 
_ cU(M? — 1) ys d= t/(M +1) de=2/c(M2—1) | = (7 — v0 1 e m 
® —_——_—_—— ie A =— —- — —(-—M)— 
02 = - [x] + [olijc'¢1) f pe se ar — D sin : 





1 tf y packet el m - . 
ap yen (1 ~ Z)[ sna — om) — 5] 301 + ge) vim ap) cam 


where \ = t/(M + 1) isa point of discontinuity. 


U(M? — 1 i are sais 
ag ot ( = t1} MHD {1}? M1) (1, > (eI) 





t/(M +1) t/(M — 1) 
ti a) -2 =s 
eo —)" “(1 —_ “) - = __ : \ for 1 + 2 ZuyZ®” 
Mv me a 0 M 


where \ = ¢/(M + 1) andA = t/(M — 1) are points of discontinuity and 
@o3 = [(M? — 1)/2/M](x?/v0) [1 — (v00/v0)] (A13) 
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In those cases vp — 14, and voo — v0. 


U(t — 4). 


1 Possio, C., L’azione aerodinamica sul profilo oscillante alle 
syelocita ultrasonore, Acta, Pont. Acad. Sci., Vol. 1, No. 11, pp. 93- 
105, 1937. i 

* Garrick, I. E., and Rubinow, S. I., Flutter and Oscillating Air- 
Force Calculations for an Airfoil in a Two-dimensional Supersonic 
Flow, N.A.C.A. T.N. No. 1158, October, 1946. 

yon Borbely, S., Aerodynamic Forces on A Harmonically 
Oscillating Wing at Supersonic Velocity (2-Dimensional Case), 
R.T.P. Translation No. 2019, British Ministry of Aircraft Pro- 
duction (from Z.A.M.M., Bd. 22, Heft 4, pp. 190-205, August, 
1942). 

‘Schwarz, L., Ebene instationare Theorie der Tragflaéche bei 
Uberschallgeschwindigkeit, Jahrb. 1943 der deutschen Luftfahrt- 
forschung, IAOIO, Stephan Geibel & Co. (Altenburg/Thur), pp. 
1-8 (available as Translation No. F-TS-934-RE, Air Materiel 
Command, Wright Field, Dayton, Ohio, March, 1947). 

5 Kussner, H. G., General Airfoil Theory, N.A.C.A. T.M. No. 
979, 1941. 

6 Watson, G. N., A Treatise on the Theory of Bessel Functions, 
2nd Ed.; The MacMillan Company, New York, 1944. 

7 Garrick, I. E., and Rubinow, S. I., Theoretical Study of Air 
Forces on an Oscillating or Steady Thin Wing in a Supersonic 
Main Stream, N.A.C.A. T.N. No. 1383, July, 1947. 


TRANSIENT REACTION OF AN 


Similarly, we can evaluate ®y, Py, and &,; by replacing t — f& by t — h. Let » = x/U(t — th) and m0 = x0/- 
A useful relation may be noted: 


Oo, = Oy + (Uxt,/2) 








AIRFOIL 





655 


(Al4) 


REFERENCES 


8 Miles, J. W., The Aerodynamic Forces on an Oscillating Airfoil 
at Supersonic Speeds, Journal of the Aeronautical Sciences, Vol. 
14, No. 6, pp. 351-358, June, 1947. 


9 yon Karman, Th., and Sears, W. R., Airfoil Theory for Non- 
Uniform Motion, Journal of the Aeronautical Sciences, Vol. 5, 
No. 10, pp. 379-389, August, 1938. 


10 yon Karman, Th., and Biot, M. A., Mathematical Methods in 
Engineering; McGraw-Hill Book Company, Inc., New York, 
1940. 

‘yon Karman, Th., Supersonic Aerodynamics—Theory and 
A pplication, Journal of Aeronautical Sciences, Vol. 14, No. 7, pp. 
373-409, July, 1947. 

12 Ackeret, J., Luftkrafte auf Fliigel die mit grésserer als Schall- 
geschwindigkeit bewegt werden, Z.F.M., Vol. 16, pp. 72-74, 1925. 

18 H6nl, H., Zweidimensionale Tragfléchentheorie im Uber- 
schallgebiet, Forschungsbericht Nr. 1903, Deutsche Luftfahrt- 
forschung (Géttingen), 1944. 


14 Temple, G., and Jahn, H. A., Flutter at Supersonic Speeds, 
Part 1. Mid-Chord Derivative Coefficients for a Thin Airfoil at 
Zero Incidence, Rept. No. S.M.E. 3314, British R.A.E., April, 
1945. 





Dear Sir: 

Professor Alan Pope (‘‘Letter to the Editor,’’ JoURNAL OF THE 
AERONAUTICAL SCIENCES, Vol. 15, No. 8, p. 496, August, 1948) 
recently gave a simplified calculation for the power-saving obtain- 
able by sucking off the boundary layer and accelerating it up to 
free-stream conditions. The result of his simplified calculation 
was a net saving of 50 per cent of the power required. 

A more realistic calculation was indicated in the article by 
Professor Goldstein (‘‘Low Drag and Suction Airfoils,’? JouRNAL 
OF THE AFRONAUTICAL SCIENCES, Vol. 15, No. 4, p. 209, April, 
1948). Carrying this through for a flat plate is easily done and 
yields much less favorable results. 

Thus, if 5 is the height of the boundary layer; wm, the free- 
stream velocity; u, the velocity in the boundary layer; and p, 
the constant density, then the drag of the plate (two sides) is 


D = 2%» f,>u(ue — u)dy 
and the power required to overcome this drag is 


P= Qpue f5°u (10 — u)dy 

If the boundary-layer air is accelerated up to free-stream condi- 
tions, the rate of kinetic energy increase will equal the power re- 
quired to accelerate the air, and the rate of momentum increase 
will yield a thrust equal to the drag. The power thus required 
assuming a lossless system) is 


P’ = 2(p/2) fo >u(ur — u*)dy 


Letter to the Editor 


It we now assume that the boundary layer profile is given by 





u/ue = (y/d)” 
Then 
P = 2puo*5(n/(n + 1)(2n + 1)] 
P’ = puo*5[2n/(n + 1)(3n + 1)] 
so that ; 
P’ Qn + 1 B 2+ (1/n) 
3n +1 3 + (1/n) 


For a crude approximation to the laminar boundary layer take 
n=1. Then P’/P = 0.75 ora power saving of 25 percent. For 
the turbulent boundary layer take n = 1/7, whence P’/P = 0.9 
or a power saving of only 10 per cent. Professor Goldstein men- 
tions comparable savings of 21 and 8 per cent, respectively. 

It is gratifying that at least at one point in the boundary-layer 
control system a net gain is possible. Probably this gain will be 
more than cancelled by duct losses, etc. 

It can further be noted that the relatively small gain obtained 
by pumping off the turbulent boundary layer is typical. Suction 
of airfoils with predominantly turbulent (but unseparated) 
boundary layers will yield relatively small decreases in the profile 
drag. The large gains generally pictured as obtainable by 
boundary-layer control are possible only by delaying transition to 
far aft chord positions. 


ABRAHAM FIuUL 
Aerodynamics Engineer 
Northrop Aircraft, Inc 























Vorticity in the Supersonic Flow About 
Yawing Cones | 


MARK M. LOTKIN* 


Ballistic Research Laboratories 


ABSTRACT 


A comparison is made between Stone’s and Sauer’s methods of 
calculation of the supersonic flow about a slightly yawing cone. 
It is shown that the neglect of vorticity, while having only a small 
effect on the ratio of shock yaw to cone yaw, may lead to erroneous 
values of the normal force coefficient. 


INTRODUCTION 


Spires METHODS dealing with the flow of air about 
a slightly yawing cone have been developed. ' We 
shall concern ourselves here with the treatment of this 
problem as given by Stone! and compare it with that 
advanced by Sauer.? Of these two methods, Stone’s is 
the more exact one, since it assumes only—among other 
things—isentropy along streamlines; whereas Sauer 
bases his deductions on the premise of irrotationality, 
and, hence, isentropy, in the whole region of flow. A 
comparison, then, is of interest mainly because it per- 
mits the estimation of the error introduced in the 
assumption of potential flow. 

Following brief outlines of both theories, the com- 
parison mentioned above is made by an evaluation of 
two quantities obtainable by both methods—namely, 
the ratio of shock yaw to cone yaw, and the coefficient 
of the normal force acting on the head of the cone. 





NOTATION 


D = diameter of the base of the circular cone 

Ky = coefficient of the normal force on the cone 

M, = Mach Number of the free stream 

N = normal force of the cone 

S = entropy of the flow 

Se = entropy of the flow in the nonyaw case 

T = absolute temperature 

X, Y,Z = components of the velocity vector g* in the 


direction of 7, 0, w, respectively, measured in 
units of ¢ 
U(0), V(@) = components of g* in the nonyaw case 


a = local speed of sound in the nonyaw case, meas- 
ured in units of ¢ 

c = speed of efflux into vacuum 

d, = variation of entropy in the plane of yaw 

p = local pressure of flow 

q*,q = local velocity, speed of flow, respectively 

r = radial distance from the vertex of the cone 

wy N = perturbations of X, Y in the plane of yaw 

a = perturbation of Z normal to the plane of yaw 

da, SK = error in a, Ky, respectively 

0 = polar angle in the system of spherical coordinates 

Os = semiapex angle of the cone 


Received May 10, 1948. 
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Ox = shock wave angle in then> yawcase . 

& = potential function of the oblique flow 

on = yaw of the shock wave 

8 = angle of attack 

y = ratio of the specific heats of air 

6 = shock wave angle 

« = angle of yaw of the cone 

m1 = perturbation of p in the plane of yaw 

v = angle of rotation of the cone in Sauer’s theory 
& = perturbation of p in the plane of yaw 

p = density of the flow 

¢ = potential function of the flow in the nonyaw case 


= azimuthal angle in the system of spherical co- 
ordinates 


STONE’S THEORY 


In his theory of the yawing cone Stone makes the 
following basic assumptions: 

(a) The Mach Number 
sufficiently large, and the semiapex angle @, of the cone 
is sufficiently small for the (conical) shock wave in the 
nonyaw case to be attached to the vertex of the cone. 

(b) The flow is steady and independent of the dis- 
tance 7 from the vertex of the cone. Then, velocity q* 
of the flow, pressure p, density p, etc., are constant 
along radial lines through the vertex. 

(c) The flow is adiabatic along each streamline. It 
follows that the entropy S is constant along each 
streamline but generally assumes different values for 
different streamlines. 

(d) The angle ¢ of yaw of the cone is so small that ¢ 
and higher powers of ¢ may be neglected. 

Since the vorticity V-q* and the entropy gradient VS 
are related by 


Af, of the free stream is 


(V-q*)-q* = TVS (1) 


where 7 denotes the absolute temperature, it is seen 
that, by virtue of assumption (c), the flow is necessarily 
rotational. 

Let now X, Y, Z be the components of the velocity q* 
in a system of spherical coordinates 7, 0, w, and let w = 0) 
be the plane of yaw (see Fig. 1). For the case of axial 
flow (« = 0), Taylor and Maccoll’ have shown that 

U(0), Y = V(@) = U'(@), (U' =dU/d0), Z = 9, 


and that U, V satisfy the ordinary differential equation 


(U" + U\(V? 
with 


— a’) = aU + V cot @) (2) 


=(Y)(y-DA-a), w= U+V? @) 
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if the velocity components are measured in terms of the 

speed c of efflux into a vacuum. For later consider- 

ations [see Eq. (16)] we mention here that alwaysV < 0. 
Further, 


So = log Po/ po” (4) 

is constant everywhere. 
Because of assumptions (b) and (d) it is reasonable 
to expect that X, Y, Z, p, p are expressible in the form 


X — U = f(6, w, €) (5) 
Y — V = g(@, w, €), etc. (6) 


where f, g are small functions. Since they must be 
periodic in w with the period 27, f, g can be expanded in 
Fourier series 


(8, , €) = = o lx, (6, €) cos nw + u,(8, €) sin nw], 

uy = 0 
g(0,w,€) = 2 0 Lyn(9, €) cos nw + », (8, €) sin nw], 

w= 0 


here Xn, Mn, Yn» Y are small for small «. The same rea- 


soning leads to 

6 — 6, = Tla,(e) cos nw + A,(e) sin nw}, Ae = 0 
for the equation of the shock wave, #, denoting the 
angle of the shock wave in the nonyaw case and a,, A, 
again being small. 

Since, further, the components of VS in the r, 6, w 
directions are S,(.S, = 0S/0r), (1/r)Sp, , [(1/r) sin 0].Se, 
it follows from assumption (c) that 

0 = g*vS = A oa YS,/r + PAE a sin 6 


Here S, = 0, and Z, S, are small if € is small, so that 
ZS,, is of the same order as e?. The result, then, is 


5S, = 0, i.e., S = S(w, €) 
whence it may be deduced, as above, that 
S — Sy = X{d,(e) cos nw + D, sin nw) (7) 
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Finally, the equation of the cone, by virtue of 
assumption (d), turns out to be 


6 = 0, + €cosw (8) 


Now the flow must be of such a nature as to have a 
vanishing velocity component in the direction normal to 
the surface of the cone. By means of the Taylor ex- 
pansion of Y about 6 = @,:Y = Y, + Y;'(@ — 4), 
Eqs. (8) and (6), this requirement may be seen to result 
in 


Ww = 2eU (9) 
Vn = O, for ~ 1,v, = Ofor all n (10) 
for 6 = @,. As shown in reference 1, these conditions 


(9) and (10), in conjunction with the equation of mo- 
tion, the continuity equation, and the equation of state, 
necessitate the vanishing of all the quantities u,, v,, and 
of Xn, Yn, etc., except for m = 1, leaving as unknowns the 
functions d;(e), ai(e), xi(€, 8), v1, %, and the density and 
pressure perturbations £ and m,. Furthermore, d,(e)/e, 
etc., turn out to be independent of e, if terms of order e’ 
are neglected. With the introduction of d;/e = d, 
a,/€ = a, etc., where d and a@ are now constants, and x, 
y, 2, £, n are functions of @, Eqs. (5), (6), etc., finally be- 
come 


X —-U =e cosa | 
Y — V = eycosa | 
Z=esinw 
p—-p =e cos w } (11) 


Pb — Po = EH COS W | 
S — S, = ed cosa 
5 — 0, = eacosw } 


The last equation, in particular, shows that the dis- 
turbed shock surface is again a circular cone of semi- 
angle @,. The yaw a of the conical shock surface 
i.e., the angle between its symmetry axis and the free 
stream velocity vector g;*—is given by 


a, = €a (12) 


THE DIFFERENTIAL EQUATIONS 


For the determination of the six quantities x, y, 2, &, 
n, d the equations of motion, continuity, and state 
[Eq. (7)] supply these five differential equations: 


y—x =0! 
Vy’ + (U" + U)y + Ve + (n’/po) — 
(Epo' / po) “a | 
V,' + (U + Veot 6)z — [n/(po sin 6) ] = 0; (13) 
y’ + (cot 6+ po’/po)y + 2x + zesc d+ 
V(E/po)’ = 0| 
(n/po) — (vé/po) = d} 


| 
o 


These equations are to be solved subject to the usual 
boundary conditions, which require the fulfillment of the 
Rankine-Hugoniot equations at the shock and the van- 
ishing of the normal velocity component at the cone. 











658 JOURNAL OF THE AERONAUTICAL SCIENCES—NOVEMBER, 


Of the five conditions thus obtainable at the shock, 
we mention here only the one containing a: 


a = —xcot 6/(U + V cot 6) (14) 


System (13) of differential equations is now reducible 
to one differential equation in x(@), as follows: Elimi- 
nation in Eq. (13.2)—by Eq. (13.2) is meant the second 
equation of the system (13)—of y and é by Eqs. (13.1) 
and (13.5) and integration show that 


(15) 


since the constant of integration may be shown to van- 
ish. Now expressing in Eq. (13.3) 7/po by means of 
Eq. (15), we get this differential equation for the 
quantity ¢ = x + zsin@: 


Vt! = Ut + [dpo/(y — 1)po] = 


Vx’ + Ux + (n/po) + dpo/(y — 1)po = 0 


Its solution is 


d : pod8 
+= 1 ( Po Sin ) be i Voo(— Vpo sin 8) /: 
(16) 


Finally, insertion of Eqs. (13.1), (13.5), (15), and (16) 
in Eq. (13.4) lead to the following equation: 


(x/d)"” + A(x/d)’ + B(x/d) + C = 





(17) 
where * 


A = coté@+ AI2U + 3V cot é@+ 
(y + LIAV(U + V cot 4)] 


B = 1 — cot? @+ A[—V cot? 6+ 
(y — 1)AU(U + V cot 8)] 
C = (1 + AV)(t/d) csc? 6 
and \ = V/E, E = a* — V*. The conditions at the 


shock may be condensed into 


x —a’ 


lee y(y —1)(U+ V cot 6) 
(3) cot 0(2U — V cot @) 
.) a 


y(y? — 1)(U + V cot 6)? 
If we introduce the expression for \ into A, B, C, we 
obtain 


A,=AE = a’ cot 0 + [2a? + (vy — 1)V?] X 
(V/E) (U + V cot @) 





(18) 





B,= BE = —a’* cot? 06+ E+ (vy — 1) X (19) 
(UV?/E) (U + V cot 6) 
C, = CE = a*(t/d) csc? 0 
and 
E(x/d)” + Ax(x/d)’ + By(x/d) + Ci = 0 (20) 


for Eq. (17). 
Once this equation has been integrated, the quantity 
d can be calculated from Eq. (9) written in the form 


1948 


d(x/d),’"= 2U, (21) 


the subscript s indicating evaluation at 6 = 6@,; then 
2, n/po, ~/po may be determined as functions of # by 
means of Eqs. (19.3), (15), and (13.5). We shall be 
mainly concerned here with the quantities a and 4,; 
these can be obtained from Eqs. (14) and (15). 


THE NORMAL FORCE ON THE CONE 
The normal force N on a finite cone is given by 


= Sp, cos 0, cos y do (22) 


where y is the azimuthal angle measured from the plane 
of yaw around the cone axis (see Fig. 1), and the in- 
tegration is to be carried out over the surface o of the 
For a 


cone. Now by Eg. (11.5), p = po + en cos w. 
point on the cone, then, 

Ds _ (po)s + €”s cos w 
or 

Ps = (fo)s + ens cos (23) 
since cos w = cos y to the accuracy considered here. 
Introducing Eq. (23) into Eq. (22) leads to 

= (1/s)rD? cot 6,€n; (24) 


D being the diameter of the base of the cone. 
The normal force coefficient Ky is defined by 


Ky = N[piq:*D? sin €]— 
By Eqs. (24) and (15), then, 


0 t 4; 
Ky = —2 Sty, of 1+ +- 
Spa gn 





— U,?d = [es 
27U, 


The ratio R = (po);/p: occurring in Eq. (25) can be com- 
puted from 
R= ((p0)s/ Pw) (Pw/ pr) 
(p0)s/Po = [(1 — U,?)/(1 — qu?) ]VO- 
and 
Po  __— (y + 1)Mi’ sin? 6, _ 
a (y¥— 1) M,? sin? 6, + 2 





SAUER’S METHOD OF POTENTIAL FLow 


As already pointed out in the introduction, the 
method developed by R. Sauer differs from the one 
described above mainly in that the former is based on 
the hypothesis of irrotationality, and, hence, isentropy 
of the flow. Another distinctive feature of Sauer’s 
approach consists in assuming the same shock surface 
for both axial and oblique flow and rotating the cone asa 
whole through a small angle in order to insure the van- 
ishing of the normal velocity component of the oblique 
flow at the cone. The other assumptions made by 
Stone—namely, steadiness of the flow, its conical char- 
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acter, and smallness of the perturbations—are left un- 
changed. 
Starting, then, from 


P(r, 0, w) = o(r, 0) + H(r, 0) cos w (26) 
where ® denotes the potential function of the oblique 
flow (again in spherical coordinates), g the potential 
function of the axially symmetric flow, and H the small 
disturbance potential; noting that ¢ = rU(@) and that 
His of the form H = rh(@), there result for the velocity 
components X, Y, Z the expressions 


X =%,= U+hcosw 
Y = @/r= V+h' cosw 
Z = ©,/rsin@ = —(sin w/sin 4h 


(27) 


If these relationships are introduced in the continuity 
equation V-(pg*) = 0, if use is made of the equation of 
state, and if terms of the order h* are neglected, the 
following equation is obtained: 


Eh” + A,h'’ + Bh = 0 (28) 
where E, A;, B, were defined in Eqs. (19). It is thus 
seen that, with respect to form, Eqs. (20) and (28) differ 
only in the term C, containing the variation of entropy. 

Now let 6 be the small angle of attack (see Fig. 2). 
As shown in reference 2 the Rankine-Hugoniot condi- 
tions specialize to 


h/8 = U tan 0 
(h/B)’ = [2(¢y — 1)U/(y + 1)] + Veot@ (29) 
at 0 = @,, to the order of accuracy considered here. 


Once Eq. (28) for the function ¢ = h/8 subject to the 
conditions (29) has been integrated, it becomes possible 
to compute the angle v through which the cone must be 
turned until the normal component d@/dn of g* = V® 
This requirement may be shown to lead to 


y = h’/2U, (30) 


vanishes. 


at 6 = 0, 





SHOCK WAVE 








SHOCK WAVE 


Fic. 2. Angles used in Sauer’s theory. 
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Since 8 here represents the shock yaw and 8 — » is the 
cone yaw, the ratio a corresponding to a = a,/e [see 
Eq. (12)] can be obtained from 


a = 





(31) 


where, by Eq. (30), 


THE NORMAL Force COEFFICIENT Ky 


Let us now calculate the normal force coefficient Ky 
based on the pressure distribution resulting from 
Sauer’s theory. 


Since 
yb/p = (*/2)(y — 1) — 9’) 
so that 
b = ('/ev)(y — lel — g*) 
and 
p, = F(1 — q,2)/- 
with 


F = (1 — qu®)-“7~ » (pu/pr)pr 
it follows that 
bs = (‘/ey)(y — 1)F(L — 92)” ~ 9 
But by Eqs. (27), 
qg.”° = X,?+ Y,2+ Z,? = U,? + 2U,8f, cos w 


Therefore, using the binomial theorem, and, again, 
cos w = cos yf, 


b: = (C/ey)(y — 1)F(1 — U,2)4% ~ P 
(po)sUs 855 cos y 


By definition [Eq. (22)], then, 


NO = — (?/s)4D? cot 9,(po)s U,8f; 
whence, by Eq. (31), 
Kyw = (1/8)R(cot 6,/q7) Usa (32) 


NUMERICAL RESULTS 


Integrations of Eq. (28) were carried out for a number 
of cases, and the values of a and Ky were deter- 
mined. Some of the results are shown in Table 1. 
The values of a exhibited there are taken from the 
M.I.T. Table No. 3. 

As may be deduced from the cases shown in this 
table, the values of a and a show good agreement, the 
error Aa = (a — a )/a being in no case larger than 
about 3 per cent. 

In view of this result it is permissible, for the following 
considerations, to put 


(33) 


a = @i® 
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TABLE 1 
0, M, a a Aa(%) 
9° 32’ 1.678 0.112 0.111 0.9 
10° ° 2.779 0.341 0.339 * 0:8 
20° 2.839 0.729 0.715 1.9 
20° 4.154 0.864 0.847 2.0 
30° 1.820 0.800 0.775 3.1 
i re) 1.055 1.041 1.3 





Now a comparison of Eqs. (11.1) and (27.1) shows that 
ex, = —£f,. Since the angle 8 corresponds to the 
shock yaw ay, it follows that 


x3 = —(B/e)§s = —(a1/€)f = —af, = — al, 


in view of Eq. (33). By virtue of this fact Eq. (32) 


becomes 
Ky = —(x/8)R(cot 0,/qi?) Usxs (34) 


Ky is thus seen to be the main term in Eq. (25) for 
Ky. We may conclude, then, that 
wr cot é,1 — U,? 
Ky — Ky = —- R—— ——d (35) 
8 qi? 27 
In case of irrotationality, and, hence, isentropy, d = 0 
in the Eq. (11.6) of state, so that, by Eq. (35), the 
values of Ky and Ky agree. For d > 0, it may be 
deduced from Eq. (35) that Ky < Ky. The error of 
Ky in that case is given by 


Ky — Ky 27U, *) si 
i Ge fe Ee l 
mle uae) +] 


(36) 


As the M.I.T. table shows, the quantity 2yU,x,/d(1 — 
U,*) increases for increasing 6,, as well as for increasing 
M, from large negative values to a value of about — 1.4. 
The error AK may thus become large. This fact has 
also been illustrated in Fig. 3, where AK has been 
plotted against Mj, for some values of 6,. If one were 
to take into account the actually existing small dis- 
crepancy in a — a, which, by Table 1, seems to be 
positive, then the values of AK would have to be multi- 
plied by 1 — Aa + (Aa/ AK); since Aa > 0, AK < 0, 





100>——F - - 


Beueunns 
boca 














FREE STREAM MACH NUMBER M, 


Fic. 3. Error of normal force coefficient due to neglect of 


vorticity. 
one would thus obtain values of the error in Ky some- 
what below those depicted in Fig. 3. The dashed 
curve in Fig. 3 gives the corrected error for @, = 20°, 
on the assumption of Aa = 2 per cent. 
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The Bending Modulus of Rupture of Round 
~ Magnesium Tubing. 


C. H. MORTENSON* 
The Dow Chemical Company 


SUMMARY 


Bend tests were conducted on 28 round tubes of FS-la and Ma 
magnesium alloy. Each alloy covered the d/t range from 20 to 
75. The results are summarized empirically, and the following 
formula is established as representing a minimum value for the 
bending strength. 


fr (min.) = 0.7 F.y(Et/ Feyd)*-* 
SYMBOLS 
f; = modulus of rupture stress 
F.y = compression yield stress 
E = modulus of elasticity 
t = wall thickness 
d@ = mean diameter 
ors = fr/Fey 
1/5s = Et/F.yd 


INTRODUCTION 

peer OF THE LACK of any previously published 

experimental data, these bending tests were con- 
ducted to establish preliminary design data for the 
bending modulus of rupture of round magnesium tub- 
ing. In order most effectively to simulate actual pro- 
duction applications and at the same time cover a 
large d/t range, these round tubes were fabricated by 
arc welding from flat magnesium sheet. The method of 
testing and evaluation of test results is similar to that 
of reference 1 for aluminum and steel. The selection 
of an empirical formula to fit magnesium data differed 
only slightly from that used in reference 1. 


PREPARATION OF SPECIMENS 


A total of 28 specimens was prepared for this test, 
three of which were extruded tubes, while the remaining 
25 were fabricated by welding from flat magnesium 
sheet. The fabricated tubes were hot-rolled from flat 
Ma and FS-1la sheet alloy and closed by are welding a 
longitudinal butt weld. The tubes were rolled so that 
the longitudinal direction of the tube was transverse 
to the rolling direction of the original sheet. All tubes 
were welded in a fixture that fitted the outside diameter. 
This resulted in the O.D. being reasonably constant 
over the thickness range of 0.040 to 0.125 in. All 
fabricated tubes were prepared with a nominal outside 
diameter of 3 in. and a length of 36 in. 
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METHOD OF TESTING 


The test setup is shown in Fig. 4. Variations were 
made in the distance between the load support points 
during the tests by increasing the length of the tubes. 
This was accomplished by the use of solid extension end 
plugs, which were used in all tubes with a d/t < 36. 
The exact distances are recorded with the test data in 
Table | and vary between 4.6D and 7.4D. The smaller 
distance (4.6D) was suitable for the larger d/t ratios 
where the bending failure was local and occurred cen- 
trally between loading points. The greater distance 
(7.4D) was provided when it was observed that the 
specimens having the smaller values of d/t were failing 
by a gradual flattening of the tube between loading 
points in addition to extreme deflection. The greater 
span between loading points permitted the tubes to 
assume their natural shape without restraint from the 
semicircular formers at the loading points. 

The extruded tubing was loaded through snug-fitting 
wooden formers in the loading hangers. The fabri- 
cated tubes were loaded through steel formers over the 
upper half of the diameter only. The width of both the 
wood and steel formers at the point where they con- 
tacted the specimens was */s in. Solid circular formers 
about | in. long and snugly fitting the inside of the tube 
were located under the loading formers in all tubes 
witha d/t > 47. 

The end support formers rotated on knife edges, 
which were, in turn, mounted on small roller bearings 
not visible in the figure. This practically eliminated 
the possibility of any axial loads in the tubes. The ex- 
truded tubing was tested with the thinnest wall in com- 
pression. The fabricated tubing was tested with the 
weld seam adjacent to the neutral axis on the compres- 
sion side. 


RESULTS 


Test results are given in Table 1. The design curves 
of the bending modulus of rupture as determined em- 
pirically from the test data, are plotted in Fig. } as a 
function of d/t for various values of Fy. 

In order to show the fit of various empirical curves 
to a nondimensional plot of the test data, Fig. 2 has 
been included. The relationship between two mean 
hyberbolic curves of the same type used in reference 1 
for steel and aluminum, as well as a mean power func- 
tion curve, is shown. A minimum value of the power 

















662 JOURNAL OF 


THE AERONAUTICAL SCIENCES—NOVEMBER, 


TABLE 1 
Tubing Bend Test Results 
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Ult. Maximum 
Mach. Moment 
P Si Steax: L Load Supported fr = MC/I l/o = Et ors =f, 
Specimen, In. d/t (In.) (In.) (In.) (Lbs.) (In. Lbs.) (Lbs. per Sq.In.) Feydm Fey Fey 
FS-la Alloy : 
0.041 K 3.05 74.4 14.15 9.81 33.77 1,235 6,930 23,800 4.66 1,253 19,000 
0.042 X 3.06 72.8 14.05 10.57 35.00 1,085 6,680 22,150 4.8 1.175 18,850 
0.052 * 3.04 58.5 14.06 10.74 35.18 1,540 9,230 25,300 5.31 1.188 21,300 
0.052 K 3.08 59.2 14.10 10.64 35.34 1,615 9,540 25,600 5.18 1.186 21,600 
0.067 X 3.07 45.8 14.09 10.74 35.45 2,050 11,980 25,470 6.75 1.186 21,500 
0.081 XK 3.06 37.8 14.04 16.25 46.42 1,540 13,980 25,350 8.31 1.190 21,300 
0.103 K 3.04 29.5 13.97 16.21 46.40 2,050 18,080 26,900 10.80 1.275 21,100 
0.100 X 3.02 30.2 14.07 16.47 46.92 2,065 18,450 28,080 11.10 1.396 20,100 
0.100 K 3.05 30.5 13.99 16.37 46.70 1,940 17,330 26,000 11.25 1.326 19,600 
0.128 X 3.02 23.6 20.15 15.79 51.65 2,925 24,500 30,400 13.33 1.408 21,600 
0.128 XK 3.06 23.9 20.10 15.74 51.40 3,025 25,200 30,400 12.56 1.345 22,600 
FS-1 Ext. 
0.125 X 2.5 20.0 18:84 15.72 48.83 2,250 19,100 36,300 15.08 1.512 
FS-1 Ext. 
0.125 K 2.5 20.0 18.95 15.48 48.66 2,125 17,800 33,900 15.08 1.412 
Ma Alloy 
0.040 K 3.05 76.2 14.09 9.98 34.06 770 4,730 16,840 7.85 1.531 11,000 
0.054 X 3.03 56.2: 14.09 10.70 35.42 935 5,950 16,050 11.00 1.493 10,750 
0.052 X 3.04 58.5 14.09 10.68 35.40 1,110 6,870 19,350 7.20 1.232 15,700 
0.053 & 3.05 57.6 14.12 10.62 35.28 985 6,170 16,910 9.93 1.458 11,600 
0.065 K 3.02 46.5 14.07 10.74 35.49 1,175 7,270 16,700 12.67 1.480 11,300 
0.065 X 3.05 46.9 14.24 10.68 35.53 1,225 7,480 16,680 12.53 1.477 11,300 
0.082 X 3.01 $6.7 14.08 16.37 66.73 1,035 9,930 18,800 15.43 1.595 11,800 
0.082 X 3.03 37.0 14.07. 10.7 35.51 1,700 10,100 18,730 15.45 1.600 11,700 
0.081 X 3.01 37.2 14.05 18.45 50.78 900 9,940 18,410 14.80 1.513 12,190 
0.102 K 3.00 29.4 14.08 16.35 46.75 1,350 12,500 19,300 21.00 1.772 10,900 
0.100 X 3.07 30.7 14.13 16.42 46.93 1,300 12,140 18,050 20.35 1.680 10,750 
0.100 X 3.07 30.7 14.11 16.40 46.87 1,200 11,300 16,680 19.90 1.516 11,000 
0.123 K 3.03 24.6 14.05 16.438 46.78 1,850 16,680 21,200 20.80 1.606 13,200 
0.123 X 3.06 24.9 20.05 15.62 52.39 1,925 16,440 20,600 20.50 1.550 13,300 
M Ext. 
0.12 5X2 2 5 20.0 18.90 15.59 48.41 1,700 14,630 27,800 © 27.00 2.320 
| i | | | | | rca. ti. . tl wee ce pe ge ag ee ee ae ee 
7] sh E a "7 
<= Frnin = TFey(Et a] cen 20 F ‘ ai 
os a a | \-238 
& Org =a(1 J o 7 
bs 40 — — ave. 8s 0 ae 
6 =s=8---- 4 
a oO -3=- —— - oe 
S - _ o =o o 
p Fey - ose 4 
as 30 ‘pepe <a Y 2” AD 
2 30,000 & 12 o, =7( a —_ 
= — Ss J 7 
S 20 |— 20,000 + x en oe 
=. "Lae erected MOOI 
8 10-— 7 t Wy = L75x_ ° FS-la 
= 6 aH / X+2.63 a 
"te We Me are ae a a ee OS a ae a : ere 
10 20 30 40 50 sO.) CD 80 ne I. rel % a: ore = 
(0) 1 
% fe) 4 : 12 16 20 24 
Fic. 1. The bending modulus of rupture of round magnesium - 5 
alloy tubing. Ss 
Fic. 2. The comparison of several empirical curves fitted to 


function curve is also plotted. The equation of this 
curve was used to establish the design curves of Fig. 1. 

An illustration of typical load deflection curves for 
the FS-la tubes with four values of d/t is shown in Fig. 
3. The significance of these curves will be brought out 
in the section titled “‘ Discussion.” 

In order to differentiate between the types of failure 
in the large and small d/t ranges, Fig. 5 has been in- 
cluded. The gradual bowing of the small d/t specimens 


the nondimensional data for FS-la and Ma alloy. 


is compared to the crippling type of failure in the walls 
of those tubes with a large d/t ratio. 


DISCUSSION 


In the analysis of the test data contained in Table 1, 
the modulus of rupture, f,, has been computed from the 
flexure formula: 
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Fic. 3. Load deflection curves for FS-la alloy tubing. 
O.D. = 3 in. 
f, = Me/I 


where M is the maximum moment supported by the 
specimen and J/C is the section modulus. The test 
results contained in Table 1 have been converted to 
nondimensional form and plotted in Fig. 2 in a manner 
similar to that of reference 1. The plot of the data in 
Fig. 2 makes use of the following parameters: 


ors = fr/ Fey 


and 
1/3, = Et/ Fade 
where 
t = wall thickness 
dm = O.D. — t (mean diameter) 
E = modulus of elasticity (6.5 X 10° Ibs. per 
sq.in.) 


Fy = compression yield strength 


The compression yield strength shown in Table 1 
was obtained from flattened specimens cut from un- 
stressed portions of the tubes after static test. The 
properties were taken in the longitudinal direction of 
the tube, which is transverse to the rolling direction of 
the original sheet. 

Three empirical curves have been fitted to the non- 
dimensional plot of the data, as shown in Fig. 2. Two 
of these curves are hyperbolas of the form y = ax/- 
(x + b) and y = a — (b/x). They were selected for 
fitting as a result of their use in the interpretation of 
the data for aluminum and steel tubing in reference 1. 
In addition a simple power function curve of the form 


y = ax” has been fitted to the data (Fig. 2) and has 
been selected to compute the design values of the 
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modulus of rupture stress as plotted in Fig. 1. This 
formula has been fitted to the data twice, once as an 
average curve through the test points and again as a 
minimum curve representing the test data and sug- 
gested as a minimum design value. 

The construction curves for the determination of the 
constants in hyperbolic and power function formulas 
will not be discussed here, since the method is simple and 
straightforward. It is apparent that the two hyperbolic 
curves fit the data about as well as the power function 
curve. The latter was selected to represent the data, 
because of its simplicity. 

The use of the results in Figs. 1 and 2 should not be 
without careful consideration of the permissible 
amount of deflection. This is particularly true when de- 
signing in the range of d/t < 35. A study of the load 
deflection curves plotted for each tube indicated that, 
for a d/t ratio of less than about 35, the deflection will 
probably become critical before the tube actually 





Fic. 4. Test setup for determining modulus of rupture. 





Typical failures of FS-1la alloy tubing. 


Fic. 5. 
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fails by bending. Fig. 3 has been included to show 
the load deflection curves for several values of d/t of 
the FS-la tubes. It is apparent in these curves that, 
for a ratio of d/t < 38, the last 10 per cent of load is ac- 
companied by approximately a 100 per cent increase 
in deflection. Obviously, these deflections are so ex- 
cessive as to be intolerable in any structure. In gen- 
eral, those tubes with a d/t > 38 failed by sudden in- 
stability of the tube walls, which produced a sharp 
ended load deflection curve similar to tube No. 5, Fig. 
3. The failures of those tubes with a d/t < 38 con- 
sisted of a general flattening of the tubes together with 
excessive deflection. Fig. 5 clearly illustrates the type 
of failures experienced in the FS-la tubes over the en- 
tire d/t range. 
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CONCLUSIONS 
The bending modulus of rupture of round magne. 
sium FS-la and Ma alloy tubing can be represented by 
the following formulas: 


f-(min.) = 0.7F,,(Et/ Fadm):*® 
flav.) = 0.8Fey(Et/ Feydm) 9-238 


The first formula represents a minimum value of the 
test data and has been used to establish the design 
curves of Fig. 1. The second formula represents a 
mean curve through the test points and may also be 
used for design purposes if the designer so wishes. 
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Letter to the Editor 


Dear Sir: 

In his ‘‘Letter to the Editor” (JoURNAL OF THE AERONAUTICAL 
Scrences, Vol. 15, No. 8, p. 502, August, 1948), F. R. Shanley 
has discussed some results in my paper “Inelastic Column 
Theories and an Analysis of Experimental Observations’”’ (Jour- 
NAL OF THE AERONAUTICAL SCIENCES, Vol. 15, No. 5, p. 283, May, 
1948). However, I feel that the points raised in the discussion, 
need further clarification. 

It would be indeed a mistake to use the coefficient of end fixity 
C in the inelastic column theory as a “multiplying factor.” In 
Eqs. (16), (18), and (19) of my paper, C was, however, not used as 
a “multiplying factor.’’ This misunderstanding can be attri- 
buted to the fact that in the well-known elastic case—i.e., in the 
case of Euler’s formula—C can actually be used as a ‘‘multiplying 
In the inelastic case, the Eq. (16) 


Per = (Cw*EI)/L? 


factor.” 


is, in fact, not sufficient in determining the buckling load P,,, 
since both P,, and EI are variables. Actually, one additional 
equation is needed to give the relationship between P,; and EI, 
and this is usually tacitly given in the form of stress-strain curve 
of the material from which this relationship can be computed. 


Poy, therefore, must be determined from two simultaneous equa 
tions, one of which is given in the form of a curve. If this faét is 
made clear, there is no danger of misinterpreting C as a “‘multi- 
plying factor’ as mentioned by Shanley. Shanley has also sug- 
gested the use of another factor ko = Lo/L; actually this amounts 
exactly to the same thing, for ko is equal to VC. 

Mr. Shanley has pointed out that the use of two tangent 
moduli instead of a simple one in the writer’s paper does not lead 
to simple computation. This is indeed so, since the two-moduli 
formula was never intended for actual computational use. The 
reason for introducing the two moduli is to explain some of the 
column behaviors that could not be explained by the tangent 
modulus theory. Mr. Shanley has contended that the deviation 
of test results from the tangent modulus formula can all be 
attributed to the initial imperfections of the column specimens. 
However, in the writer’s paper, in using the two-moduli formula, 
one can go one step farther by proving that, even if perfect 
column specimens are possible, there would still be deviations 
from the tangent modulus formula. 

Cui-TEH WANG 
Associate Professor of Aeronautical Engineering 
New York University 
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Unsteady One-Dimensional Flows with Heat 
Addition or Entropy Gradiénts’ 


A. KAHANE? ano LESTER LEES}? 


Princeton Unwersity 


SUMMARY 


The present interest in aerodynamic propulsive devices, such as 
the pulse-jet and the ram-jet, has stimulated recent research in 
the field of combustion. One important aspect of the combustion 
problem considered herein is the unsteady gas dynamic effects 
associated with the burning process. 

Nonlinear differential equations of the Riemann type are de- 
rived for the solution of problems involving the propagation of 
one-dimensional waves in flows in tubes of slowly varying cross 
section with heat addition or entropy variation. These equations 
can be solved graphically or numerically with the method of 
characteristics; the method is described in the Appendix. The 
transient flows arising when heat is added to a section (as in a 
combustion chamber) of an initially isentropic flow in a tube have 
been calculated. The results of this calculation afford an insight 
into the gas dynamic aspects of intermittent heat addition in a 
flowing gas, as well as to the apparently anomalous behavior at 
sonic velocity of steady gas flows with heat addition. The equa- 
tions derived are well suited to the calculation of the nonlinear 
pulse-jet cycle with heat addition, once a suitable model of com- 
bustion has been selected. 


INTRODUCTION 


— PROPAGATION OF PLANE WAVES due to combus- 
tion was theoretically investigated by MacDonald’ 
and others with application to the pulse-jet cycle. 
The Lagrangian form of the equations of motion were 
used; analytic solutions were obtained by linearizing 
through the use of the small perturbation method. A 
disadvantage of the linearization, however, is that waves 
of finite ‘amplitude are essentially nonlinear and, with 
linearization, phenomena such as the distortion of the 
wave form and the formation of shock waves are pre- 
cluded. In the present report, nonlinear differential 
equations of the Riemann type are derived for the solu- 
tion of problems involving the propagation of one- 
dimensional waves in flows in tubes of slowly varying 
cross-sectional area with heat addition or entropy gra- 
dients. These equations can be solved graphically or 
numerically with the method of characteristics; the 
method is described in the Appendix. The transient 
flows arising when heat is added to a section (as in a 
combustion chamber) of an initially isentropic flow in a 
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tube have been calculated. The results of this calcula- 
tion afford an insight into the gas dynamic aspects of 
intermittent heat addition in a flowing gas as well as to 
the mechanism of the readjustment of the steady flow 
in a ram-jet when the quantity of heat added to the flow 
is changed. 

The most complete analysis of the pulse-jet cycle up 
to the present time has been made by Schultz-Grunow.? 
His analysis, however, does not consider the effect of 
heat input, the combustion process being simulated 
merely by an isentropic pressure rise taking place in the 
combustion chamber, with all the resulting wave propa- 
gation processes also considered as isentropic. Al- 
though these simplifying assumptions lead to a quick 
method of calculating the effect of various parameters 
on the pulse-jet cycle, the method obviously cannot give 
any insight into the effects of heat addition on the 
cycle, nor can values of thrust be calculated. The 
equations derived herein make it possible to calculate 
the nonlinear pulse-jet cycle with a full consideration of 
heat addition, and it should be possible to obtain values 
of thrust from such calculations. 


SYMBOLS 


a = velocity of sound, ft. per sec. 

= cross-sectional area of tube, sq.ft. 

= specific heat at constant pressure, B.t.u. per slug per °F. 
specific heat at constant volume, B.t.u. per slug per °F. 
length of heating chamber, ft 

static pressure, lbs. per sq.ft. 

Riemann variable, P = u + [2/(y — 1)]a, ft. per sec. 
= Riemann variable, Q = u — [2/(y — 1)]a, ft. per sec. 
= quantity of heat added per unit mass, B.t.u. per slug 
gas constant, B.t.u. per slug per °F. 

entropy, B.t.u. per slug per °F. 

time, sec. 

absolute temperature, °R. 

= particle velocity, ft. per sec. 

= distance along tube, ft. 

= ratio of specific heats, C,/C, 

= mass density of fluid, slugs per cu.ft. 
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DERIVATION OF CHARACTERISTIC EQUATIONS 


The characteristic equations for the propagation of 
waves in flows with heat addition or entropy variations 
will now be derived. The flows will be assumed to take 
place in tubes of slowly varying cross section and will 
therefore be assumed one-dimensional. It will further 
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be assumed that the gas constant, R, and the specific 
heats, C, and C,, are constants. 
The equations of motiof are: 


Dynamical equation 


Ou Ou  10p 
= ant ee Ow 1 
ry + u ae + ee 0 (1) 
Continuity equation 
Op . 1 O(puA) 
See oft co ee” ie 2 
ot A Ox . @) 
Energy equation 
dQ* 1S dT d, 
dQ* _ OF i OM on. (3) 
dt dt dt pdt 
Equation of state 
b = pRT (4) 


From Eqs. (3) and (4) the relation between the pressure 
and density can be found to be 


p = constant X pve?’ (5) 


and by definingt 





a = V yp/p 
Eqs. (1) and (2) can be transformed to 
Ou Ou 2 Oa a? oS 
= aon ll woman il ane ss eopeecenenioms ty 1 
” baat cig BS inet (y — 1)C, ox ed 
2 oa 2 Oa Ou y ads 
y —10 ga ae? ie y¥-—1C,d 
uadA 
aaa 2 
A dx > 


where 
dS/dt = (0S/dt) + u(dS/dx) 


is the rate of change of entropy following a fluid particle. 
Addition and subtraction of the above equations give: 


ra) ( a 2a )+ aynere re) ( a 2 ) 
— — a) — —a) — 
ot . y-1 . Ox ¥ y-1 


a? oS y @ds 





y¥—-DGe 7-168. A a = 
The above equations are essentially wave equations and 
indicate that disturbances are propagated downstream 
and upstream at the velocity a relative to the stream 
velocity (the upper and lower signs denote downstream 
and ypstream wave motions, respectively). 

It will be noted that for isentropic waves in constant 
cross-section tubes, the equations reduce to that of 





t The quantity a is equal to the value of the velocity of sound if 
the compression and expansion processes occurring in the sound 
wave are isentropic. If the passage of the sound wave in itself 
causes heat addition, so that the compressions and expansions of 
the sound wave aré nonadiabatic, then a is not exactly the veloc- 
ity of sound.® 
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Riemann ;*’ > because of heat addition and area change, 
therefore, a single progressive wave will be continually 
transmitted and reflected as it passes down the tube, 
and, of course, its reflections will in turn be transmitted 
and reflected. For steady flow, the equations can be 
reduced to those of Chambré and Lin.® 

Since in most cases analytic solution of Eq. (6) is ex. 
tremely difficult, the numerical or graphical calculation 
procedures of the well-known characteristic method 
must be used. Such calculations are carried out on the 
time-distance or t-x plane. The equations of the char- 
acteristic lines are 


(dx /dt); 
(dx /dt)» 


ut+a 


u—a 


If the variables P and Q are introduced (Riemann and 
others used the symbols ‘‘r” and “‘s’’), such that 


P=u+ [2/(y — l1)Ja 
Q=u— [2/(y — Je 


then the total change of these variables with time as one 
proceeds along characteristic lines is 


6P oP 4 of (=) 

dt ot ow \dt/, 

60 a2, a (2) 

dt ot em\dt/, 

and from Eq. (6) 

— 5 I dy Ms 
4 G@=-DGe' +-168 A&ese 
POET A... LER ES I 
dt (y—1)C,d%& y—1C,dt'Adx “ 


Thus for a given tube the characteristic network may 
be constructed as long as the entropy variations are 
known. A method that may be used in constructing 
the network is described in the Appendix. The method 
of characteristics is also described in the Shock Wave 
Manual’ and by Sauer.® 


APPLICATIONS 


The equations derived in the preceding section have 
many interesting applications in the field of gas dy- 
namics. Some possible applications include the calcu- 
lation of the pulse-jet cycle, once a suitable model of 
combustion has been assumed, and studies of the sta- 
bility of the flow pattern in the ram-jet. A study of the 
stability of the ram-jet operation might include the gas- 
dynamic effect of a wave disturbance upon the initial 
steady-flow combustion pattern or the effect of disturb- 
ances generated by “rough burning’ (simulated by 
small variations of the quantity of heat input) on the 
motion of the diffuser shock wave in a supersonic ram- 


t 6/dt indicates the total derivative along a characteristic line. 
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A problem that is treated in this report is the effect 
of addition of heat to an initially isentropic flow in a 
tube. The processes occurring in such a configuration 
gre similar to those that occur during intermittent com= 
pustion in a flowing gas, as well as being similar to those 
that occur during steady-flow combustion when the 
quantity of heat addition is suddenly changed. 


THE ADDITION OF HEAT TO AN INITIALLY ISENTROPIC 
FLow 


Of interest is the transient phenomena occurring 
when heat is suddenly added to an initially isentropic 
fow in a constant area tube of infinite length. Con- 
jder a tube in which the flow initially is at a Mach 
Number of 0.80. Heat will be assumed to be added 
between stations A and B, at a rate such that 


(dQ*/dt)/T = dS/dt = constant = K 


within the heating chamber. The value of dS/dt was 
chosen as a constant to simplify the calculations. 

The transients can now be calculated with the charac- 
teristic method utilizing Eqs. (7a) and (7b). The 
entropy may be determined. in the following manner: 
a plot of the actual particle paths can be constructed 
along with the characteristic calculation on an x-t plane 
the equation of a particle path line is dx/dt = u). The 
entropy at-any point is found by performing the follow- 
ing integration along the particle path 


S—S, = Si ((dQ*/dt)/T)dt 


Once the entropy distribution is known, values of 
)S/dt and 0.S/dx at a given point can be estimated. 

A great simplification in the present calculation can 
be made by the following observations. If the charac- 
teristic x-t plane is divided up into regions bounded by 
the lines x = A and x = B (the entrance and exit of the 
heat chamber) and the particle path lines of the particles 
that are initially atx = A andx = B at the time ¢ = 0 
(see Fig. 1), then, for the prescribed heat addition 
(dQ*/dt)/T = K, it can be ascertained that the entropy 
variations in the various regions are approximately 


. bm = 0 ) 
o0S/dx = 0 | 
1 (§0S/at = K | 
(OS/dx = 0 | (8) 
my J05/at = 0 
(oS/dx = K/u | 
: jdS/dt =0 | 
nN \oS/dx = —K/u J 


The above formulas were obtained by assuming that 
all the particle path lines in the x-t plane deviate only 
slightly from straight lines all parallel to each other. 
Asan example, consider the region numbered II. Since 
the heat is first added at ¢ = 0, at any given constant 
value of ¢ within region II, each particle has increased 
its entropy by the same amount under the assumptions 
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Fic. 1. Regions of the characteristic plane for which approxi- 
mate entropy derivatives are determined. 
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Fic. 2. Characteristic diagram for the addition of heat in a tube 


to an initially isentropic flow of Mp = 0.8, KL/aoCp = 0.36 


made. Thus, clearly, 0S/Ox = 0. Then, since 


dS/dt = (0S/dt) + u(OS/dx) = K 
it follows that 
0S/dot = K 
A similar reasoning applied to the other regions will give 


the remaining equations. 


These approximations to the entropy variations were 
used in calculation of the transient flow. Use of the 
approximation may also be regarded as equivalent to 
slightly varying the addition of heat during the tran- 
sient phase so that the expressions (8) are satisfied. 


The characteristic x-t plane for this calculation is 
shown in Fig. 2; the main results are illustrated in Fig. 
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Distribution of u/ao and a/ao for several values of the 
time parameter dot/L. 


Fic. 3. 


3, in which distributions of u/a) and a/ay are plotted 
against «/Z for several values of the time parameter 
aot/L. These graphs indicate that, with the addition of 
heat, compression waves are propagated both upstream 
and downstream of the heat chamber. It may be 
observed that the upstream moving compression wave 
is rapidly steepening into a shock wave. For the infi- 
nite tube the flow would always be unsteady because of 
the upstream and downstream moving waves. The 
transition to a steady flow in a tube of finite length will 
be discussed in a later paragraph. 

The calculation carried out affords an insight into the 
processes occurring during intermittent combustion in 
an initial steady flow. It is evident that in combustion 
of this type the pressure is raised in the combustion 
chamber, this pressure increase causing the propagation 
in both directions of compression waves. The magni- 
tude of the pressure rise in the combustion chamber is 
proportional to the rapidity of combustion—i.e., the 
rate of release of energy Q*. Thus, for large values of 
Q* the thermodynamic combustion process approaches 
the constant-volume process (9 = constant) obtained 
in the Otto-cycle; whereas for small values of Q*.the 
burning process approaches the so-called constant pres- 
sure process obtained with steady-flow combustion as 
in the ram-jet and turbojet. 


The Transition to the Steady-Flow State 


As pointed oyt earlier in this paper, no steady state 
flow will ever occur in an infinitely long tube when heat 
is being added to an initially isentropic flow. Consider 
now that the tube is of finite length, such as a ram-jet 
moving at a flight Mach Number of M = 0.8. (See 
Fig. 4.) When the heat is added to the flow, the com- 
pression waves that are propagated upstream and 
downstream will be reflected at the entry and exit of the 
tube and will continually reflect and refract until a new 
steady flow is established. The new steady flow will 
depend only on the ram stagnation pressure, the quan- 
tity of heat added by the burner, and the exit pressure 
(assuming frictionless flow) and can easily be calculated 
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by the equations presented by Chambré and Lin’ o¢ 
Hicks.’ Since the exit pressure is atmospheric and since 
the steady addition of heat causes a pressure drop across 
the burner, it is evident that the pressure upstream of 
the burner is higher than that occurring in the com. 
pletely adiabatic flow. Thus, the velocity ahead of the 
burner and, therefore, the mass flow through the system 
are reduced. Also, it may be observed that, since the 
stagnation pressure of the flow is reduced in passing 
through the combustion chamber, the exit Mach Num. 
ber is less than that of the completely adiabatic flow. 
A plot of flow parameters for steady state flows with 
several values of the quantity of heat added is presented 
in Fig. 4. It is evident that increasing the quantity of 
heat added will cause a readjustment of the flow with a 
further reduction of mass flow through the system. 

Consider now a ram-jet propelled at supersonic speed 
with a supersonic diffuser at the inlet (Fig. 5). In this 
configuration the mass flow is fixed once the supersonic 
flow is established in the inlet of the diffuser, and a 
normal shock wave will be present in the diverging por- 
tion of the diffuser. If the quantity of heat added is 
increased, a new steady flow will be established with 
the shock moved to a more forward position. 

The reason for this can easily be seen by considering 
the diverging part of the diffuser as analogous to the 
diverging part of a Laval nozzle. Since the mass flow 
through the system is fixed, as long as sonic velocity 
exists at the throat of the exhaust nozzle, it follows that 
the Mach Number at the burner exit must remain con- 
stant. Thus, increasing the quantity of heat added 
will reduce the Mach Number 1/7, at the entry of the 
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Fic. 4. The steady flow obtained in an initially subsonic flow for 
several values of the quantity of heat added. 
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burner. Since the flow is adiabatic upstream of the 
burner, the mass flow per unit area depends only on the 
stagnation density behind the shock; the stagnation 
temperature, which is constant; and the Mach Number 
at the burner entry—i.e., m = ps VT sf(M,4). There- 
fore, a reduction in 4, must be accompanied by an in- 
crease in ps behind the shock, and the shock wave is 
weaker and occurs farther upstream. Additional in- 
crease in the quantity of heat added will finally move 
the shock to a more forward position, until eventually it 
will reach the converging portion of the diffuser where 
it is unstable’ and it will move out in front of the ram- 
jet tube as a detached shock. Now the internal flow is 
subsonic, and the effect of further addition of heat will 
be the same as that of the subsonic jet. 


The conclusions of the present study also give an in- 
sight into the apparently anomalous behavior at sonic 
velocity of a steady gas flow in a straight tube with a 
steady rate of heat addition. From the equations of 
motion for steady flow it is concluded that the gas can- 
not absorb heat at sonic velocity. Chambré and Lin* 
conclude that a “‘sharp front’’ will result if heat is'added 
to the gas at or beyond the point at which sonic velocity 
is attained. From the present study it appears that 
when a quantity of heat, greater than the “critical’’ 
amount that would bring the gas to sonic velocity, is 
added, compression waves traveling in both directions 
are generated initially, and a new equilibrium flow is 
established. 
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The steady flow obtained in an initially supersonic flow 
for several values of quantity of heat added. 
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Fic. 6. Graphical method of constructing characteristics. 


Appendix 


A GRAPHICAL CALCULATION PROCEDURE 


The method of characteristics is a graphical or numer- 
ical procedure for integrating partial differential equa- 
tions of the hyperbolic type. It is assumed here that 
all quantities including u, a, and S are knowh at the 
points 1 and 2 in the x-/ plane and that the values of 
these quantities are to be found at a third point desig- 
nated 3. The rate of heat release is also assumed known 
throughout the x-/ plane, or 


O* = O*(x,!) 


To calculate the point 3, the characteristic differential 
Eqs. (7a) and (7b) are written as difference equations 
and are solved by substitution in them of mean values 
of the quantities between 1 and 3 and 2 and 3, which are 
referred to as the quantities at I and II, respectively. 
Inasmuch as the quantities at 3 are unknown as yet, 
the computation is of necessity an iterative procedure, 
the first approximation of which is as follows (see 
Fig. 6): 

(1) The characteristic lines are drawn through points 
1 and 2 at angles 


ll 


cot! (14 + a) 
cot! (a2 — uz) 


ay 


a2 


The intersection of these segments will be considered 
the first approximation to the point 3. 

(2) The values of dS/dt and 0S/Odx are to be found 
at the points I and II which bisect the characteristic 


segments 1-3 and 2-3, respectively. The value of 
dS/dt may be found from 
dS (y — 1)C,(dQ*/dt) 





a = (dO*/dt)/T = 


a? 


Since the value of a at I and II is not yet known, the 
values at 1 and 2 may be used. * The value of 0.S/0x 
may be estimated from the definition of a partial de- 
rivative, which for the point I is 
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lim Slr + (Ax/2), 4] — Slr — (Ax/2), t] 
a,-7*0 Ax 





In estimating 0S/dx finite values of Ax must be used. 
The values of S used in the formula are determined by 
integrating along suitable particle path lines ustng the 
equation 


S— Sy) = Si’ ((dQ*/dt)/T]dt 


(3) Eqs. (7a) and (7b) are utilized to determine P3 
and Q;. They may be written as follows: 


P ia | —_— (=) ins Y a (3) a 

SL - DG els ¥—-1 CN, 

1 log A 
Uj; (4) | (ts; — h) 
dx I 
ait =) - arr (3) 

ee, eee emer. eens | no) ee csi == 

Cs Qs E- _ va II oy am 1 C dt Il T 


dlog A 
Uy1Ay1 ( e—) | (ts — te) 
x Jy 


uy = (tt + u3)/2 
Uyy = (ue + U3) /2 . 


where 


ete. Since “3; and a3 are unknown as yet, in the first 
approximation the values of wu and a at 1 and 2 are sub- 
stituted for their values at I and II, respectively. 
Values of u3 and a; are then found from 
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u = (Ps + Qs) /2 
as = [(y — 1)/4](Ps — Qs) 


The second approximation consists of constructing new 
characteristic segments using mean values of u and a 
between 1 and 3 and 2 and 3, finding a better approxi- 
mation to the entropy derivatives, and recalculating P, 
and Q; using the mean values of u and a along the seg- 
ments as well as the more accurate time change. The 
method usually converges rapidly, and for most applica- 
tions the first approximation is sufficiently accurate. 
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SUMMARY 


The problems of avoiding airflow separation, air-frame shaking, 
and control surface buffeting at transonic flight speeds are dis- 
cussed. A design parameter, identified as the angle between the 
aft surface of a body and the line of flight is proposed as a useful 
tool for avoiding these difficulties and is correlated with flight 
and wind-tunnel test data. The correlation showed that, for 
preliminary design purposes, the proposed parameter represents 
a satisfactory index of the Mach Number at which separation 
phenomena will occur on bodies and lifting surfaces. For wing- 
body junctures and tail intersections it is suggested that separa- 
tion Mach Numbers be considered greater than the critical 
Mach Numbers by 0.10. 


INTRODUCTION 


O* HIGH-SPEED AIRCRAFT many new aerodynamic 
phenomena have occurred. Disturbances have 
been encountered by airplanes diving to speeds that 
have never before been reached. In reviewing the 
nature of the disturbances encountered at these high 
speeds, we see that they have taken various forms. 
At first it was usually a longitudinal disturbance as 
evidenced by a change in trim, elevator control position, 
or elevator force, and these phenomena are at present 
reasonably well understood. Such changes in stability 
and control usually resulted from changes in the lift 
curve, pitching moment curve, and down-wash. A 
summary of information on this phase of the problem is 
presented in reference 1. 

More recently, difficulty has been experienced in 
directional stability and lateral control at high speeds. 
The directional difficulty is a type of dynamic instability 
which is sometimes caused by a combination of pressure 
distribution changes and separation on the vertical tail 
which influence the hinge moments of the rudder. The 
lateral control difficulty experienced is a type of con- 
trol surface oscillation, or one-degree flutter, sometimes 
called “‘buzz,’’ a phenomena that has been discussed in 
various papers. These problems are being investigated 
by research organizations, the military services, and 
manufacturers of high-speed aircraft. 

Another problem, which may be related to the above 
difficulties but in some respects is independent, is be- 
ginning to appear. This problem is the ‘occurrence of 
airflow separation and buffeting at high Mach Numbers, 
causing shake of the air frame or buffeting of a control 
surface. It is this problem that is to be discussed, and 
an attempt will be made to provide a means of avoiding 
this type of difficulty on high-speed aircraft. 
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Separation | 


VERNON OUTMAN* anp ARTHUR A. LAMBERTTt 
McDonnell Aircraft Corporation 


THEORETICAL CONSIDERATIONS 


From the practical point of view a design parameter 
is needed which, when correlated with test data, will 
serve as a criterion for avoiding separation and, hence, 
buffeting at transonic speeds. Before discussing a de- 
sign parameter, however, it may be well to examine the 
conditions that may lead to separation at high Mach 
Numbers. An adverse pressure gradient permits sepa- 
ration, and thus the magnitude of such gradients that 
exist on a body will be of cardinal importance in estab- 
lishing a design parameter with which to correlate the 
Mach Number of separation. In this regard, the dif- 
ference in the effect of stream filament contraction on 
the pressure gradients for subsonic and supersonic 
flows must be considered, since both types of flow are 
encountered locally on an airfoil traveling at high 
transonic speeds, and -these effects will have a direct 
bearing on the magnitude of separation that may 
develop. In subsonic flow, a tendency for the stream 
filaments to contract, due to a wedge-shaped boundary- 
layer growth, will result in a tendency to decrease the 
local adverse pressure gradients and thus allay further 
thickening of the boundary layer. In supersonic flow, 
however, the opposite is true. Stream filament contrac- 
tion due to a wedge-shaped boundary-layer growth re- 
sults in increased adverse pressure gradients and, 
therefore, in a greater tendency for the boundary layer 
to thicken or separate. Hence, when a region of super- 
sonic flow exists on an airfoil or similar aerodynamic 
body, disturbances in the boundary layer resulting from 
the occurrence of shock waves, combined with a large 
adverse pressure gradient, will result in a considerable 
boundary-layer growth, possible airflow separation, and, 
in general, a high degree of turbulence. 

As the Mach Number of the principal flow approaches 
unity, the main shock waves move to the trailing edge 
and the unstable tendency of the boundary layer dis- 
appears, as is shown by high-speed schlieren photo- 
graphs. Thus, severe buffeting may be expected to 
begin at free stream Mach Numbers somewhat above 
the critical value for the airfoil and disappear at ap- 
proximately sonic speed. 

High-speed wind-tunnel tests offer a good means of 
indicating regions of transonic separation when the 
model is properly tufted. Usually, however, several 
months are required to build a high-speed model, since 
in most cases they must be made of steel. Therefore in 
the interests of expediency, it is usually necessary to 
settle the design before the wind-tunnel test results are 
available. Because of this the designer needs a tool 
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to lay out the airplane to avoid transonic separation. 
Such a rough tool is indicated by the above rudimentary 
considerations. 


PROPOSED DESIGN PARAMETER 


The single design parameter that offers the best 
index of the adverse pressure gradient existing on a 
body in potential flow is the slope of the surfaces aft of 
the position of maximum thickness. For a given type 
of body, this afterbody angle can also serve as an index 
of the body’s critical Mach Number. Because of this 
interrelation, a design parameter identified as the angle 
between the surface of the aft portion of a body and the 
line of flight is suggested. 

The relation between the suggested design parameter 
and boundary-layer conditions at high Mach numbers is 
illustrated by the pressure distribution diagrams pre- 
sented in Figs. 1 and 2. Fig. 1 presents the variation 
of surface pressure coefficient with Mach Number and 
chordwise position for the 66,2-215 (a = 0.6) airfoil 
which has an afterbody angle at zero lift (@) of 10.7° 
as measured on the upper surface between the 70 per 
cent chord point and the trailing edge. Referring to 
the figure, it is found that the upper surface pressures 
over the aft 30 per cent of the chord are essentially 
constant up to and including a Mach Number of 0.752. 
At M = 0.779 and above, drastic changes in the pres- 

,Sure distribution over the afterbody occur which are 
not explained by potential flow theory and, hence, can 
be attributed to boundary-layer conditions. There- 
fore, the possibility of considerable turbulence of the 


airflow over the upper surface and airflow separation at 
these Mach Numbers is indicated. The experimental 
critical Mach Number for the airfoil at this angle of 
attack is M = 0.67. 

Fig. 2 presents*similar data for the 0015 airfoil where 
6 = 8.8°. For this airfoil, the upper surface pressures 
over the aft 30 per cent of the chord are essentially 
constant up to and including M = 0.805. At M = 
0.837 the aft pressures undergo a considerable change. 
The critical Mach Number for this condition is 
M = 0.70. 

Buffeting and airflow separation is related to these 
pressure distribution changes as will be shown by a 
correlation with the angle @ of the Mach Number for 
airflow separation as determined from flight and wind- 
tunnel test data. 


FLIGHT AND WIND-TUNNEL EXPERIENCE WITH 
TRANSONIC SEPARATION 


Various examples of transonic separation have been 
accumulated and the data correlated in Fig. 3. The 
phenomena, as observed on several airplanes in flight 
and some wind-tunnel models, are discussed below. 


(A) From Flight Tests 

This airplane has a local section which had a rather 
high thickness ratio and the aft portion of the wing 
was steep. At a Mach Number of 0.63 the airplane 
experienced severe buffeting, which restricted the pilot 
from going to a higher Mach Number. Examination 
of low-speed wind-tunnel data showed that, in the 
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region of the wing described above, there was one tuft 
that was slightly disturbed.- Tufts were placed on the 
airplane in the region in question and observed from 
an accompanying “hase” plane. The region was 
found to be much more violently disturbed than indi- 
cated by the low-speed model tests. The fairing was 
installed which reduced the slope of the aft portion of 
the wing in this region considerably. Most of the fair- 
ing was placed on the top surface with the intent to 
take out some of the camber in this region. The 
thickness ratio was changed from 22 to 20 per cent, 
but the slope of the top surface relative to the wing 
was reduced from 24'/,° to 9°. These two points are 
located in Fig. 3. 


(B) From Flight Tests 

This airplane, with a fairly thick wing, experienced 
air-frame shaking, changes in‘lateral trim, and rolling 
instability, which prevented the pilot from going over 
a Mach Number of 0.75. Again, tufts were placed on 
the wing to discover the critical region of flow separa- 
tion, and a fairing was installed to improve this region. 
In this case it was not possible to make an extensive 
change without rebuilding the wing, and only a slight 
gain was realized as indicated on Fig. 3. 


(C) From Flight and Wind-Tunnel Tests 


This airplane began to experience buffeting at a 
Mach Number of 0.70, and the pilot was not able to 
go beyond a Mach Number of 0.75. Tests were made 
in the wind tunnel to determine thé critical region, 
and two modifications were made. The first consisted 
of extending the chord in the critical region, reducing 
the thickness ratio, and reducing the trailing edge angle 
from 21° to 15°. Another modification consisted of 
adding a fairing mostly to the top surface of the wing, 
giving the section a reflex. This modification also re- 
duced the thickness ratio somewhat and reduced the 
trailing-edge angle from 21° to 20°. This point is 
marked C-3. These two modifications show remarkable 
agreement with the curve as drawn in Fig. 3. If these 
data were to be correlated on a thickness ratio basis or 
based on the total trailing-edge angle, the gain in 
critical Mach Number would not be consistent. 


(D) From Flight Data 

This airplane was prevented from going to higher 
Mach Numbers because of buffeting reported by the 
pilot. 


(E) From Flight Data 


This airplane encountered buffeting, buzz, and general 
instability at this Mach Number. It has been flown to 
a higher Mach Number, but more serious difficulties 
have been encountered. 








(F) From Wind-Tunnel Data 


This point was determined from tuft studies made on 
a wing of 12 per cent thickness at high Mach Numbers, 
The tufts began to be disturbed at slightly above a 
Mach Number of 0.75, and the separation is substan- 
tially complete at a Mach Number of 0.82. 


(G) From Wind-Tunnel Data 


This point is from a 10 per cent wing in which tuft 
tests indicated a disturbance starting at approximately 
0.80 Mach Number and increasing until separation is 
substantially complete at point 0.86 Mach Number. 


(H) From Wind-Tunnel Data 


Tuft studies indicate separation starting above 
Mach Number 0.85 and separation being substantially 
complete at 0.91. 

A study of tuft action on wind-tunnel models tested 
at high speeds indicates that the Mach Number for 
transonic separation on three-dimensional bodies will 
be 15 per cent higher than that shown in Fig. 3. This 
factor is in agreement with *the theoretical critical 
Mach Number difference between two- and _ three- 
dimensional ellipsoids. 

A comparison of the values of the separation Mach 
Numbers given in Fig. 3 with general data on the 
critical Mach Numbers of airfoils shows that transonic 
separation in general occurs at approximately 0.10M 
above the wing section critical Mach Number. This 
rough rule may be applied to predict separation Mach 
Numbers for wing-body intersections and tail junctures. 
Under certain conditions, however, the results may be 
slightly conservative, particularly for airfoils of the 
more conventional type where the position of maximum 
thickness is located well forward. For unfilleted junc- 
tures, using the assumption that the local velocities 
are additive, reasonable estimates of critical Mach 
Numbers can be obtained. When extensive filleting is 
used, low-speed wind-tunnel tests are necessary to 
determine the critical Mach Number and, hence, the 
approximate separation Mach Number. 


CONCLUSION 


It is concluded that for preliminary design purposes 
the afterbody angle of a lifting surface or body repre- 
sents a satisfactory index of the Mach Number at which 
airflow separation and buffeting will occur. 

For wing-body junctures and tail intersections it is 
suggested Mach Numbers be considered greater than 
the critical Mach Numbers by 0.10. 
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SUMMARY 


The variational principle in hydrodynamics was first studied 
by Hargreaves,! who-has shown that, when the hydrodynamic 
equations are satisfied, the integrand in the variational principle 
is a linear function of the pressure. For steady, irrotational 
flows, Bateman?:? has shown that the integrand is the pressure 
only, which can be written as a function of the velocity potential 
¢. Adirect method to obtain approximate solution in the case of 
subsonic flow of a compressible fluid can therefore be formulated 
following the Rayleigh-Ritz method instead of solving the original 
nonlinear differential equations. In applying this method to the 
case of a compressible fluid passing an arbitrary airfoil, however, 
there are two difficulties: first, since the fluid region is infinite, 
Bateman’s integral is not directly applicable; second, it is not 
easy to assume the potential function satisfying the boundary 
conditions. The nonlinear problem of compressible flow passing 
an arbitrary airfoil has been formulated in this paper. It is 
proved that, by conformal transformation, the velocity potential 
can be set up to satisfy the boundary conditions in the general 
case, and, hence, the Rayleigh-Ritz method can be applied. 
It is also shown that, by using the Theory of Residues in the 
integration and a modified Crout!® method in solving the resulting 
simultaneous equations, the computational labor can be mate- 
rially reduced. Since the variational principle is valid for sub- 
sonic, sonic, and supersonic flows, the possibility of applying it 
to the study of mixed flow is indicated. A numerical example of 
the nonlinear problem of the compressible flow around a circular 
cylinder is carried out, and the results compare favorably with 
those found by the existing methods. 


INTRODUCTION{ 


r THE THEORETICAL STUDY of subsonic flow of a com- 
pressible fluid, there are at least three different 
analytical methods for the approximate solutions of the 
complicated nonlinear differential equations. The 
methods are: the development of the velocity potential 
¢ into power series in the Mach Number, the develop- 
ment of ¢ in powers of a thickness parameter, and the 
hodograph method. The method of developing ¢ into 
power series in the Mach Number is generally known 
as the Rayleigh-Janzen Method. Its modifications and 
extensions bear the name of Poggi, Imai and Aihara, 
ioe February 17, 1948. Revised and received September 

* The author is indebted to Prof. Th. von Karman, K. O. 
Friedrichs, and C. C. Lin for their many inspiring discussions. 
To Profs. R. von Mises and H. S. Tsien, the author is indebted for 
their valuable criticisms. 

t Associate Professor of Aeronautical Engineering, Daniel 
Guggenheim School of Aeronautics, College of Engineering. 


t For a more detailed discussion of the many methods men- 
tioned, the reader may refer to reference 5, 
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and others. By this method solutions of first approxi- 
mation have been carried out in the case of flow passing 
a circular cylinder, elliptic cylinder, Joukowsky airfoil, 
and a sphere. Second approximation solutions have 
been carried out in the case of a circular cylinder and a 
sphere; third approximation, in the case of a circular 
cylinder.'? Computation beyond the third approxi- 
mation is practically unmanageable. In the case of 
flow passing a given arbitrary airfoil, the method is so 
tedious that it is almost impossible to carry out the 
computation beyond the first approximation. Since 
the solution is expressed in terms of power series of 
Mach Number, obviously for high-speed airflow, the 
method does not give good approximate results. 

The method of developing ¢ in powers of a thickness 
parameter gives the so-called Prandtl-Glauert rule as 
its first approximation. An iteration process for com- 
puting higher order approximations has been developed 
by Ackeret. Kaplan has carried out the calculations 
up to the third approximation in the case of flow past a 
bump. The computation proves to be laborious. 

The hodograph method transforms the exact non- 
linear equations of motion into exact linear equations, 
thus greatly simplifying the differential equations. 
However, in doing so, the boundary conditions are 
made more difficult to fulfill, The present-day applica- 
tion of the method is essentially based on two simplify- 
ing assumptions. First, y, the ratio of specific heat at 
constant pressure to that at constant volume, is as- 
sumed to be —1. Second, in the application of the 
well-known rule of Karman-Tsien, a distortion of the 
shape of the airfoil is assumed to be only slight and is 
In explaining the assumption of y = —1, 
is usually 


disregarded. 
the isentropic. pressure-density relation 
plotted with pressure and volume as the coordinate 
axes. In this case, the relation appears to be a straight 
line with y = —1. However, if they are plotted with 
pressure and density as coordinate axes, the discrep- 
ancy is much more amplified—i.e., the isentropic curve 
and the curve with y = —1 will have the second deriva- 
tive in different signs (Fig. 4). On the other hand, the 
distortion of the airfoil shape can now be calculated by 
the recent methods developed by Lin’ and Gelbart."* 
It can be shown that the distortion increases with the 
Mach Number and is actually large at large Mach 
Numbers. If an arbitrary airfoil is given, it is almost 
impossible to compute, by the hodograph method, the 
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flow past the same airfoil in the compressible plane at 
various Mach Numbers. 

In view of these drawbacks of the existing methods, it 
does not seem to be useless to make a study of another 
method based on an entirely different approach— 
namely, the variational method. The variational 
method has the advantage of being a direct method— 
i.e., by means of which one may be able to calculate the 
compressible flow passing a given airfoil at different 
Mach Numbers. Computation can also be made to 
obtain a solution of higher accuracy without additional 
mathematical difficulty. In the variational method, 
instead of solving the nonlinear differential equations of 
motion, a variational integral is formulated and 
approximate solutions are calculated by the Rayleigh- 
Ritz procedure. Braun,* in his dissertation at Gét- 
tingen under Prandtl, used the variational principle to 
study a linearized problem, but unfortunately the 
method was erroneously carried out. The nonlinear 
problem of a compressible flow past a given arbitrary 
airfoil has been formulated in the present paper. A 
numerical example is carried out for the case of flow 
past a circular cylinder, and methods for simplifying 
the computational work are indicated. The results 
compare favorably with Taylor’s® experimental results 
and the third approximation by the Rayleigh-Janzen 
method computed by Imai.'* Since the variational 
integral is valid for subsonic, sonic, as well as supersonic 
flows so long as the flow is irrotational, it is believed 
that the method can be used to investigate the mixed 
subsonic and supersonic flows before the appearance of 
a compression shock. This case is now being studied, 
and the result will be reported in a subsequent paper. 


THE VARIATIONAL PRINCIPLE 


It was first shown by Hargreaves! that, when the 
hydrodynamical equations are satisfied, the integrand 
in the variational principle is a linear function of the 
pressure. A variational principle for an inviscid com- 
pressible fluid was derived by Bateman.’ For steady, 
irrotational flows, Bateman has shown that the inte- 
grand is just the pressure, and he has further proved 
that the “pressure energy” is a maximum and the vari- 
ational problem is regular when the velocity every- 
where is below the local velocity of sound. 

Bateman’s work, in a simplified form, can be proved 
as follows. By denoting the velocity vector by gq, the 
Cartesian coordinates by x, x2; 3, and the corresponding 
velocity components by 1, %, and 3, the continuity 
equation for a steady flow is 


div (pg) = 0 
or 
: 
Y, @/rx:) (oni) = 0 (1) 


where’ p is the density. The condition of irrational 
flow is 


rot 9g = 0 
or 
(Ov3/Ox2) — (Ove /Ox3) = 0| 
(Ov, /O0x3) - (Ov3/0x;) = 0> (2) 
(2ve/2a:) — @v1/2xs) = Of 


With Eq. (2), the Euler equations of motion can be 
written as 

(p/2)(0g?/Ox;) = —Op/dx;, (i = 1, 2, 3) (3) 
where g? = v,’ + vw” + v;” and p is the pressure. For 
ideal gases, the pressure-density relation is of the follow- 
ing form: ° 


p= A+ Bp’ (4) 


where A, B, and y are constants. Eqs. (1), (2), (3), and 
(4) describe the steady irrotational motions of a com- 
pressible fluid. Bateman’s problem is to prove that 
the first variation of the integral 


JS, p dr = max. ; (5) 


is equivalent to Eq. (1) when Eqs. (2), (3), and (4) are 
satisfied, where dr is the elementary volumes of the 
fluid and the integration is to be extended to the whole 
volume of the fluid. 


Eqs. (2) are:satisfied if we introduce the velocity po- 
tential ¢ so that v, = 0¢/0x;. Since dp/dp = a*, where 
a is the velocity of sound, Eq. (3) can be integrated to 
give 

a? + [(y — 1)/2] g? = ao? + [(y — 1)/2] qo? (6) 
where the subscript 0 refers to the quantities at in- 
finity. From Eq. (4), we have 

dp/dp.= Byp’~' (7) 
Again noting dp/dp = a’, we obtain, from Eq. (6), 
Byp’~' = ao + [(y — 1)/2] @? — @*) 8) 
Combining Eqs. (4) and (8) gives 
p=A+ BOVEY yVO-Y fae + [(y — 1)/2] X 
et + hirow 


Since the constants are irrelevant to the variational 
problem, the variational integral (5) can be written as 


S” {ao? + (Cy — 1)/2] (qo? — 9?) }/0- dr = max. 


; (9) 
To simplify the integral, write 
dm? = [2/(y — 1) Jac? + qo? 
and the variational integral becomes 
Si (Qn? — 92)/0-) dr = max. (10a) 
or 
BK, (qn? — g2)/O-9 dr = 0 (10b) 
Since 
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jollowing the usual tensor notation that a repetition 
of the subscripts indicates the summation, Eq. (10b) can 
be written into the following form: 


iS, [Qm” ad 


Carrying out the first variation and dropping out the 
constant, it is obtained that 


Sh. lan? — (06/2x;) (06 /Ox,) 1/7 -9OG/Ox,)5 X 
(Od/Ox;) dr =0 


(0¢/dx))(0¢/dx,)]7/-Y dr = 0 (10c) 


(11) 
Combining with Eq. (8) and neglecting the constant 
coefficient, Eq. (11) is then 

S, p(0¢/dx;)(0/dx,)(6¢)dr = 0 (12) 
Itis obvious that 


o °¢ do O 2 ( ae 
5 =p——(8 to — lo» — 
dx, (- . Ox; ‘ Ox; Ox; vedi Ox, , 


and, therefore, Eq. (12) can be written as 


> (ose 28) — 56 2 (o S* ) | 
6g — bo - — lr = 0 
SN 2 (se % ; hat weer 
or P 
"2 ra) (» 2) 
— dr — 6 lr = 0 
[2 (2) : £ "ie aad 


Applying Gauss’s theorem to the first integral, we have 


re) 28) j Op 
dr = — 5g — dS 4 
Sz (+3 * r 3” ay (14) 


where 0 /On is derivative in the inner normal direction 
and S is the surface that encloses the volume V, and the 
integration is taken over the surface S. With the rela- 
tion (14) established, Eq. (13) becomes 


Od f o ( 2e) 7 
sae a mead — — d => ‘ 
fovs je dS F bo my ;, T 0 (15) 


In fluid dynamics, the boundary conditions are of the 
following two types: (1) if the boundaries are stream 
surfaces, 0¢/On = 0 on the boundary; (2) if the bound- 
ary is not a stream surface, the value of ¢ is prescribed 
and é@ = 0. Asaresult, pid(0¢/On) = 0 on the bound- 
ary, and, if the domain is finite, the first integral of 
Eq. (15) vanishes. Since V is arbitrary, the second in- 
tegral then gives us the continuity equation 


(0/Ox;) [p(0d/ox;)] = 0 


and it thus completes the proof. It is to be noted that 
the variational principle is valid no matter whether the 
flow is subsonic, sonic, or supersonic; however, the 


(13) 


integral can be proved to be maximum only in the case 
of subsonic flow. 


COMPRESSIBLE FLUID 







APPROXIMATE SOLUTIONS BY RAYLEIGH-R1ITz METHOD 


A great advantage in transforming differential equa- 
tions into variational integral is that one can then use 
the Rayleigh-Ritz Method to obtain expressions closely 
approximating the solution of the boundary value prob- 
lem. The Rayleigh-Ritz Method consists, essentially, 
of choosing a sequence of functions that all satisfy the 
boundary conditions but with undetermined constants, 
such as a;, 0;,, ¢;...andz = 1, 2,...,-Nin the functions. 
Let I be the variational integral. Then, by substi- 
tuting these functions into J and setting the first vari- 
ation of J to be zero, we obtain a set of simultaneous 
algebraic equations as follows: 


ar " ol da ol 
da, ' «= ac, 
er 


Solving these simultaneous equations, the constants 
can thus be determined. 

If the sequence of functions chosen is complete and 
the integral has a minimum or maximum, it can be 
proved that the solution so obtained is a minimizing 
or maximizing sequence to the original variational 
problem. An enumerable sequence of functions con- 
tained in the class R is said to be complete if every 
function in R can be approximated by a finite linear 
combination of the functions with preassigned accuracy. 
Although the method was originally applied by Ritz 
to a linear problem, the proof is also valid for non- 
linear problems. 

The Rayleigh-Ritz Method has been widely used in 
the field of Theory of Elasticity because the boundary 
conditions in these cases are relatively simple and func- 
tions that satisfy these conditions are not too difficult 
to choose. In fluid dynamics, the situation is differ- 
ent; the boundary conditions—e.g., flow around an 
arbitrary airfoil—are difficult to prescribe in simple 
analytical form. This is perhaps one of the reasons why 
the variational method has not been used very much 
in fluid dynamics. It will be shown in this paper that 
such difficulty can be overcome by using the idea of 
conformal mapping. . 


THE VARIATIONAL INTEGRAL IN INFINITE DOMAIN 


In applying the variational integral (5) to the case 
where the fluid extends to infinity, two difficulties arise. 
First, the integral becomes infinity, and it may be 
argued that a maximum or minimum is meaningless in 
this case. Secondly, the integral (14) may not be zero, 
although the integrand may approach zero at infinity, 
since S now is an infinite surface. Let us begin with the 
first difficulty. In the case of incompressible flow, it can 
be easily seen that Bateman’s variational principle is 
equivalent to Kelvin’s well-known theorem by mini- 
mizing the kinetic energy. In applying Kelvin’s 
theorem to the infinite domain, it is a usual practice® to 
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find the velocity potential of a corresponding flow that 
is at rest at infinity, and the actual velocity potential of 
the fluid is the sum of the velocity potential so obtained 
and that due to the uniform motion. In such a way, the 
variational integral becomes finite. This method, how- 
ever, cannot be applied to the compressible flow, since 
the principle of superposition is now no longer appli- 
cable. 

It can be noted that the addition or subtraction of an 
invariant function—i.e., a function whose variation is 
zero—is irrelevant to the variational principle. In other 
words, after carrying out the integration of the integral 
(5), all terms that do not contain the undermined con- 
stants may be subtracted from .the integral without 
affecting the variational principle. The resulting in- 
tegral may then be finite. In the subsonic flow, it is 
convenient, as will be shown later, to assume the ve- 
- locity potential as consisting of two parts, ¢: and ¢2, ¢; 
being the velocity potential of the corresponding in- 
compressible flow and ¢, the correction part due to 
the compressibility effect. The terms in p(¢) which do 
not contain the undetermined constants are equal to 
b(¢1), since ¢; is a definite function, 


5 S,p(¢1) dr = 0 
and the variational integral (5) can be replaced by 


Si[p(¢) — p(¢:)]dr 


* (16) 


max. 














1948 


With ¢; and ¢2 as assumed in the next section, the vag}. 
ational integral (16) is actually finite, provided that th 
integration is first carried out with respect to 6. This 
latter restriction is due to the fact that the integral j, 
improper, and one may arbitrarily choose the order 9 
integration in such cases. 

So far as the second difficulty is concerned, let ys 
again take the case with ¢; and ¢2 as assumed in the net 
section. Without dropping the constants, the fir, 


variation of Eq. (16) is equal to the expression (1j 
2 wa =) 
multiplied by ape it . We have then 
SNe 


+ = " =, 5) , oe 1 9 -1 
(Cle O)-C% 
— 1\/4a-0 L\qVG -» 
- (ai) fer oro] 
2Y 


- 
- l 
bb (6A;, cos nO + 6B,, sin n0) + 0 (4) 
n=1 rx 





a lie 


bp = 


In treating the problem where the fluid extends to 
infinity, we may first consider the fluid enclosed in a 
circular boundary and then extend the boundary to 
infinity. We therefore have: 


-U eos + 0(%) 


Se 
r2 


op 


on 


_o% 
or 


where 0(1/r?) represents terms of the order 1/r* ot 
higher. Substituting into the integral (14), taking 
y = 2, and integrating, we have 


9 ay — 1) 6 = 2s a 
E 1) f pod 2° ae | 
y¥-1 @=0 On r= 


—4(qn? = U?) Ur6A 


(17 


Similar results can be obtained by taking y = 1.5and 1.33; 
the significances of taking y to be these values will be 
discussed later. It is obvious that the expression (17) 
is not zero unless A; isa known constant. But thereis 
no reason why, in the general case of fluid dynamics, one 
should expect Ay to be any given value. If integral 
(14) is not zero, then the first variation of Eq. (16) isnot 
equivalent to the original differential Eq. (1). It is 
clear that Bateman’s principle must be modified when 
it is applied to the fluid extending to infinite region. 
The necessary modification can be obtained by sub- 
tracting from the integral (16) a term that is equal to 
Eq. (17), if ¢: and ¢: are so assumed. Then the final 
corrected variational integral in this case is as follows: 


I = Sy [b(o) — plod ]dS + 4(Qn? — U2)UrAn = 


maximum 


(18 


Note that the first variation of Eq. (18) is now exactly 
equivalent to the differential Eq. (1). If ¢: and ¢, ate 
assumed in other forms, similar terms can be worked 
out and should be subtracted from Eq. (16). 
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Braun, in his analysis of the linearized problem of 
compressible flow passing a circular cylinder, used 
Bateman’s integral without such corrections. His 
results are therefore erroneous. However, because of 
the mistake that he made subsequently in his calcula- 
tions, his final results turned out to be reasonable. 
present writer has checked Braun’s work and found 
that the coefficient A; is a large value with an opposite 
sign as given by Braun if using Bateman’s original 
integral. The results do not check at all with those 
found by Taylor and Imai. 


COMPRESSIBLE FLOW PASSING A CIRCULAR CYLINDER 


The boundary conditions for the two-dimensional 
potential flow about a body of arbitrary shape, in terms 
of polar coordinates, are as follows: (1) Atr = ©, 
g= Uord0¢/dr = U cos @ and (1/r)(0¢/00) = —U sin 
#; (2) on the surface of the body, 0¢/0n = 0, where n 
jsthe normal. To satisfy the boundary conditions, it 
is convenient to consider ¢ as consisting of two parts, 
¢, and de, where ¢; satisfies the above-mentioned con- 
ditions (1) and (2) and @» satisfies the condition (2) 
but at r = ©, O¢o/Or = (1/r) (O¢2/00) = 0. It is 
obvious that @ = ¢; + ¢@z2 satisfies the given boundary 
conditions. 

Let us consider first the case of a compressible flow 
passing a- circular cylindrical body. We can assume 
¢, and ¢» as follows: 


¢,; = Ul[r + (6?/r)] cos 6 — (x/2r)0 (19) 


M wN 
= DL sy [A mnRm(r) cos nO + 


m=1 n= 
BinnRm (r) sin n8) (20) 


where 
Rn(r) = (6+ 3 /mr™) — [b"+3/(m + 2)r" 77] 


in which 0 is the radius of the circular cylinder and x is 
the strength of circulation, An, and Bm, are undeter- 
mined constants. The boundary conditions are seen 
to be satisfied. 


Two-DIMENSIONAL COMPRESSIBLE FLOW PASSING A 
Bopy OF ARBITRARY SHAPE* 


For cylindrical bodies of arbitrary shapes, the circle 
can be mapped into any desirable shape by conformal 
mapping. If the mapping is conformal, even though 
the function ¢ is not harmonic, it can be proved that 
the boundary conditions are preserved—i.e., 


0¢/0r = U cos 6 
and 

(1/r) (0¢/00) = Usin@é, atr=o 
and 0¢/On = 0 on the surface of the body. 


* The author is indebted to Professor Lin who helped in form- 
wlating this section. 
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The proof can be established as follows: Referring 
to Fig. 1, let x = x(§, n) and y = y(é, »). Then 
(21) 


(2%) « 2) 
on = \dx/ \On Ov/ \On 
= V2) +R) 
On’ 0£/\Oon’ On/ \On’ 


where is the normal in the z-plane and mn’ is the nor- 

mal in the ¢-plane. Let s be the direction of tangent 

in z-plane and o be the direction of tangent in {-plane. 
From Fig. 2, it is seen that 


dx = ds cos a 





dx/ds = cos a 
also 
dy = dncosa 
dy/dn = cos a 
Therefore, 
dx/ds = dy/dn (23a) 
Similarly, 
dy ’ dx : dy dx 3b) 
— = sin a, — = —sin a, —=/-—-— (2 
ds on dn oo ds. dn , 
¥=|,.405 
/ 
/ ¥=2 
‘ 
/ 
é 
/ Y s-| 
YY 
¢ Po: fo) 
. c/ 
7 
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In the {-plane, we can establish the relations 
dt/da = dn/dn’, dn/do = —dé/dn' 


Substituting the relations (23a), (23b), and (24) 
into Eqs. (21) and (22), we obtain 


(24) 


Op. _ 0¢\/Oy o¢ Ox 

on - (22) + (*2) (=) (25) 
op - _ (2¢\(2n\ | (24\/0¢ aid 
ie (32)(22) ig (52)(5%) 638) 


Since 


2-@G)-@@) 
ox o¢/ \Ox On/ \dx ips 
o¢ 


dy - (535) 7 (2)(22) (28) 


by combining Eqs. (27) and (28) with Eq. (25), we 


ls 


have, P 
OF _ De (_ by DE , Br 2) 
on dé Os Ox r Os Oy t 
Del ele BR 
On Os oe Os 


If the mapping is conformal, (x, y) and (é, 7) must 
satisfy the following Cauchy-Riemann conditions. 


at/dx = On/dy,  dn/Ox = — dE/Dy 
Substituting Eqs. (30) into Eq. (29) and noting that 


ar ~ (Ss)(ae) + (SE)(3) 
a (SNS) + GS) 


it is obtained 
oe 7 ~ (ae) ae) + (S5)(3) > 
on a¢/\ds/ + \dn/\as) = 

(SE) (Se) (Gi) + SE)(5e)(G) = (Se) (G3) 

de ds) + \On/\do/\as) ~ \an'/\as 

The differential quotient do/ds in the above relation is 
in the form of total differentials because, for any 
mapping, ¢ is a function of s only. 

Denote the mapping function by ¢ = f(z). Then we 
have do/ds = f'(z). For conformal mappings, f’(z) # 
0. Therefore, if 0¢/0n in the z-plane is zero along the 
circle, 0¢/0n’ in the {-plane must also be zero along the 
arbitrary shape that is mapped from z-plane by the 
conformal mapping. - 

We have shown that, if the assumed potential ¢ for a 
compressible flow around a circle in the z-plane can 
be mapped conformally onto the ¢-plane, even though 
¢ is not harmonic, the boundary conditions are still 


preserved. Thus, if the mapping function that trans- 
forms a circle into an airfoil in an incompressible flow 


(30) 


(31) 


(32) 


is known, the expression for ¢ satisfying the boundary 
conditions in compressible flow around the same air. 
foil can be easily assumed. 

With ¢ thus assumed, we can substitute the expres. 
sion into the variational integral, differentiate the jn. 
“tegral with respect to the unknown constants, and set 
the final expression equal to zero. Thus, we have as 
many equations as unknown constants, and _ these 
constants can be determined by solving the equations 
simultaneously. 

In the actual computation, however, instead of sub- 
stituting the mapping function into ¢ directly, some 
simplification can be worked out. Let the velocity and 
potential function in z-plane and {-plane be denoted by 
q(z), o(z), and g(f), $(f), respectively. Then, accord- 
ing to the Theory of Complex Variables, 


stn ett) Ea bY) {free 
q(f)? -| de + On = 2 ae 


Pea +N -n 


If we denote the polar coordinates in the 2-plane asr 
and @ and those in the ¢-plane as R and ®, we have 


SS (aq 


{-plane 





mn? — g(f)*]”°"” RdRd® = 


ae ye” la 
ne dr d@ (33) 
vie Ee /dz|? fl eat ” 


s-plane 
1 = [HOT 4 [2] 
ate)’ = | Ox | r oy 


APPLICATION TO THE MIXED FLOW PROBLEMS 


where 





In establishing Ritz’s proof, the integral must have 
either a maximum or minimum, and it has long been 
regarded that the Rayleigh-Ritz Method cannot be 
applied if the boundary value problem does not lead 
to an equivalent extremum problem. However, Cour- 
ant and Friedrichs’ have recently pointed out that the 
Rayleigh-Ritz Method can still be applied in such a 
case. The mathematical significance in this case lies 
in the fact that we are seeking to satisfy the differential 
equation, as well as the boundary condition, not com- 
pletely but in the mean. 

Therefore, even though the variational integral (5) 
does not have an extreme value when the flow is not 
entirely subsonic, as long as the flow is irrotational 
or before the occurrence of a compression shock, a 
good approximation can be obtained by means of the 
variational method if there exists a solution in the mixed 


‘flow region at all. 


NUMERICAL EXAMPLE 


To illustrate the method, a numerical example has 
been carried out for the case of a compressible flow 
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passing a circular cylinder with unit radius without 
arculation. Referring to Fig. 3, since there is no circu- 
lation, the flow must be symmetrical with respect to @ 
or the horizontal axis, and, therefore, B,,, = 0. The 
horizontal velocity components must also be equal at 
points symmetrically situated with respect to the ver- 
tical, as well as the horizontal, axes such as points A, 
B,C. Asaresult A,,,, must be zero when m, n are even 
integers. 

The potential function @ can now be written as fol- 


lows: 


3 1 l 1 
o=U iar rt ae - — —]cos6 + 
r 3r? 
1 1 1 
As (? = 4) cos 36 + ra = +) cos 6 + 


1 1 ] 
ta( 3 _ +) cos'36 + Ais [ - - 53) cos 50 + 
" r 


3r? br 3r 
1 
51 a” 
or 


The variational integral in polar coordinates for two- 
dimensional flow is as follows: 


eda AP ihe ; oe" (: 2 ae 
Jr=1 iz E (2 r 00 * 


r dr d6+ 4(gm? al U)UrAyt = 0 (35) 





+) cos @ (34) 
r 


The integral can be integrated analytically when y/ 
(y — 1) isaninteger. For different values of y/(y — 1), 
the corresponding values of y are as follows: 


v/(¥ =~ ]) J 
1/, —] 
2 2 
3 1.5 
4 1.333 


The value of y for air is 1.405, which is between the 
values of 1.5 and 1.333. If the calculation is carried 
out at y = 1.333, a good approximate solution can be 
obtained. In the hodograph method used by Chaplygin 
and von K4rman-Tsien, the value of y is taken to be 
equalto —1. The variational integral in this case then 
becomes the familiar problem of determining the mini- 
mum surface in the Calculus of Variations. As pointed 
out by von Karman and Tsien, the velocity in the flow 


field would not differ too much from the undisturbed - 


velocity in the infinity—.e., only a small portion of the 
Pvs. p curve will be used. If one makes the curve with 
y = —1 tangent to the curve with y = 1.405 at the 
point that represents the conditions at infinity and uses 
the curve with y = —1, the error probably would not 
be too great. In our case, we may follow the same 
method and use y = 2. The advantage of using y = 
2 instead of y = —1 can be easily seen in Fig. 4. 
Although the curve with y = —1 and y = 1.405 may 
be made tangent to each other at (po, ao), their second 
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derivatives have opposite signs. Whereas by making 
the curves with y = 2and y = 1.405 tangent to each 
other at (0, po), there are small portions of the curves 
more or less close together, and therefore the curve y = 
2 provides a much better approximation. By using y = 
3 or y = 4, astill better approximation can be obtained; 
however, the amount of labor in carrying out the in- 
tegration and solving the simultaneous equations would 
also be materially increased. 

Substituting @ as given by Eq. (34) into Eq. (35) 
taking four constants As, As, Ags, Ass, and differen- 
tiating the resulting integral with respect to the con- 
stants, we have 


ol/0Au = 0, O1/0A13 = 0, 
OL/0As33 = 0, O1/0A15 = 0, 


WRda = ol (36) 
dT/dAs = Of 


II 


Integrating these equations, we obtain, finally,* 
(2.07407¢n" — 7.31852U*)An + 7.06667 U7Aj3 + 
(0.325926¢»,2 — 1.03704U?)As: + 0.915555 U?A33 + 
(0.1442869n? — 0.340136U")An = Ay 


7.06667 U?Au + (13.3333qm? — 44.2667 U*)Ais + 
1.20000U?A 3 + (1.600009,,? — 5.97334U*)Ags + 
32.4444U°Ass, + 0.440816U°As = K, 


(0.325926gm? — 1.03704U")An + 1.20000U?A13 + 
(0.09007419,7 — 0.207407 U?)A si + 0.192254U?A 33 + 
(0.0411428¢,? — 0.0810885U?)As = Ks 


0.915555 U*An + (1.60000g,,? — 5.97334U")Ai3 + 
0.192254 U?As + (0.277333qm? — 1.01181 U?)A33 + 
4.92444U?Ai5 + 0.0813062U?A 5) = Kg 


32.4444 U7A,3 + 4.92444U7A33 + (35.4444¢,,? — 
120.948U?)Ays = Ks 


(0.1142869,,2 — 0.340136U2)An + 0.440816U2A13 + 
(0.0411428q,,? — 0.0810885U?)As + 0.0813062U*A3s 
+ (0.02183099,.? — 0.0360739U2)Am = Ke 


where Ki, Keo, K3, K4, Ks, Kg are sums of nonlinear terms 
in Ay, Ai3, Asi, Ags, Ais, and As; and their complete ex- 
pressions are given in Appendix I. Denote A,’ = 
A,/aand K,’ = K,/a,’, for M = 0.4, we have 


3.309044 n’ + 1.13067Ais’ + 0.538074Ax’ + | 
0.146489A 33’ + (0)A15’ + 0.192435A sy’ = Ky’ | 


1.13067An’ + 21.7173A13’ + 0.192000An’ + 
2.50027A33’ + 5.19111A15’ + 0.0705306A 51’ 
= K,’ 


0.538074An’ + 0.192000Ai3’ + 0.161375A3)’ 
+ 0.0307606As3’ + (0)As’ + 0.0758944Am’ | (87) 
= K,’ 

0.146489A n’ + 2.50027A13’ + 0.0307606A 5)’ 
+ 0.437150As3’ + 0.787911A1s’ + | 
0.0130090A 51’ - Ky’ 


(0)Au’ + 5.19111Ais’ + (0)An’ + | 
0.787911A33’ + 58.0883A15’ + (0)An’ = Ks’ | 





* The numerical example was carried out by Messrs. Harry H. 
Hilton and G. V. R. Rao. 
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TABLE 1 
Velocity Distribution at @ = 2/2 from Various Methods 
* 4q Taylor’s 
a 2 ——————-————--———_——Variational Method a Rayleigh-Janzen Experimental 
ie One Two Three Four Five Six Method—3rd Approx. Results, 
> term terms terms terms terms terms y¥=2 y=1.4 vy = 1.407 
1.0 2.138 2.286 2.284 2.298 2.333 2.334 2.303 2.284 ‘<2 
1.04 2.057 2.199 2.198 2.210 2.244 2.244 2.215 2.197 2.13 
: Bp 1.805 1.874 1.924 1.931 1.959 1.959 1.934 1.920 1.89 
1.5 1.523 1.607 1.609 1.610 1.629 1.629 1.611 1.603 1.60 
2.0 1.298 1.348 1.351 1.349 1.360 1.360 1.347 1.345 1.35 
3.0 1.133 1.157 1.159 1.157 1.162 1.162 1.156 1.155 1.18 
0.192435A 1’ + 0.0705306A13’ + 0.0758944A 3’ (37 to four significant figures after decimal point, have the 
+ 0.0130090A 33’ + (0)Ais’ + pear following values: 


0.0413830As)’ = Ke’ | 


With Eqs. (37) so arranged, it is seen that the Crout 
method!® for solving a system of linear simultaneous 
equations can be used, and the constants Ay’ may be 
written in terms of K,’, Ko’, K3’, K4’, Ks’ and K,’ as 


follows: 
Ax)’ = 1.15571.Ky’ a 0.0185510 Ko’ 7 
9.60031K3’ + 0.0273402K,’ + 
0.00128699K,’ + 12.2554K.’ 
Ay3’ = --0.0185512K,’ + 0.135726K.’ + 
0.0111499K3’ — 0.770094K,’ — 
0.00168374K;’ + 0.0765767K.’ 
Asx’ = —9.60031K,’ + 0.0111473K.,’ + 


125.199K3’ — 0.155007K,’ + 
0.00110634K;’ — 184.936K,’ 


A33’ =:0.0273412K,’ — 0.770094K.2’ — 
0.155021K3’ + 6.75435Ky’ — 
0.0227960K5’ — 0.653607 K¢’ 


Aj’ = 0.00128700K,’ — 0.00168373K’ te 
0.00110620K3’ — 0.0227960K4’ + 
0.0176748K;’ + 0.00202234K,’ 
Agi’ = 12.2554Ky’ + 0.0765810K2’ — 
184.936K3’ — 0.653631K,’ + | 
0.00202211Ks5’ + 306.415Ke¢’ | 


+ (38) 





where it may be mentioned again that K,’, K,’, K;’, 
K,', Ks’, Ke’ contain nonlinear terms in Ay)’, A13’, Aa’, 
A33', Ais’, and Ag’. To solve for these constants, a 
method of successive approximations may be used. 
First, assume the value of A’s to be zero in the K’s; 
then A };’, Ais’, Agi’, As3’, Ais’ and As,’ can be obtained 
from Eq. (38). This‘is actually the solution for the 


linearized problem. Use the values of A’s so deter-: 


mined to compute K’s, and solve for these constants 
again. This then gives us a second approximation. 
Repeat the cycle until the desired accuracy is obtained. 
To make the convergence more rapid, after the third 
approximation, instead of following the method of suc- 
cessive approximation rigidly, one may extrapolate the 
results and use the values of A’s so obtained to repeat 
the process. 

For M = 0.4, Eqs. (38) are solved by the method of 
successive approximation and the constants, accurate 


An’ = 0.1038, Aj3’ = — (0.02264 An’ = — 0.06156 


A33' = — 0.03727, Aj,’ = 0.003156, Asi’ = 0.03036 


and the potential function is, therefore, equal to 


1 1 
$ = 0.4 (- cL *) cos 8 + 0.1038 (; a *) cos #@ — 
r r 


ao 
l Pa ei 1 
0.02264 |r + — ] cos 36 —0.06157 | — — =—]cosé@- 
r 3r* 5r5 
oe 1 4 
0.03727 |—- -— = cos 39 — 0.003156 x 
3r® 575 


l 1 1 
(, + -)cos50+ 0.03036( -_— +) cos § 
r 5r° ir? 


To compare with Taylor’s® experimental results and the 
third approximation of the Rayleigh-Janzen Method by 
Imai,!? the velocities at 9 = 7/2 are computed and are 
tabulated in Table 1. The velocity at 6 = 7/2 can be 
computed from the following equation: 


'9_ _10 )/ - ( 5) 
U a wert at 


3.8986 ( : 4) 4.8998 (2, * +) 
; 3r4 sr k 5ré Ss oT 


To get an idea of the convergence of the results, the 
velocities at these poinfs are also computed by taking 
the series with only one term—i.e., term containing 
Ay—two terms—i.e., An, Ais—and so forth until six 
termsare taken. It isseen from the results that the con- 
vergence of the series is satisfactory., It may perhaps 
be mentioned that if only An is used, the first equation 
in the system of Eqs. (37), which is actually 0//0Aun = 
0, must be used in solving for An. If two constants 
are to be used, say An and A;;, then the equations 
d1/0An = 0 and dJ/0A;: = O must be used in solving 
for these constants. The values of An are different 
when they are obtained from different calculations. 

The values of A’s so computed are as follows: 

(1) With term containing An only, we have Au’ = 
0.8307. Até@ = 2/2, the velocity is equal to 


: 1 2077 (4 - 3.) 
U = (1 a *) + 0.2077 2 374 





(2) Wi 
Ay’ = 0 
f= T/ 2 i 


(5) Wi 
and As, 
Ay’ - 


.003315¢ 
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THEORY 





(2) With terms containing Ay and Aj;, we obtain 
Ay’ = 0.09430 and Aj3’ = —0.02580. The velocity at 
j= 7/2 is equal to 


q ‘) i: (? ‘) 
S = ~ 1. —-_— — 
U (1 © ape LG ~ 


(3) With terms containing Ay, Ais, and Aj, we have 
Ay’ = 0.1026, Ais’ = —0.02581, and Amn’ = —0.04693. 
The velocity at 6 = 7/2 is equal to 


(iA) tose 3) 
2=(1+ —] + 0.4502 |— — —] - 
U fe r? 3r4 


0.1173 (2 +) 
ae 3r4 5r6 


(4) With terms containing only Ayn, Ais, As, and 
Ay, we obtain Ay’ = 0.1018, Ais’ = —0.02225, An’ = 
-0.04381, and A33’ = —0.03156. The velocity at 06 = 
7/2 then equals 


] 1 l 1 ; 

t. (14+ —)+ 04214(-- si) + 

l i Yo 3r4 
0.1271 (. : +) 
= 3r4 5r6 


(5) With terms containing only Au, Ais, As, Ass, 
and Ay;, we have Ay,’ = 0.1025, Ai3’ = —0.02265, 
Ay’ = —0.04297, Ags’ = —0.03722, and Ai,’ = 
(0033155. ‘The velocity at 6 = 2/2 is equal to 


q 1 a l 1 
~=(1+ *) + 0.4654 (- _ =) + 
U yr? r2 3r4 
Dae oa 
0.1717 374 + zy8 


DISCUSSION OF THE ACCURACY OF THE RESULTS 


If one can take infinitely many terms of a complete 
set to represent ¢, it can be.proved that Ritz’s procedure 
will give an exact solution. When finite number of 
terms are taken, the result is only approximate. Com- 
paring our results with those of Imai (also by taking 
y = 2), the variational method yields higher values of 
the velocity. However, it is to be seen (Table 2) that, 
in the Rayleigh-Janzen method, the values of the veloc- 
ities increase with the increase of the number of 
approximations. This suggests that probably the re- 
sults given by the third approximation in the Ray- 
kigh-Janzen method are still too small in comparison 
with the exact values. In this respect, the higher 
values given by the present method is toward the cor- 
rect side. 

On the other hand, the magnitude of the velocities 
given by the present method also increases with the 
lumber of terms taken. This again indicates that our 
results are probably smaller than the exact values. 
The effect of using y = 2 instead of y = 1.4 is, however, 
to make the velocities larger. This can be seen from 
the results of Rayleigh-Janzen in using y = 2 and y = 
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TABLE 2 
Velocity Distribution at 6 = 2/2 from Various Approximations 
of Rayleigh-Janzen Method (M = 0.4, y = 2.0) 
‘om 

\U . 

o% Zero First Second Third 
b ™\_ Approx. Approx. Approx. _— Approx. 
1.0 2.000 2.187 2.262 2.303 
1.04 1.925 2.104 2.176 2.215 
1.2 1.694 1.845 1.903 1.934 
1.5 1.444 1.552 1.592 1.611 
2.0 1.250 1.315 1.337 1.347 
3.0 1.111 1.142 1.152 1.156 





1.4 as tabulated in Table 1. It seems to suggest that 
these two effects may compensate each other to some 
extent on the final results. 


NOTES ON COMPUTATIONS 


In carrying out the integration of Eq. (36), much 
labor required in the computation can be reduced by 
using the following simplifications: First, for the inte- 
gration with respect to 6, the integral is of the type 


D F(cos 0, sin 0)d6 


where the integrand is a rational function of cos @ and 
sin 6 finite on the range of integration. The integral 
can be evaluated by writing e” = z, where z is a complex 
variable. Since 


1 l , 

= ° ol 4 _ eu_ @~! 

cosO= 5 (+2 ), sin @ = 5. tz g—*) 
the integral takes the form JeS(2)dz, where S(z) is a 
rational function of z finite on the path of integration 
C, the circle of radius unity whose center is the origin. 
By the Theorem of Residues,'' the integral is equal to 
2ri times the sum of the residues of S(z) at those of 
its pole which are inside that circle. Because the inte- 
gration of Eq. (36) involves many trigometric func- 
tions, much time can be saved and also confusion can 

be avoided by using this theorem. 


In the integration with respect to 7, the same integral 
may occur many times in different equations. It would 
be convenient to list all the different integrals on an- 
other sheet, and each integral is evaluated only once. 
Since the limits of integration are from r = 1 tor = © 
and since the integrand consists of terms of the type 
Si® (dr/r"), n > 1, which is simply 1/(m — 1) after 
integration, the result can be written down without 
thinking. 
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APPENDIX I 


Expressions of K,, K:, K;, Ki, Ks, Ke 


0 et ed hd a ha 


0.00698418UA5:? — 3.7926UA Ais + 0.512169U AAs — 0.633228U'A A353 + 0.186243 U AAs — 
0.580741 UAyAs + 3.14921 UA13A33 — 18.4593 UA13A15 — 0.195374U'A13A51 — 0.119789UAg1A33 + 
0.0331367 UAsAs — 2.79450UA33A15 — 0.0467426UA33:A 51 — 0.412699A117Ai3 + 0.118731A1;7Ag + 
0.0793648A 1:7A33 + 0.0466519A 17A51 + 3.24233A13?A11 + 0.602752A13"Agi + 7.72275A13"Ai5 + 
0.235043A 13"A 51 + 0.0276543A 912A 11 + 0.000211740A 217A 13 — 0.000251150As1:7A33 + 0.00380477A 517A5, + 
0.0993296A 332A 11. + 0.0206111A332A 31 + 0.225960A 332A 1s + 0.00869874A 332A 51 + 8.91429A152An + 
1.64995A15"Asi + 0.642938A)5"As + 0.00574819A 517A + 0.00276237A 51:7A13 + 0.00197180A51:7Ag + 
0.000545640A 5:7A33 + 0.203175An* + 0.00261873A* + 0.000361590A 5:3 — 0.0948150An AAs + 
1.07768A 1A 13A33 — 4.52064A 11A13A15 — 0.0228974A 1A 13A51 — 0.0218696A 11431433 + 0.0242076A Andy 
— 0.680635A 1:A33A15 — 0.00669122A 1:A33A51 + 0.210793A13A31A35 — 0.616721A13A3:A15 + 
0.00548246A13A9As + 2.63154A13A33A15 + 0.0853833A13A 33451 — 0.184894A1345451 — 
0.0823982A 3:4 33A15 + 0.000880140A 3:4 33451 — 0.0190871A33A15A 51 


Ky = —7.06667U* — 1.89630UA1,? — 0.0334396 UA3,;? + 0.00163270U'A;,? + 19.7333 UA Ai; — 0.580741 UA As 
+ 3.14921UAnA33 — 18.4593UAnA1s — 0.1953874UAnAs + 2.78857 UAi3As1 + 50.8444UA13A15 + 
0.912108UA43As1 ob 0.481269 UA3A33 _ 3.02138UA3:A15 =— 0.0132427 UA3:A51 ++ 7.63937 UA33A 1s = 
0.163629UA33As, — 1.07719UAyAm + 3.24233A17A13 — 0.0474075A1:"As1 + 0.538838A1:7A33 — 
2.26032A 11?Ais — 0.0114487A1:2As1 + 6.82287A13°A33 + 0.000211740A5;°An + 0.143916A5:7Ai3 + 
0.0261427A317A33 + 0.0203168A3:7A15 + 0.00178155Ag:7As + 1.20889A33"A13 + 75.3439A157A13 + 
12.5681A45"A33 + 0.00276237A5:?Au + 0.0343677A5:7Ai3 + 0.0115756A517Asi + 0.00657200A5:7A33 + 
0.00883425A 5:27A15 — 0.137566A un? + 13.6381A13? + 0.000838230A;3:? + 0.0744644A3;? + 
0.000233450A 5° + 1.20550A1A13As + 15.445541A13A15 + 0.470087A1A13A8 + 0.210793AnAsAs - 
0.616721AyAgAis + 0.00548246AyA9A51 + 2.63154AnA33A415 + 0.0853833A1A3A51 + 
0.184894A 1;A15A 51 + 2.69545A 134314 15 + 0.132909A 134 31A 51 +t 1.008874 134 15A 51 + 0.526137A3A33As » 
0.0248924A 3A 33A 51 + 0.00699350A 31:A15A51 + 0.218654A33A15A 51 


Ks = 0.631111U* + 0.256085UAn? + 1.39429UA,;? + 0.00495237UA3,;? + 0.0424126UA33" + 3.90349UA);* + 


0.000407200 UAs:? — 0.580741 UA1A13 + 0.0842327 UAyAs — 0.119789 UAA33 + 0.0331367 VA nAa - 
0.0668791 UA3A31 + 0.481269 UA43A 33 ev 3.02138 UA33A15 —_ 0.0132427 UA33As51 baad 0.0196402UA3As 7 
0.00326554UAs1As — 0.489791UA33A15 — 0.00644570 UA33As — 0.0474075A 1:7A13 + 0.0276543A1:2An - 
0.0109348A 1:7A33 + 0.0121038A117As) + 0.602752A13?An + 0.143916A13"Ag) + 1.34773A132Ar5 + 
0.0664544A 13"A51 + 0.00785619A3:2An + 0.00251469A 3:7A13 + 0.000646110A5:7A33 + 
0.00153984A3:7As: + 0.0206111A332A11 + 0.00556006A 332As1 + 0.0502731A33"A1s + 0.00276736A33°Aui + 
1.64950A 1:7Ai5 + 0.390548A15"A 3, + 0.180134A15’As1 + 0.00197180A51:7An1 + 0.00115756A517A13 + 
0.000882780A 517A +- 0.000329100A512A33 + 0.0395768A 1:3 + 0.000953470A3:° + 0.000177840A5:° + 
0.000423480A 1:4 134A 31 + 0.210793A 11A 13A 33 x 0.616721A1;A13A15 + 0.00548246A 1:4 13A 51 kame 
0.000502300A 1:4 314 33 + 0.00760954.A nAgAsi = 0.0823982A 1:4 93A15 + 0.000880140A 1:4 33A 51 oe 
0.0522855A 1A 31A3 _ 0.0406336A 134A 31A 45 + 0.00356310A13A3As1 os 0.526137A 313A 33A 15 + 0.0248924 x 
A 3A33A51 + 0.00699350A 13A 15451 — 0.00112654.431:A33A 15 + 0.000955200A 314 33A 51 + 0.00438570A s3A sds 
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— 0.800000 U* — 0.316614UA,;? — 0.00982009 UAs,;? — 0.000742200UA5,? + 3.14921UAnAis — 
0.119789UAyAs + 0.551112UAyA33 — 2.79450U An Ais — 0.0467426UAnAs + 0.481269UApAn + 
7.63967 UA33Ais + 0.163629UA13As: + 0.0848252UAsA33 — 0.439791 UAsi:Ais — 0.00644570UAg:A5: + 
1.19873 UA33A15 + 0.0287347UA33A5, — 0.147785UAAs + 0.538838A1:2Ai13 — 0.0109348A1:°An + 
0.0993296A 1:°A33 — 0.340318A17A15 — 0.00334561A112As + 1.20889A13"A33 — 0.000251150A3:?An + 
0.0261427A 3:°A1s + 0.00556006A 3:2A33 — 0.000563270A3:7Ais + 0.000477600A3:7A 51 + | 
0.223393A 33"A1s + 12.5681A18*A13 + 2.27753A152A33 + 0.000545640A 512A + 0.00657200A517Ai3 + 
0. 000329100A 5;2A 3) + 0. 0015 oz 2112A517A33 - 0. 00238452A 517A 1s = 0. 0264549A 1,3 — 3. 27429A;3 a 
0.000215370As:* + 0.0143717A33* + 0.0000720001A 5° + 0.210793A1A13Am + 2.631544 1A13Ai5 + 
0.0853833A 1A 13A51 + 0.041222141,A31A33 — 0.0823982A1A3A15 + 0.000880140A1A3Am + 
0.451920A11A33A15 + 0.0173975AnA33A8 — 0.0190871AnAisAn + 0.526137A13;AnA15 + 
0.0248924A 13;A31As: + 0.218654A 134 15d 51 + 0.100546A9A33A15 + 0.00553472A9AgAn + 
0.00438570A 3A1sAs1 + 0.0453130A33:4 An 


K; = 25.4222 UA1s? ate 0.599365 UA33? _ 18.4593 UAy,A = 2.79450 UA A333 + 34.8148 UAyA Tees 3.02138 UA A 31 


+ 7.63937 UA13A33 — 1.07719UA13As1 — 0.439791 UA1A33 + 7.80699UAg1A15 — 0.147785UA3An + 
2.56097 UAysAs — 2.26032A1:7Ais — 0.340318A 11:7A33 + 8.91429A 11:7Ais + 7.72275A13"An + 
1.34773A13"Aa + 75.3439A13"Ais + 0.504433A13?A51 — 0.0203168A5:7Ai3 — 0.000563270A3:2A33 + 
0.390548A 3:°Ais + 0.225960A 33"A 11 + 0.0502731A 337A 31 + 2.27753A55°A1s + 0.0026565A 33"A 5) + 
0.00883425A 5:7A1s3 + 0.00238452A 5:7A33 + 0.0930836A 51?A15 + 104.102A15* — 0.616721A1A13As1 + 
2.63154A 11A13A33 — 0.184894A11A13A51 — 0.0823982A11A31A33 + 3.29990A1A31A15 — 0.0190871A1A33A51 
+ 1.285864 1A15As1 + 0.526137A13A31A33 + 0.00699350A 13A 1A 51 + 25.1361A13A33A15 + 0.218654A13A 55- 
Asi + 0.00438570A3:A33A 51 + 0.360269A3:A 5A 51 


7UA3,? + 0.0143673 UA3;? +1.28048UA,;? — 
0.0000465000 UA 5;* — 0.195374UA1Ais + 0.0331367 VA Aa — 0.0467426 UVAnA33 + 0.0139684 UA As — 
0.0132427 UA13As; + 0.163629 UA13A33 — 1.07719UAisA1s5 + 0.00326540 UA13As, — 0.00644570 UA31A33 + 
-0.000814400UA31As — 0.147785UA33A15 — 0.00148440UA33As — 0.0114487Ay)7Ai3 + 
0.0121038A1,?As1 — 0.00334561A 1:?A33 + 0.00574819A17As1 + 0.235043A13"A11 + 0.0664544A13°As) + 
0.504433A 137A 15 -- 0.0343677A 137A 51 a 0.00380477A3:2An + 0.00178155A3;7Ai13 + 0.000477600A 31°A 33 + 
0.000882780A 3:2As: + 0.00869874A 332A 11 + 0.00276736A53"As1 + 0.0226565A53*A1s + 0.00152112A53°A 5; 
+ 0.642928A 15°Ai + 0.180134A15’Agi + 0.0930835A 157A si + 0.00108477A51*An + 

0.000700350A 5:7Ai3 + 0.000533520As:7As1 + 0.000216000A 51:7A33 + 0.0155506A 1)? + 0.000513280A5:* + 
0.000112860A 51° + 0.00548246A 11:A13A 1 + 0.0853833A 11:A13A33 — 0.184894411A13A15 + 
0.00552474A1A13sAs: + 0.000880140A1A3A33 + 0.00394360AnA3As — 0.0190871A1A3A15 + 

0. 001091284 11:4 33A 51 + 0. 02489244 13A 3:A 33 + 0. 006993: 50A 134 3:A 15 4 0.00231: 512A 13A3 1A51 ee 

0.218654A 13433415 + 0.0131440A 13433451 + 0.0176685A13A1 sAsr + 0.00438570A 3:A33A15 + 

0. 0006: 58200A 31A 33A 51 a 0. 00476904A 334 15A 51 
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The Theoretical Induced Deflection Angle in 
Cascades Having Wall Boundary Layers 


GEORGE F. HAUSMANN*. 
United Aircraft Corporation 


ABSTRACT 


A theoretical analysis is presented which employs a lifting-line 
theory assuming incompressible, nonviscous flow to determine the 
induced deflection angle in two-dimensional cascades having wall 
boundary layers. Two correction factors are derived for the de- 
flection angle of the idealized infinite cascade: one that accounts 
for the induced deflection due to the trailing vortex system in the 
wall boundary layer, and the other that accounts for the de- 
crease in total circulation of the bound vortex system. 

Plots are included which show the effect of blade length, bound- 
ary-layer thickness, and the number of blades in cascade on the 
deflection angle of several typical cascade configurations. 


SYMBOLS 
0, 0’ = arbitrary coplanar points equidistant upstream and 
downstream from the lifting line 
h = normal distance from 0 and 0’ to the lifting line 
r = distance from lifting-line element to 0 and 0’ 
ds = lifting-line length element 


0,§,7,¢ = angles associated with points 0 and 0’ and the lifting 
line (see Fig. 2) 


l = lifting-line length or cascade blade length, in. 

lx = circulation 

To = free-stream circulation of a finite cascade blade or 
lifting line 

y = distance from tunnel wall to trailing vortex origin, 
in. 

q = induced downwash of bound vortex system 

w = induced downwash of trailing vortex system 

Cr = cascade lift coefficient 

L = cascade lift, Ibs. 

T = blade spacing, in. 

b = blade chord, in. 

B = air angle, deg., referred to cascade line 

V = velocity, ft. per sec. 

p = density, slugs per cu.ft. 

Cr = axial velocity component, ft. per sec. 

Cu = velocity component parallel to cascade, ft. per sec 

65 = boundary-layer total thickness, in. 

6* = boundary-layer displacement thickness, in. 

K = circulation correction factor 

n = number of half-span lengths along, lifting line, or 
cascade blade, image system 

S = blade relation within tunnel referred to tunnel 
centerline 

FP = blade or image number referred to tunnel centerline 

B = number of lifting lines (oda), or blades, in cascade 

R = (P — S)/B 

Subscripts 
1 = conditions at cascade inlet 
2 = conditions at cascade outlet 


m = mean conditions where cot Bm = (cot 8; + cot B2)/2 


Received September 24, 1947. Revised and received June 28, 
1948. 
* Aeronautical Research Engineer, Research Department. 


Superscripts 
’ = finite cascade characteristics corrected for induced down- 
wash of trailing vortex system 
” = finite cascade characteristics corrected for induced down- 
wash of trailing vortex system and modified downwash 
of bound vortex system 


INTRODUCTION 


iw LACK OF CORRELATION of deflection angle data 
from various cascade tunnels is primarily attributed 
to variations in the ratio of blade length to wall bound- 
ary-layer thickness, assuming that Reynolds Number 
and other such effects are equivalent. As a means of 
comparing cascade data from different tunnels on a 
common basis, it has been proposed to correct these 
data to “‘infinite’’ conditions by consideration of the in- 
duced velocities of the modified vortex system that is 
inherent with cascades having wall boundary layers. 
In theory, such corrected data could then be adapted to 
axial-flow compressor design by applying the induced 
effects that are to be expected in the elements of a 
rotating blade row. 

Although a theoretical analysis of this nature does 
not contemplate the complicated secondary flows that 
are associated with retarding cascades, it is felt that the 
form of the variables in the derived equations will be 
helpful in compiling more accurate empirical relation- 
ships for axial-flow compressor design. 


THEORETICAL ANALYSIS 


The vortex system of an arbitrary two-dimensional 
cascade with side-wall boundary layers is shown in Fig. 
1. The cascade vanes are replaced by lifting lines that 
induce components upstream and downstream and 
parallel to the cascade plane. As a means of simplify- 
ing the analysis, the end walls are assumed to be curved 
to match the streamlines that would be associated with 
a cascade having a prescribed circulation and an infinite 
number of lifting lines. By this expedient, the images 
of the bound vortex system have no reversal of sign, 
and the downwash expression may be readily derived. 

The trailing vortex system shown in Figs. la and 1b 
is formed by the gradient of circulation in the low veloc- 
ity, side-wall boundary layer at the blade extremities. 
At any arbitrary downstream station where the exit 
streamlines are straight, these vortices have images if 
the two dimensions of a plane that is normal to the exit 
air-stream direction. In order to satisfy the boundary 
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conditions, the vortex images that extend in the cascade 
direction have a cyclic variation of sign that depends on 
the number of blades in the cascade configuration. The 
side-wall images of the trailing vortex system are posi- 
tive in sign as are the side-wall images of the bound 
vortex system. 


Bound Vortex Correction 


The presence of the side-wall boundary layers results 
in a decrease in circulation. at the blade extremities 
which is reflected in the infinite number of images along 
the lifting line. In a cascade this cyclic pattern of the 
circulation produces a downwash angle that differs 
somewhat from the downwash resulting from a cascade 
having infinite lifting lines with constant circulation. 
The magnitude of this correction is evaluated from the 
following analysis. 


Consider the line vortex diagram for an airfoil of 
length, /, in a tunnel bounded by end walls having 
boundary layers of thickness, 6, as shown in Fig. 2. 
The circulation around this line has a cyclic variation 
that repeats from + © to—©. OQ and 0’ are arbitrary 
points equidistant upstream and downstream from the 
cascade, hf is the normal distance from 0 to the line 
vortex system, and 7 is the distance from 0 to any ele- 
ment, ds, on the line. 


Now, referring to the vector diagram in Fig. 3, 
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point, 0, resulting from an element, ds, of a straight line 
vortex of length, AB, is 


Ids . 
dq = —— sin @ = 
4nr* 


x S odd 
4xh fe 


(1) 
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The total induced downwash of the line vortex is 
obtained by integrating from ¢ = —[(m/2) — &] to 
@ = +[(2/2) — y] so that 

gq =“(I'/4xh) (cos — + cos y) (2) 

Considering, now, each cascade lifting line and each 
of the wall images to be elements of an infinite line 
vortex, the induced downwash is 


qgq= JS (Tds/4ar?) sin @ 


or, with reference to Fig. 2, 
/ *31/ 
SB: atts 31/2 oe 
= ——— gift ——~ sin 9 + 
q —1/2 4qr? 1/2 Amr? 


—1/2 Ww. 
I Ids . in, 
sin a 

—31/2 4mrr? 


Assuming that / is large compared with //2, so that 
sin 0/4rr? may be considered constant over each ele- 


ment of integration, 
: } 31/2 
sin 6 | 
+ —| fi Ids 
Amr? | 5 31/2 S 1/2 


Noting that each integration extends over the limits 
of one blade length and substituting KT,/ for each 
integrand, 





(3) 


sin 6 | 


“ } 1/2 
- i Ids 
Amr? |5a1/2 J-i/2 


(4) 





= | 


sin 6 | sin 6 
ee | rata i ae | 





g = KT + 
q 4rr? ls =| /2 , 
where K = St I'ds/T)l and Ty is the value of the circu- 
lation outside of the boundary layer. Since the sum- 
mation of the sin 6/42r? terms extends to infinity in both 


directions from the measurement plane, then 


gq = J" (KT» sin 6/4nr?)ds 


4nrr | s =31/2 


(6) 
Comparing this expression with Eqs. (1) and (2), it is 
evident that Eq. (2) may now be written as 


g = (KT/4mh)(cos — + cos y) (7) 
Now, since £ = —aand y = +7, 
q = KIT)/2rh (8) 


From the law of Kutta-Joukowsky, the circulation may 
be expressed as 


L = pV,Il (9) 
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where JV, is the mean velocity through the cascade. 
The cascade lift is given by the familiar expression 


L = '/2pVn?Cybl = prCzVm(cot Bi — cot B2) (10) 


where 7 is the blade gap and the angles are self-explana. 
tory in Fig. 1. 

Now, solving for the circulation in terms of the cas. 
cade lift coefficient, 


yl = L/p Vol = VinCzb/2 (11) 


By assuming a turbulent boundary-layer velocity profile 
where V/V, = (y/6) "kK may be evaluated as follows: 
ja Tds — Vo(l — 28) + VnCrb Ar® (y/6)'"dy 


K=-22'2 
Tol I'd 





(12) 


where 6 is the boundary-layer thickness and V,,C,b/2 = 
To. Then, 

- _ Tod — 206 + (14/s)500 é 
= = 1 = SS 
lo 41 
From the downwash derivation of Kutta-Joukowsky, 
the total downwash, or vector difference of velocities, 
parallel to an infinite cascade of infinite length lifting 
lines is 





K (13) 


AC, = To/r = C,(cot Bi; — cot Be) (14) 


From the preceding analysis on the cyclic variation of 
circulation along the infinite lifting line, it is evident 
that the corrected downwash resulting from the bound 
vortex system of a two-dimensional cascade having wall 
boundary layers is 


AC, = To/Kr = C,(cot B; — cot B2)/[1 — (6/4/)] (14) 


Trailing Vortex Correction : 

The induced downwash parallel to the cascade which 
results from the trailing vortex system shown in Figs. la 
and 1b is derived on the basis of Eq. (2), where 


q = (T'/4mh)(cos — + cos y) 


Since the induced downwash due to trailing vorticity is 
considered at an infinite distance downstream from the 


cascade § = —zand y = wso that 
g = T'/2ah (16) 


Applying this general expression to the cascade of lifting lines having wall boundary layers such that the strength 
of trailing vorticity is (dI'/dy)dy, the induced downwash, w, resulting from the complete image system is 


9 


1 ° ar 


(1/2) — y 





w= — 


On expansion, this becomes 


8[(12/4) — y? — 7? sin? Be] 


(1/2) + y be. a] 


2 2 
2r Jo dy y (>, ye (1/2) + 9 ia — y]? + 7? sin? 6, [(2/2) + y]? + 7? sin? By! 


(17) 





:. se ( 4 \ 
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At the blade center-span where 0 < y? < 


am+o R=o g=(B—1)/2 


INDUCED DEFLECTION 


6° is small compared with /*/4, the y® terms may be omitted, and the 
downwash resulting from the complete trailing vortex system may be expressed as a sum in the following manner: 


2n + 1)%(/2/4) — (BR 


ANGLE IN CASCADES 6x9 








27? sin? Bo 





erty be 


ee, R=—o gs=-—(B-— 


where the m limits refer to the number of half span 
lengths in the lateral direction, the S limits refer to the 
relation of blades in the tunnel, and the R limits refer 
to the images in the longitudinal, or cascade, direction. 
If P is the blade, or image, number referred to the tunnel 
centerline, R is related by the expression 
R= (P — §$)/B 

where B is the number (odd) of lifting lines in the finite 
cascade. 

No simple solution for this triple summation has been 
found, nor has its convergence been definitely estab- 
lished. The form of this equation is, however, similar 


Therefore, 


n= o@ 


> 


w= (VC b6* /27) 


Substituting Eq. (10) for (V,,C,b/2), 


w 76*(cot Bi - — cot rt Be) > : > | 
Cy 


Application of Corrections 


The induced velocity vectors that are evaluated froin 
the equations derived in the preceding analysis are 
applied to the finite cascade velocity vector diagram 
to determine the equivalent infinite cascade perform- 
ance. 

Referring to the vector diagram in Fig. 4, Vi and §; 
are the inlet velocity and direction to the actual cascade 
and also to the corrected infinite cascade. V2 and p, 
refer to the velocity vector at the exit of the arbitrary 
cascade having wall boundary layers. V2’ and 2’ 
describe the exit velocity vector that has been modified 
by the application of the tangential component of the 
induced downwash, w, resulting from the trailing vortex 
V2.” and 8,” describe the velocity vector that 
also includes the change of the downwash, AC,, result- 
ing from the nonuniformity of circulation around the 
lifting lines. V2" and #6” thus define the outlet velocity 


system. 


of the arbitrary cascade having wall boundary layers’ 
but corrected to infinite conditions. 





1, [2m + 1)?(/?/4) +.(BR 


Se? (dV /dy)ydy = ( 


n=—-o R=—ow g=—(B-1) 


(19) 


ar 
—1)* bj — ydy 
0 ay 


S) 
S)? 


+/+ 


27? sin? Bo]? 


to the expression derived by Glauert' for the induced 
downwash of a single airfoil bounded by tunnel walls. 
Glauert’s equation has been empirically correlated and 
is generally accepted as authoritative. A numerical 
substitution of values in Eq. (19) indicates convergence, 
although this test cannot be considered decisive. 

The value of the integral in Eq. (19) may be estab- 
lished by again assuming a turbulent boundary-layer 


velocity distribution where V/V, = (y/6)'”. Then 
r =a | 9 VinCyb(y/8)' ly 
and 
w(dT'/dy) = (*/2) VinCrby'"6—' 
VinC,b/2)(} (20) 


and since the displacement thickness, 6*, of a seventh root boundary-layer velocity distribution is equal to '/s6,- 


=(B—1)/2 


: ae 


(21) 


° 


@ S=(B—1)/2 


(2n + 1)?(7* 
(Qn + 142 


}) bead (BR +s 2a * sin® * Ba 
t) + (BR + s)? by? sin? Bo]? 


© a te ae 


From the vector diagram and the preceding analysis, 


cot 8, = C,,/C,, cot Bs = Ci2/C- 


cot Be’ = Cy’/C,, cot B” = Cyw"/C; 
sin Bs = (Cw’ — Cy)/w 
Substituting 
cot Bo’ = (wsin Bs + Cyw)/C- , (23) 
and 
AC, = (Cur — Ca’)/K = Cu — Cr” 


Substituting for AC, from Eq. (15), 





cot 6; — cot Bs — (w/C,)sin Bs 
cot Be” = cot 6; — Br — : A [C, B (24) 
1 — (6/41) 
DISCUSSION 


The induced deflection angles of several typical cas- 
cade configurations have been calculated and are plotted 
in Figs.5and6. For the example, a cascade inlet angle, 
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Ai, of 30°, outlet angle, 82, of 90°, and a gap-chord ratio 
equal to 2 have been chosen. The circulation required 
to produce the 60° turning is independent of vane 
camber or cascade stagger angle. 

The form of the linear cascade variables in Eq. (22) 
excludes the use of nondimensional parameters so that 
the examples have been calculated in terms of absolute 
cascade dimensions. The summation of terms from 
Eq. (22) has been carried out in symmetrical order and 
includes all downwash terms that are greater than 1 per 
cent of the induced downwash resulting from the vane 
at the tunnel centerline. 

Referring to Fig. 5, it is seen that the number of vanes 
in the finite cascade has little effect on the induced de- 
flection angle of the arbitrary configuration. It should 
be realized, however, that the inlet walls to the cascade 
are assumed to be curved to match the inlet streamlines 
so that the upstream potential flow field is inherently 
independent of the number of vanes in the cascade row. 

The curves in Fig. 5 also show the effect of blade 
length on the induced deflection angle of a cascade that 
has a given blade spacing (2 in.), chord (2 in.), and 
boundary-layer total thickness (0.24 in.). It is inter- 
esting to note that an increase in blade length results in 
a reversal in the sign of the induced deflection angle so 
that a particular combination of blade length and 
number of blades actually results in a zero correction for 
the induced deflection angle of the specific finite cascade. 
An aspect ratio parameter has not been applied to the 
curves in Fig. 5, since the blade length and spacing in 
Eqs. (22) and (24) are independent variables for a 
given aspect ratio and gap-chord ratio. 

The effect of wall boundary-layer thickness on the in- 
duced deflection angle of an arbitrary cascade is shown 
in Fig. 6. The curve is not extended beyond a value of 
5* equal to 0.10, since the simplifying assumption used 
in obtaining Eq. (18) limits the analysis to small values 
of boundary-layer thickness. 

The principal limitation of the application of the 
theoretical correction to the induced deflection angle 
appears to be in the use of a realistic value of wall bound- 
ary-layer thickness and velocity distribution. Al- 
though the theory assumes a lifting line in place of the 
cascade vane, the boundary layer through an actual 
cascade changes shape and thickness along the vane 
chord, so that true representative boundary-layer condi- 
tions are difficult to establish. It is also probable that 
the complex, secondary flow vortices within the blade 
passages have considerable influence on the induced 
downwash. In spite of these shortcomings, it is felt 
that the theory presented in this paper could be useful in 
establishing a first approximation to the induced effects 
in cascades, 
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Surface-Pressure Gradient and Shock-Front 
Curvature at the Edge of a Plane Ogive 


with Attached Shock Front 


M. M. MUNK* ano R. C. PRIM* 
Naval Ordnance Laboratory 


SUMMARY 


The flow of air at high Mach Numbers past a sufficiently sharp 
plane ogive is characterized by the presence of a curved shock 
front attached to the edge of the ogive. This paper presents 
and carries out a simple method for computing the exact surface- 
pressure gradient and shock-front curvature at the edge of the 
ogive 

A comparison is made between the exact edge-pressure gradient 
and that given by two approximate methods of computing the 
pressure distribution along an ogive. 

The numerical results, presented graphically, cover all ogives 
with attached shock fronts for Mach Numbers up to 13. 


INTRODUCTION 


Sn FLOW OF A PERFECT GAS at high Mach Numbers 
past a sufficiently sharp plane ogive is characterized 
by the presence of a curved shock front attached to the 
edge of the ogive. Even when the vorticity and non- 
uniformity of entropy resulting from the curvature of 
the shock front are negligible, the detailed computation 
of the flow field behind the shock front requires a 
laborious step-by-step calculation for each special case. 
Hence, the application of simpler approximate methods 
to the computation of pressure distributions on bodies of 
this type is highly desirable. The degree of accuracy to 
be expected from the application of any of these 
methods to a specific case is, however, difficult to ascer- 
tain. 

This paper presents and carries out a simple method 
for computing the exact surface-pressure gradient at the 
edge of the ogive. This information is useful both in its 
own right and as a means of appraising the reliability of 
results obtained by approximate methods. 

The curvature of the attached shock near the edge is 
alsocomputed. This information is of assistance in the 
initial steps of detailed computations of flows of this 
type by characteristic methods and for interpretation of 
experimental results. 


THE ‘‘WEDGE + PRANDTL-MEYER” AND ‘“‘TANGENT- 
WEDGE”’ APPROXIMATIONS 


An obvious and reasonable approximation to the 
pressure distribution on plane-ogival bodies can be 
obtained from the well-known wedge and Prandtl- 
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Meyer solutions. The shock front is taken as that fora 
flow at the given Mach Number impinging on a wedge 
of the same included angle as the edge angle of the 
ogive. The flow along the curved wall is then com- 
puted as a simple Prandtl-Meyer flow, with initial con- 
ditions determined by the straight shock front estab- 
lished as above. 

A second plausible approximation procedure is to 
assign as the pressure at each point on the ogival surface 
the surface pressure in a flow at the given Mach Number 
against a wedge tangent to the ogive at the point in 
question. 


CORRECTION OF THE ‘‘WEDGE + PRANDTL-MEYER”’ 
APPROXIMATION AT THE OGIVE EDGE 


The approximation involved in the “Wedge + 
Prandtl-Meyer’’ method lies in the neglect of the inter- 
action between the shock front and the successive dis- 
turbance lines originating along the ogive surface. 
Actually, this interaction gives rise to reflected families 
of disturbance lines which progressively weaken the 
shock front so that its inclination decreases with dis- 
tance ‘along the shock, approaching asymptotically the 
Mach angle of the undisturbed flow. (At points re- 
mote from the body, the difference between the shock 
angle and the Mach angle decreases more rapidly than 
the inverse of the distance. Otherwise, the wave drag 
can be shown to be infinite. ) 

By limiting consideration to a region sufficiently near 
the ogive edge that the pressure, Mach angle, shock 
angle, and streamline inclination and the first deriva- 
tives of these quantities do not differ appreciably from 
their initial values, the initial rates of change of these 
quantities will be evaluated. This amounts to per- 
forming a linearization that approaches exactness as the 
régime considered shrinks to the edge itself. 

Fig. 1 represents the flow régime near the ogive edge. 
The ogive wall OBC is inclined at the angle 6 to the on- 
coming flow. The rate of change of the inclination 6 
with distance along the wall is determined by the given 
wall curvature. The shock front OAD is inclined at the 
angle a to the oncoming flow. Its curvature must be 
determined so as to be consistent with the shock rela- 
tions and with the strengths of the basic and reflected 
disturbances arising from the curvature of the wall. 
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BD is a member of the basic family of disturbance lines 
originating at the wall; AB and DC are members of 
the family of disturbance lines reflected from the shock 
front. The “strength” of these disturbance lines may 
be defined in terms of the rates of change of pressure 
along members of each family. That is, a measure of 
the intensity of the basic family of disturbance lines 
originating at the wall (e.g., BD) is the rate of change of 
pressure along a disturbance line of the reflected family 
—i.e., Op2/0C;. Similarly, the rate of change of pres- 
sure along a basic disturbance line (Op2/0C,) is a meas- 
ure of the intensity of the reflected family of disturb- 
ance lines. These rates of change of pressure along the 
two families of disturbance lines and the curvature of 
the shock front may be considered as the three basic un- 
known quantities. The conditions necessary for their 
evaluation are furnished by the two characteristic 
equations relating the changes in pressure and changes 
in flow direction along disturbance lines of the same 
family and by the shock relations, which relate the 
changes in pressure and changes in direction along the 
shock front. 

The procedure will be to impose the conditions that 
(1) the pressure difference between points A and D 
immediately behind the shock front must be the same 
when computed from the shock relations along AD or 
from the characteristic relations along ABD and (2) 
that the change in flow direction computed along BDC 
be the same as that given by the prescribed curvature 
along BC. 

Letting P = p2/p; and applying the first condition, 
we obtain 


dP oP aP 
-_ = , = 64 = = = _—- 
in iad (5 ) (6p — 64) i. (1) 


But, . 
bp — 54 = (06/0Ci)a + (06/0C2)c (2) 


and from the characteristic equations 


08/0C, = (QP/OC,)(d8/dP)¢, (3) 


oF (B) (Ss) 4) 
OC. C2\dP/o,  2C2\dP/o, 
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From Eqs. (1)—(4), 

ey eM a 
oC; dé /,\dP/¢, (<) 


OP a 
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(a) 
dé /.\dP/¢ 


c/a = sin (u + a — 8)/sin (u — a + 8) (6) 


(5) 
but from Fig. 1, 


where yu is the Mach angle behind the shock front, 
Hence, 


= (ea) 
OC, dé s dP C1 sin (u — a = 6) 


ee ~sin(v— ata) 
dC> dé ],\dP/¢, 


Now applying the second condition and proceeding asin 
Eqs. (2)—(4), we obtain 
) 


But from Fig. 1, 





dé = oP 
C1 


De — dp = —K,d = ‘oo a en oe 
oe ‘= = 


where K,, is the curvature of the wall. 


e/d = c/d = 1/2 cosu (9) 
Hence 
oP oP e 2K,, COs pu (10 
oC, 3G. —«(d8/dP)c, mM) 
and 
dP oP oP 
d=—e+ =e 


dw oo 
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Now concerning the shock curvature K,, 
dé 
da ds 
= <= —lUcU TC Ol (12) 
ds (=) 
da/, 
But from Eqs. (2)-(4), 
a ds, *) oP oF «| 
— = — = —— —a1 — 
= ae dP/¢\0C, 2G 
whence 
dé _ (=) sin (u — a+ 4) . 
ds \dP CO: sin 2u 
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(dé/da), 


sin 2u 
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oes . a ee aeeeremnonete a, Lae 
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Now Eggs. (7) and (10) determine 0P/O0C; and 0P/O0C2 
m terms of (OP/dé),, (dé/dP)¢,, wu, a, 6, and Ky. By 
specifying 6 and K,, (fixing the problem geometrically), 
4, a, and P can be readily obtained from the well- 
known solutions for the supersonic wedge flow as func- 
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tions of the Mach Number of the oncoming flow. The 
basic characteristic equation yields for (d5/dP)¢, 
(dé/dP)c, = sine2u/2yP (15) 


where y is the adiabatic exponent, and (dP/dé), can be 
found from differentiation of familiar relations for in- 
clined shock fronts 


(dP /d5), = (OP/0a);,/(06/0a), (16) 
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By use of these relations, then, from Eqs. (11) and (14), 
dP/dw and K, are determined as functions of M, 6, ang 
K,. The values of surface-pressure gradient and shoe} 
curvature thus computed may be considered as approx. 
imate values that approach the exact values as the 
So | régime considered shrinks to the edge of the ogive.* 

The results of these computations (for air y = 1.405 
are presented graphically in Figs. 2, 3, 4, and 5. 
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« approximation, the surface-pressure gradient at the 
ds " edge is given by a 
: (2) ee | 
¥ = — dw) ew + pm” dé Jo dw ie (d5/dP)c, oe 
‘ <i yet be = ge oon as where (dé/dP) ., is given by Eq. (15). By use of Eq. 
i, 3 4 (17) and the exact value of dP/dw computed in the 
Peer eae or previous section, a correction factor \ (a function of ¥ 
Fic. 6. and 6) can be computed such that ghich re 
J/dP\.- the attac 
where (dP/dw) exact = » (=) (18 Since 1 
dw) uw + pm | 
ae ~ sla . ; both gre< 
a ( “7 Mi)! sin? a — 1) Similarly, the surface-pressure gradient at the edge | to make 
Pies I given by the “‘tangent-wedge”’ approximation of Mach 
and , (=) a (=) a (2) as 
6 t 4 [2/(y — 1)] + Mi? sin? a \ dw) wp wr dS/,dw “\ dé /, : 
= a-— tan > ; rage 7 " —— 
ly + 1)/2(y — 1)]M;? sin 2a where (dP/d6), is computed from Eq. (16), can be used 
to compute a correction factor V (a function of 1M; and 


6) such that 
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The results of computation of \ and W are presented 
graphically in Figs. 6, 7, 8, and 9. d 
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It will be observed that limiting lines appear in the 
graphs of the results because steady flow with attached 
shock is possible only in a certain régime of combina 
tions of Mach Numbers and wedge angles. For any ur 
given value of half-wedge angle 6 (< sin~' 1/7) there 
exists a minimum Mach Number for which detachment 
of the shock does not occur. Also, for any given Mach 2 
Number, there exists a maximum half-wedge angle for 
which detachment of the shock does not occur. The 
maximum possible half-wedge angle for which flow with C 
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Y 4M, =|2 
Lo 3 attached shock is possible is 
5 i 
a 6 = sin~! 1/y (6 = 45.38° for y = 1.405) 
wa 5° 10° 15° 20° 25° 30° 35° * For a more general derivation of relations involving curved 
OGIVE ANGLE (8) shocks see ‘‘On Curved Shock Waves,” by T. Y. Thomas, Journal 
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which requires an infinite Mach Number to maintain approximation method should be safe to use. From 
the attachment. Figs. 6 and 8, however, it would appear that generally 
(18)] Since the correction factors \ and W take on values the “tangent-wedge’’ approximation is better at P 
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Announcing a New Department .. . eaders MUM | 
(20 oe : 
It has been suggested by several members of the I.A.S. Editorial Committee that the Journal open 
“in a special department to fill the need for the more immediate communication of essential results of research 
y ° A ° ° ° 
in the greatly expanded and accelerated fields of the aeronautical sciences. Therefore, it has been 
decided to establish a new department, the ‘Readers’ Forum,”” which will include letters to the editor 
for the quick publication of (1) brief reports of important findings in aeronautical sciences and (2) dis- 
cussions of papers published in the Journal. 
in the 
ached Specifications for Submission of Material 
ibina- (1) All letters should be accompanied by a suggested short title, giving some indication to the reader of the subject 
* any under discussion. 
there (2) These letters should not exceed 800 words in length. 
iment (3) Publication will be completed 6 to 8 weeks after receipt of the material. 
Mach - (4) The Editorial Committee will not hold itself responsible for the opinions expressed by the correspondents, 
: and no proofs will be sent to the authors. 
le 2 (5) All letters should be double-spaced. 
The 
7 with It is hoped that the new department may be started with the January, 1949, issue of the Journal 
Contributions should be sent to: 
) Editor, Journal of the Aeronautical Sciences 
2 East 64th street 
— New York 21, New York 
ourne: 
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